# HG changeset patch # User davidvanzessen # Date 1386678176 18000 # Node ID a5c224bb0be5f18f4c492562308e8848611efa29 # Parent 8d83319a0f3dcd126210b2ec8fce681c38b537a4 Uploaded diff -r 8d83319a0f3d -r a5c224bb0be5 RScript.r --- a/RScript.r Tue Dec 10 05:53:08 2013 -0500 +++ b/RScript.r Tue Dec 10 07:22:56 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"])