Page tree
Skip to end of metadata
Go to start of metadata

See: http://www.stat.berkeley.edu/users/terry/Classes/s246.2006/Week9/Wk9L2.pdf

Code:

 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)
}
  • No labels