From 94cc97d0c906a42fc954d1f238ff699705f1a1fb Mon Sep 17 00:00:00 2001 From: Bo Li Date: Fri, 18 Mar 2011 15:37:35 -0500 Subject: [PATCH] a bug in rsem-plot-model is fixed --- rsem-plot-model | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) 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]]) -- 2.39.2