Displaying 1 result from an estimated 1 matches for "errcol".
Did you mean:
errcode
2011 May 19
3
A better way to do this
...min(dat2[,molecule],na.rm=T)+4,max(dat2[,molecule],na.rm=T)+1))
all.sub.id<-dat2$Subject.ID
if (group=='control'){
col=c('blue')
}else{col=c('red')}
if(group2=='case'){
col2=c('red')
}else{ col2=c('blue')}
uniq.sub.id<-unique(all.sub.id)
errcol<-c()
for (i in 1: length(uniq.sub.id))
{
sub<-dat2[grep(uniq.sub.id[i],dat2$Subject.ID),]
sub<- sub[order(sub$Age.at.Sample.Collection),]
sub<-sub[sub[,molecule]>cutoff,]
sub.id<-uniq.sub.id[i]
if (dim(sub)[1]<=1) errcol<-c(errcol, sub.id)...