Mercurial > repos > davidvanzessen > argalaxy_tools

comparison report_clonality/RScript.r @ 39:2e0a7c35082e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 22 Jan 2016 10:54:01 -0500
parents fd55088e48d2
children d97e1421aa86
comparison
equal deleted inserted replaced
38:539618b939bc 39:2e0a7c35082e
562 clonalFreqCount$WeightedCount = clonalFreqCount$Count * clonalFreqCount$Weight 562 clonalFreqCount$WeightedCount = clonalFreqCount$Count * clonalFreqCount$Weight
563 563
564 ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "clonaltype")]) 564 ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "clonaltype")])
565 ReplicateReads = data.frame(data.table(ReplicateReads)[, list(Reads=.N), by=c("Sample", "Replicate")]) 565 ReplicateReads = data.frame(data.table(ReplicateReads)[, list(Reads=.N), by=c("Sample", "Replicate")])
566 clonalFreqCount$Reads = as.numeric(clonalFreqCount$Reads) 566 clonalFreqCount$Reads = as.numeric(clonalFreqCount$Reads)
567 ReplicateReads$squared = ReplicateReads$Reads * ReplicateReads$Reads 567 ReplicateReads$Reads = as.numeric(ReplicateReads$Reads)
568 ReplicateReads$squared = as.numeric(ReplicateReads$Reads * ReplicateReads$Reads)
568 569
569 ReplicatePrint <- function(dat){ 570 ReplicatePrint <- function(dat){
570 write.table(dat[-1], paste("ReplicateReads_", unique(dat[1])[1,1] , ".csv", sep=""), sep=",",quote=F,na="-",row.names=F,col.names=F) 571 write.table(dat[-1], paste("ReplicateReads_", unique(dat[1])[1,1] , ".csv", sep=""), sep=",",quote=F,na="-",row.names=F,col.names=F)
571 } 572 }
572 573