a bug in rsem-plot-model is fixed

diff --git a/rsem-plot-model b/rsem-plot-model
index 333a9e7b2e9b2ba1371f06f56e52d3b6cbbf6a87..e9225eddf41466424d265cc9b55ef7ebf96b7aa4 100755 (executable)
--- 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]])