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

Changed VCF subsetting procedure. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4742 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
kiran 2010-11-28 00:46:29 +00:00
parent f8a3bd7243
commit 9cca14acc5
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
R/VariantReport/VariantReport.R
View File

@ -249,8 +249,14 @@ plot.alleleCountSpectrum <- function(eval, novelty_name="all", col=c("#FF6342",
); );
} }
eval.getMetricsByAc <- function(eval, jexl_expression, novelty_name="all") { eval.getMetricsByAc <- function(eval, jexl, novelty="all") {
piece = eval$MetricsByAc[which(eval$MetricsByAc$evaluation_name == "eval" & eval$MetricsByAc$comparison_name == "dbsnp" & as.character(eval$MetricsByAc$jexl_expression) == as.character(jexl_expression) & eval$MetricsByAc$filter_name == "called" & eval$MetricsByAc$novelty_name == novelty_name),]; piece = subset(eval$MetricsByAc,
evaluation_name == "eval" &
comparison_name == "dbsnp" &
as.character(jexl_expression) == as.character(jexl) &
filter_name == "called" &
novelty_name == novelty
);
} }
plot.titvSpectrum <- function(eval, novelty_name="all", col=c("#FF6342", "#FF9675", "#5C92A4", "#88EEFF", "#55BBFF")) { plot.titvSpectrum <- function(eval, novelty_name="all", col=c("#FF6342", "#FF9675", "#5C92A4", "#88EEFF", "#55BBFF")) {