Median polish to consolidate quantitations

 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)
}