# HG changeset patch # User davidvanzessen # Date 1386678145 18000 # Node ID 3ce07f5889ad2bef34d6b7d4fc2a6016dad1f4df # Parent 8ba0fd5b03a17e913ab2f7bd325b836a89930dd1 Uploaded diff -r 8ba0fd5b03a1 -r 3ce07f5889ad RScript.r --- a/RScript.r Tue Dec 10 05:52:37 2013 -0500 +++ b/RScript.r Tue Dec 10 07:22:25 2013 -0500 @@ -93,7 +93,7 @@ setwd(outDir) -write.table(PRODF, "allUnique.tsv", sep="\t",quote=F,row.names=T,col.names=T) +write.table(PRODF, "allUnique.tsv", sep="\t",quote=F,row.names=F,col.names=T) pV = ggplot(PRODFV) pV = pV + geom_bar( aes( x=factor(reorder(Top.V.Gene, chr.orderV)), y=relFreq, fill=Sample), stat='identity', position="dodge") + theme(axis.text.x = element_text(angle = 90, hjust = 1)) @@ -236,17 +236,17 @@ clonalityFrame = PROD clonalityFrame$ReplicateConcat = do.call(paste, c(clonalityFrame[c("VDJCDR3", "Sample", "Replicate")], sep = ":")) clonalityFrame = clonalityFrame[!duplicated(clonalityFrame$ReplicateConcat), ] - write.table(clonalityFrame, "clonalityComplete.tsv", sep="\t",quote=F,row.names=T,col.names=T) + write.table(clonalityFrame, "clonalityComplete.tsv", sep="\t",quote=F,row.names=F,col.names=T) ClonalitySampleReplicatePrint <- function(dat){ - write.table(dat, paste("clonality_", unique(dat$Sample) , "_", unique(dat$Replicate), ".tsv", sep=""), sep="\t",quote=F,row.names=T,col.names=T) + write.table(dat, paste("clonality_", unique(dat$Sample) , "_", unique(dat$Replicate), ".tsv", sep=""), sep="\t",quote=F,row.names=F,col.names=T) } clonalityFrameSplit = split(clonalityFrame, f=clonalityFrame[,c("Sample", "Replicate")]) lapply(clonalityFrameSplit, FUN=ClonalitySampleReplicatePrint) ClonalitySamplePrint <- function(dat){ - write.table(dat, paste("clonality_", unique(dat$Sample) , ".tsv", sep=""), sep="\t",quote=F,row.names=T,col.names=T) + write.table(dat, paste("clonality_", unique(dat$Sample) , ".tsv", sep=""), sep="\t",quote=F,row.names=F,col.names=T) } clonalityFrameSplit = split(clonalityFrame, f=clonalityFrame[,"Sample"])