Displaying 1 result from an estimated 1 matches for "log2_".
Did you mean:
log2
2011 May 19
3
A better way to do this
...b[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)
if (dim(sub)[1]>1)
{
par(new=TRUE)
plot(log2(sub[,molecule])~sub$Age.at.Sample.Collection,
ylab=paste('Log2_',molecule,sep=''),xlab="Age at Sample",pch=1, ylim=ylim,
xlim=c(0,25),main=paste(group,'-',molecule))
mod<-loess(log2(sub[,molecule])~Age.at.Sample.Collection,
na.action=na.exclude, data=sub)
pred<-predict(mod)
lines(pred~Age.at.Sample.Collectio...