See: http://www.stat.berkeley.edu/users/terry/Classes/s246.2006/Week9/Wk9L2.pdf
Code:
| Code Block |
|---|
finalans<-data.frame()
for (subs in unique(final$V13)) {
m<-medpolish(final[which(final$V13==subs),c(3,7)],trace.iter=FALSE);
msd<-medpolish(final[which(final$V13==subs),c(4,8)],trace.iter=FALSE);
myadd<-data.frame(gene=I(subs),cntmean=m$col[1]+m$overall,mutmean=m$col[2]+m$overall,cntsd=msd$col[1]+msd$overall, mutsd=msd$col[2]+msd$overall)
finalans<-rbind(finalans,myadd)
}
|