X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=rsem-plot-model;fp=rsem-plot-model;h=c15a75f230c018da6ffc0cde6171f6d1f6fa0686;hb=de840718636408cffb3ac1b614bb26a5ad5dec91;hp=7466eea2398be524e3f073b1ac7d914b3c6a8d83;hpb=0f1ac7d42823cf0f47d25f4f0a04620357421589;p=rsem.git diff --git a/rsem-plot-model b/rsem-plot-model index 7466eea..c15a75f 100755 --- a/rsem-plot-model +++ b/rsem-plot-model @@ -82,10 +82,10 @@ if (model_type == 1 || model_type == 3) { if (sum(c(vecA, vecC, vecG, vecT)) < 1e-8) next x <- c(x, (i - 1)) - peA <- c(peA, ifelse(sum(vecA) < 1e-8, NA, -10 * log(1.0 - vecA[1]))) - peC <- c(peC, ifelse(sum(vecC) < 1e-8, NA, -10 * log(1.0 - vecC[2]))) - peG <- c(peG, ifelse(sum(vecG) < 1e-8, NA, -10 * log(1.0 - vecG[3]))) - peT <- c(peT, ifelse(sum(vecT) < 1e-8, NA, -10 * log(1.0 - vecT[4]))) + peA <- c(peA, ifelse(sum(vecA) < 1e-8, NA, -10 * log10(1.0 - vecA[1]))) + peC <- c(peC, ifelse(sum(vecC) < 1e-8, NA, -10 * log10(1.0 - vecC[2]))) + peG <- c(peG, ifelse(sum(vecG) < 1e-8, NA, -10 * log10(1.0 - vecG[3]))) + peT <- c(peT, ifelse(sum(vecT) < 1e-8, NA, -10 * log10(1.0 - vecT[4]))) } matplot(x, cbind(peA, peC, peG, peT), type = "b", lty = 1:4, pch = 0:3, col = 1:4, main = "Phred Quality Score vs. Observed Quality", xlab = "Quality Score", ylab = "Observed Quality")