Question: qQCReport Missing Plots Described in Documentation
gravatar for Dario Strbenac
7 weeks ago by
Dario Strbenac1.4k
Dario Strbenac1.4k wrote:

According to the documentation, there are eight kinds of plots generated by qQCReport. However, I only see four.

testFile <- system.file("extdata", "ex1.bam", package="Rsamtools")
qQCReport(testFile, "testQC.pdf")

Can the function's result be changed to match the documentation? Also, each plot is very small and takes up a small part of a large A4 PDF page, which affects plot legibility. QuasR version 1.22.1 is used.

quasr • 89 views
ADD COMMENTlink modified 6 weeks ago by Hotz, Hans-Rudolf400 • written 7 weeks ago by Dario Strbenac1.4k
Answer: qQCReport Missing Plots Described in Documentation
gravatar for Hotz, Hans-Rudolf
6 weeks ago by
Hotz, Hans-Rudolf400 wrote:

Hi Dario

The number of plots you get depends on the input, e.g. if you run qQCReport on a qproject generated from paired-end fastq files, you will get 8 plots.

The plot sizes have been optimised to display several samples per page.

Hope this helps/clarifies Regards, Hans-Rudolf

ADD COMMENTlink written 6 weeks ago by Hotz, Hans-Rudolf400

But the documentation ?qQCReport provides a list of seven plots which may be generated for any BAM file and an eighth plot which is for BAM files containing paired reads. The test file I used has paired reads

samtools view -c -f 1 Rsamtools/extdata/ex1.bam

So, I expect to get eight plots but I only get four. Could you run the test code to check how many plots you obtain? The documentation of qQCReport and its output seems inconsistent.

ADD REPLYlink written 6 weeks ago by Dario Strbenac1.4k

Yes, I agree the documentation is misleading (or even partly wrong). I will talk to the QuasR developers and try to get the documentation improved for of qQCReport.


ADD REPLYlink modified 6 weeks ago • written 6 weeks ago by Hotz, Hans-Rudolf400

Hi Dario

The documentation and also plotting issues have been addressed in QuasR 1.23.17 (will become available in BioC-devel in a day or so).

The documentation has been simplified to make it clear that some plots will not be available if the necessary information is missing (e.g. base quality plots when working with a fasta file). Furthermore, qQCReport now has a new argument a4paper (a4paper = TRUE by default) that allows you to switch off the page-optimised layout for mapping and uniqueness plots.

Thank you for your report.


ADD REPLYlink written 4 weeks ago by Michael Stadler330
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 246 users visited in the last hour