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:
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:j...@joepking.com> j...@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@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.