Question: qQCReport Missing Plots Described in Documentation
0
gravatar for Dario Strbenac
7 months ago by
Dario Strbenac1.5k
Australia
Dario Strbenac1.5k wrote:

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

library(QuasR)
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 • 161 views
ADD COMMENTlink modified 7 months ago by Hotz, Hans-Rudolf400 • written 7 months ago by Dario Strbenac1.5k
Answer: qQCReport Missing Plots Described in Documentation
0
gravatar for Hotz, Hans-Rudolf
7 months ago by
Switzerland
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 7 months 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
3307

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 7 months ago by Dario Strbenac1.5k

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.

Hans-Rudolf

ADD REPLYlink modified 7 months ago • written 7 months ago by Hotz, Hans-Rudolf400
1

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.

Michael

ADD REPLYlink written 7 months ago by Michael Stadler330
Please log in to add an answer.

Help
Access

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