v1.1.5.b

[rsem.git] / rsem-plot-model

diff --git a/rsem-plot-model b/rsem-plot-model
index c0eaa1b0b70faa4175478830b8e5d26fd4eff61c..333a9e7b2e9b2ba1371f06f56e52d3b6cbbf6a87 100755 (executable)
--- a/rsem-plot-model
+++ b/rsem-plot-model
@@ -96,10 +96,10 @@ if (model_type == 1 || model_type == 3) {
     vecG <- as.numeric(list[[3]])
     vecT <- as.numeric(list[[4]])
     if (sum(c(vecA, vecC, vecG, vecT)) < 1e-8) break
-    peA <- c(peA, ifelse(sum(vec) < 1e-8, NA, 1.0 - vecA[1]))
-    peC <- c(peC, ifelse(sum(vec) < 1e-8, NA, 1.0 - vecC[2]))
-    peG <- c(peG, ifelse(sum(vec) < 1e-8, NA, 1.0 - vecG[3]))
-    peT <- c(peT, ifelse(sum(vec) < 1e-8, NA, 1.0 - vecT[4]))
+    peA <- c(peA, ifelse(sum(vec) < 1e-8, NA, (1.0 - vecA[1]) * 100))
+    peC <- c(peC, ifelse(sum(vec) < 1e-8, NA, (1.0 - vecC[2]) * 100))
+    peG <- c(peG, ifelse(sum(vec) < 1e-8, NA, (1.0 - vecG[3]) * 100))
+    peT <- c(peT, ifelse(sum(vec) < 1e-8, NA, (1.0 - vecT[4]) * 100))
   }
 
   x <- 1 : length(peA)