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

diff --git a/rsem-plot-model b/rsem-plot-model
index 333a9e7..e9225ed 100755
--- a/rsem-plot-model
+++ b/rsem-plot-model
@@ -23,7 +23,8 @@ std <- sqrt(weighted.mean((x - mean)^2, y))
 plot(x, y, type = "h", main = "Fragment Length Distribution", sub = paste("Mean = ", mean, ", Std = ", std), xlab = "Fragment Length", ylab = "Probability")  
 
 # mate length distribution
-bval <- as.numeric(readLines(con, n = 1)[1])
+if (model_type == 0 || model_type == 1) bval <- as.numeric(readLines(con, n = 1)[1]) else bval <- 1
+
 if (bval == 1) {
   list <- strsplit(readLines(con, n = 2), split = " ")
   vec <- as.numeric(list[[1]])