Mercurial > repos > davidvanzessen > argalaxy_tools
changeset 39:2e0a7c35082e draft
Uploaded
author | davidvanzessen |
---|---|
date | Fri, 22 Jan 2016 10:54:01 -0500 |
parents | 539618b939bc |
children | a2fdcf1a2e50 |
files | report_clonality/RScript.r |
diffstat | 1 files changed, 2 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/report_clonality/RScript.r Fri Jan 22 10:44:50 2016 -0500 +++ b/report_clonality/RScript.r Fri Jan 22 10:54:01 2016 -0500 @@ -564,7 +564,8 @@ ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "clonaltype")]) ReplicateReads = data.frame(data.table(ReplicateReads)[, list(Reads=.N), by=c("Sample", "Replicate")]) clonalFreqCount$Reads = as.numeric(clonalFreqCount$Reads) - ReplicateReads$squared = ReplicateReads$Reads * ReplicateReads$Reads + ReplicateReads$Reads = as.numeric(ReplicateReads$Reads) + ReplicateReads$squared = as.numeric(ReplicateReads$Reads * ReplicateReads$Reads) ReplicatePrint <- function(dat){ write.table(dat[-1], paste("ReplicateReads_", unique(dat[1])[1,1] , ".csv", sep=""), sep=",",quote=F,na="-",row.names=F,col.names=F)