X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=rsem-plot-model;fp=rsem-plot-model;h=333a9e7b2e9b2ba1371f06f56e52d3b6cbbf6a87;hb=e123915d329682a1712bf77204e556baadb69a2f;hp=c0eaa1b0b70faa4175478830b8e5d26fd4eff61c;hpb=a42d21510e502fe5cdd21a5d7ad023b33d2e6b99;p=rsem.git

diff --git a/rsem-plot-model b/rsem-plot-model
index c0eaa1b..333a9e7 100755
--- 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)