search for: probes_blood_sub

...This can solve my problem that I need to set the NA values manually and reset j in for-loop every time error returned. Doing this manually would be very tedious. m1<-matrix(ncol=3,nrow=0) for (j in 1:1000){ fit<-lme(Expression~1,random=~1|MonkeyID, data=mydata1,subset=Probes==as.character(Probes_Blood_sub[j,"probe"])) ci<-intervals(fit)$reStruct$MonkeyID^2 m1<-rbind(m1,ci)} Thanks,Sityee _________________________________________________________________ The i’m Talkathon starts 6/24/08.  For now, give amongst yourselves. [[alternative HTML version deleted]]