From 9cca14acc5a7fae25f21657406f9e9e2e2f78463 Mon Sep 17 00:00:00 2001
From: kiran <kiran@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Sun, 28 Nov 2010 00:46:29 +0000
Subject: [PATCH] Changed VCF subsetting procedure.

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4742 348d0f76-0448-11de-a6fe-93d51630548a
---
 R/VariantReport/VariantReport.R | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/R/VariantReport/VariantReport.R b/R/VariantReport/VariantReport.R
index 525b7ea4d..7f83eacfd 100644
--- a/R/VariantReport/VariantReport.R
+++ b/R/VariantReport/VariantReport.R
@@ -249,8 +249,14 @@ plot.alleleCountSpectrum <- function(eval, novelty_name="all", col=c("#FF6342",
     );
 }
 
-eval.getMetricsByAc <- function(eval, jexl_expression, novelty_name="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),]; 
+eval.getMetricsByAc <- function(eval, jexl, novelty="all") {
+    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")) {