Re: [R] using latticeExtra plotting confidence intervals

From: Joe P King <jp_at_joepking.com>
Date: Sun, 13 Jun 2010 22:23:13 -0700

So if the categories was race (A and B) and I had a male and female nested in each group, that would give me 4 different data points but I want two data points within each panel.

Joe King
206-913-2912
jp_at_joepking.com
"Never throughout history has a man who lived a life of ease left a name worth remembering." --Theodore Roosevelt

-----Original Message-----
From: Deepayan Sarkar [mailto:deepayan.sarkar_at_gmail.com] Sent: Sunday, June 13, 2010 10:15 PM
To: Joe P King
Cc: r-help_at_r-project.org
Subject: Re: [R] using latticeExtra plotting confidence intervals

On Sun, Jun 13, 2010 at 10:10 AM, Joe P King <jp_at_joepking.com> wrote:
> I am wanting to plot a 95% confidence band using segplot, yet I am wanting
> to have groups. For example if I have males and females, and then I have
> them in different races, I want the racial groups in different panels. I
> have this minor code, completely made up but gets at what I am wanting, 4
> random samples and 4 samples of confidence, I know how to get A & B into one
> panel and C&D in to another but how do I get the x axis to label them
> properly and have it categorized as two. I am not sure what to have to the
> left side of the formula. This is the example code:

(1) Your code results in

> length(categories)

[1] 2
> length(mu)

[1] 4

which makes the formula mu~ci.upper+ci.lower|categories meaningless.

(2) You are effectively plotting xyplot(mu ~ mu | categories), with additional confidence intervals in one direction. I'm sure that's not what you want, but it's not clear what it is that you do actually want.

-Deepayan

>
>
>
> library(lattice)
>
> library(latticeExtra)
>
>
>
> sample1<-rnorm(100,10,2)
>
> sample2<-rnorm(100,50,3)
>
> sample3<-rnorm(100,20,2)
>
> sample4<-rnorm(100,40,1)
>
> mu1<-mean(sample1)
>
> ci.upper1<-mu1+2*2
>
> ci.lower1<-mu1-2*2
>
> mu2<-mean(sample2)
>
> ci.upper2<-mu2+2*3
>
> ci.lower2<-mu2-2*3
>
> mu3<-mean(sample3)
>
> ci.upper3<-mu3+2*2
>
> ci.lower3<-mu3-2*2
>
> mu4<-mean(sample4)
>
> ci.upper4<-mu4+2*1
>
> ci.lower4<-mu4-2*1
>
> categories<-c("A","B")
>
>
>
> mu<-cbind(mu1,mu2,mu3,mu4)
>
> ci.upper<-cbind(ci.upper1,ci.upper2,ci.upper3,ci.upper4)
>
> ci.lower<-cbind(ci.lower1,ci.lower2,ci.lower3,ci.lower4)
>
> segplot(mu~ci.upper+ci.lower|categories, centers = mu, horizontal=FALSE)
>
>
>
> I also tried this
>
>
>
> seq1<-seq(1,4,1)
>
> segplot(seq1~ci.upper+ci.lower|categories, centers = mu,horizontal=FALSE)
>
>
>
> but it also gives poor x axis, I know this is probably an elementary problem
> but any help would be greatly appreciated.
>
>
>
> Heres my data structure, sorry for bombarding with all the code.
>
>
>
> structure(c(9.85647167881417, 50.1856561919426, 19.8477661576365,
>
> 39.8575819498655, 13.8564716788142, 56.1856561919426, 23.8477661576365,
>
> 41.8575819498655, 5.85647167881417, 44.1856561919426, 15.8477661576365,
>
> 37.8575819498655), .Dim = c(1L, 12L), .Dimnames = list(NULL,
>
> c("mu1", "mu2", "mu3", "mu4", "ci.upper1", "ci.upper2", "ci.upper3",
>
> "ci.upper4", "ci.lower1", "ci.lower2", "ci.lower3", "ci.lower4"
>
> )))
>
> -------------------------------------------
>
> Joe King, M.A.
>
> Ph.D. Student
>
> University of Washington - Seattle
>
> Office: 404A
>
> Miller Hall
>
> 206-913-2912
>
> <mailto:jp_at_joepking.com> jp_at_joepking.com
>
> -------------------------------------------
>
> "Never throughout history has a man who lived a life of ease left a name
> worth remembering." --Theodore Roosevelt
>
>
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon 14 Jun 2010 - 08:31:01 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 14 Jun 2010 - 09:10:29 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive