From: Bo Li <bli@cs.wisc.edu>
Date: Tue, 26 Apr 2011 03:17:02 +0000 (-0500)
Subject: fixed a bug in rsem-plot-model
X-Git-Url: https://git.donarmstrong.com/?p=rsem.git;a=commitdiff_plain;h=de840718636408cffb3ac1b614bb26a5ad5dec91

fixed a bug in rsem-plot-model
---

diff --git a/rsem-plot-model b/rsem-plot-model
index 7466eea..c15a75f 100755
--- a/rsem-plot-model
+++ b/rsem-plot-model
@@ -82,10 +82,10 @@ if (model_type == 1 || model_type == 3) {
 
     if (sum(c(vecA, vecC, vecG, vecT)) < 1e-8) next
     x <- c(x, (i - 1))
-    peA <- c(peA, ifelse(sum(vecA) < 1e-8, NA, -10 * log(1.0 - vecA[1])))
-    peC <- c(peC, ifelse(sum(vecC) < 1e-8, NA, -10 * log(1.0 - vecC[2])))
-    peG <- c(peG, ifelse(sum(vecG) < 1e-8, NA, -10 * log(1.0 - vecG[3])))
-    peT <- c(peT, ifelse(sum(vecT) < 1e-8, NA, -10 * log(1.0 - vecT[4])))
+    peA <- c(peA, ifelse(sum(vecA) < 1e-8, NA, -10 * log10(1.0 - vecA[1])))
+    peC <- c(peC, ifelse(sum(vecC) < 1e-8, NA, -10 * log10(1.0 - vecC[2])))
+    peG <- c(peG, ifelse(sum(vecG) < 1e-8, NA, -10 * log10(1.0 - vecG[3])))
+    peT <- c(peT, ifelse(sum(vecT) < 1e-8, NA, -10 * log10(1.0 - vecT[4])))
   }
 
   matplot(x, cbind(peA, peC, peG, peT), type = "b", lty = 1:4, pch = 0:3, col = 1:4, main = "Phred Quality Score vs. Observed Quality", xlab = "Quality Score", ylab = "Observed Quality")