From: Bo Li Date: Mon, 18 Apr 2011 14:48:10 +0000 (-0500) Subject: RSEM v1.1.8 X-Git-Url: https://git.donarmstrong.com/?p=rsem.git;a=commitdiff_plain;h=8184b0257104954c1da1d42484f5bf2d2c308c5a RSEM v1.1.8 --- diff --git a/README.md b/README.md index 2bb6015..f8c4e27 100644 --- a/README.md +++ b/README.md @@ -148,7 +148,8 @@ The plots generated depends on read type and user configuration. It may include fragment length distribution, mate length distribution, read start position distribution (RSPD), quality score vs observed quality given a reference base, position vs percentage of sequencing -error given a reference base and histogram of read alignments. +error given a reference base and histogram of reads with different +number of alignments. fragment length distribution and mate length distribution: x-axis is fragment/mate length, y axis is the probability of generating a fragment/mate with the associated length @@ -158,7 +159,7 @@ Quality score vs. observed quality given a reference base: x-axis is Phred quali Position vs. percentage sequencing error given a reference base: x-axis is position and y-axis is percentage sequencing error -Histogram of read alignments: x-axis is the number of alignments a read has and y-axis is the number of such reads +Histogram of reads with different number of alignments: x-axis is the number of alignments a read has and y-axis is the number of such reads. The inf in x-axis means number of reads filtered due to too many alignments ## Example diff --git a/rsem-plot-model b/rsem-plot-model index 8ff626a..7466eea 100755 --- a/rsem-plot-model +++ b/rsem-plot-model @@ -124,6 +124,6 @@ if (model_type == 1 || model_type == 3) { close(con) pair <- read.table(file = cntF, skip = 3, sep = "\t") -barplot(pair[,2], names.arg = pair[,1], xlab = "Number of Alignments", ylab = "Number of Reads", main = "Histogram of Read Alignments") +barplot(pair[,2], names.arg = pair[,1], xlab = "Number of Alignments", ylab = "Number of Reads", main = "Histogram of Reads with Different Number of Alignments") dev.off()