zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Updated Tearsheet with by sample QC metrics, bugfix for misnamed variables 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@5754 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
corin 2011-05-04 16:57:48 +00:00
parent 72a07e4553
commit 23efd66d31
1 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
R/DataProcessingReport/Tearsheet.R
View File

@ -1,7 +1,6 @@
#tearsheet (for cori's use) #tearsheet (for cori's use)
#New tearsheet generator #New tearsheet generator
.libPaths('/humgen/gsa-firehose2/pipeline/repositories/StingProduction/R/') #uncomment .libPaths('/humgen/gsa-firehose2/pipeline/repositories/StingProduction/R/')
#.libPaths('~/Documents/Sting/R/')
suppressMessages(library(gplots)); suppressMessages(library(gplots));
suppressMessages(library(ReadImages)); suppressMessages(library(ReadImages));
@ -19,7 +18,7 @@ cmdargs = gsa.getargs(
doc="Creates a tearsheet" doc="Creates a tearsheet"
); );
read.delim(tsv)->settable read.delim(cmdargs$tsv)->settable
squids<-unique(settable[,1]) squids<-unique(settable[,1])
@ -53,14 +52,13 @@ tearsheet<-function(){
#plot title bar #plot title bar
par(mar=c(0,0,0,0)) par(mar=c(0,0,0,0))
plot(drop) plot(drop)
text(155, 50, title, family="serif", adj=c(0,0), cex=3, col=gray(.25)) text(155, 50, cmdargs$title, family="serif", adj=c(0,0), cex=3, col=gray(.25))
# Project summary # Project summary
projects = paste(squids, collapse=", "); projects = paste(squids, collapse=", ");
# used_samples = nrow(tsv); used_samples = nrow(settable);
used_samples=6 #comment
unused_samples = 0; unused_samples = 0;