Dear all I would like to ask you about possible bug in using bwplot (from lattice) together with Hmisc library attached. I found it in my actual data, but here is a toy example. It appears only when some levels are missing.
library(lattice) library(Hmisc) # preparing data x1<-rnorm(10,5,1) x2<-rnorm(10,5,5) x3<-rnorm(10,1,1) x4<-rnorm(10,1,5) x<-c(x1,x2,x3,x4) x<-c(x,x+5) f2<-as.factor(rep(c("a","b","a","b"),c(20,20,20,20))) f3<-as.factor(rep(rep(c("c","d"),each=10),4)) f1<-as.factor(rep(c(1,2),c(40,40))) data<-data.frame(x,f1,f2,f3) rm(x,f1,f2,f3) attach(data) # 1st plot bwplot(f1~x|f2) # 2 is bigger than 1 **OK** # 2nd plot bwplot(f1~x|f2*f3) # again 2 is bigger than 1 **OK** # removing some levels detach("data") data<-data[-c(1:10,21:30),] attach(data) # 3rd plot bwplot(f1~x|f2) # 2 is bigger than 1 **OK** # 4th plot bwplot(f1~x|f2*f3) ### OOOPS 1 is bigger than 2 ***FALSE*** after detaching Hmisc detach("package:Hmisc") #5th plot bwplot(f1~x|f2*f3) # again 2 is bigger than 1 **OK** The levels for factor f1 are in the OOOPS case reversed, maybe it has something to do with factor redefinition as stated when attaching Hmisc library. Is it a bug or am I doing something wrong? R 1.7.1 WNT Package: Hmisc Version: 1.6-1 Date: 2003-06-21 Package: lattice Version: 0.7-14 Date: 2003/06/07 Thank you for any response. Petr Pikal [EMAIL PROTECTED] [EMAIL PROTECTED] ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help