Demo entry 6661438

facets

   

Submitted by anonymous on Nov 22, 2017 at 03:31
Language: RConsole. Code size: 379 Bytes.

args <- commandArgs(T)
library(facets)
rcmat = readSnpMatrix(args[2])
set.seed(args[1])
xx = preProcSample(rcmat)
oo=procSample(xx,cval=150)
fit=emcncf(oo)
x<-data.frame(purity=fit$purity,ploidy=fit$ploidy)
write.table(x,args[3],row.names=F,quote=F,sep="\t")
write.table(fit$cncf,args[4],row.names=F,quote=F,sep="\t")
pdf(args[5])
plotSample(x=oo,emfit=fit)
dev.off()

This snippet took 0.00 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).