Entering edit mode
While generation of quality metrics using ChIPQC for files without peak data has been resolved (ChIPQC without peak calls), I am unable to save HTML reports for those ChIPQC runs.
ChIPQCreport
fails, apparently due to attempting to invoke plotCorHeatmap
, despite the absence of peaks:
> experiment Samples: 2 : ENCFF863PSQ ENCFF631ENA Tissue Factor Control Replicate ENCFF863PSQ DOHH2 CTCF ENCFF631ENA 1 Reads Map% Filt% Dup% ReadL FragL RelCC SSD RiP% ENCFF863PSQ 12409070 100 0 0 36 252 5.230 0.176 NA ENCFF631ENA 23751742 100 0 0 36 229 0.973 0.112 NA > ChIPQCreport(experiment, facet=F) Saving 7 x 7 in image Saving 7 x 7 in image Error in plotCorHeatmap(object, attributes = c(facetBy, colourBy)) : No peaks to plot. In addition: Warning messages: 1: In grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : X11 used font size 34 when 45 was requested 2: In grid.Call.graphics(L_rect, x$x, x$y, x$width, x$height, resolveHJust(x$just, : semi-transparency is not supported on this device: reported only once per page 3: In grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : X11 used font size 34 when 45 was requested 4: In grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : X11 used font size 34 when 83 was requested 5: In grid.Call(L_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : X11 used font size 34 when 66 was requested > traceback() 5: stop("No peaks to plot.") 4: plotCorHeatmap(object, attributes = c(facetBy, colourBy)) 3: plotCorHeatmap(object, attributes = c(facetBy, colourBy)) 2: ChIPQCreport(experiment, facet = F) 1: ChIPQCreport(experiment, facet = F) > sessionInfo() R version 3.3.0 (2016-05-03) Platform: x86_64-pc-linux-gnu (64-bit) Running under: CentOS release 6.5 (Final) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats4 parallel stats graphics grDevices utils datasets [8] methods base other attached packages: [1] ChIPQC_1.8.7 DiffBind_2.0.2 [3] SummarizedExperiment_1.2.3 Biobase_2.32.0 [5] GenomicRanges_1.24.2 GenomeInfoDb_1.8.1 [7] IRanges_2.6.1 S4Vectors_0.10.1 [9] BiocGenerics_0.18.0 ggplot2_2.1.0 loaded via a namespace (and not attached): [1] edgeR_3.14.0 [2] splines_3.3.0 [3] TxDb.Hsapiens.UCSC.hg18.knownGene_3.2.2 [4] TxDb.Mmusculus.UCSC.mm10.knownGene_3.2.2 [5] gtools_3.5.0 [6] assertthat_0.1 [7] latticeExtra_0.6-28 [8] amap_0.8-14 [9] RBGL_1.48.1 [10] Rsamtools_1.24.0 [11] Category_2.38.0 [12] RSQLite_1.0.0 [13] backports_1.0.2 [14] lattice_0.20-33 [15] limma_3.28.10 [16] digest_0.6.10 [17] RColorBrewer_1.1-2 [18] XVector_0.12.0 [19] checkmate_1.8.0 [20] colorspace_1.2-6 [21] Matrix_1.2-6 [22] plyr_1.8.4 [23] GSEABase_1.34.0 [24] chipseq_1.22.0 [25] TxDb.Hsapiens.UCSC.hg38.knownGene_3.1.3 [26] XML_3.98-1.4 [27] pheatmap_1.0.8 [28] ShortRead_1.30.0 [29] biomaRt_2.28.0 [30] genefilter_1.54.2 [31] zlibbioc_1.18.0 [32] xtable_1.8-2 [33] GO.db_3.3.0 [34] scales_0.4.0 [35] brew_1.0-6 [36] gdata_2.17.0 [37] TxDb.Rnorvegicus.UCSC.rn4.ensGene_3.2.2 [38] BiocParallel_1.6.2 [39] annotate_1.50.0 [40] GenomicFeatures_1.24.2 [41] survival_2.39-4 [42] magrittr_1.5 [43] systemPipeR_1.6.2 [44] fail_1.3 [45] gplots_3.0.1 [46] hwriter_1.3.2 [47] GOstats_2.38.0 [48] graph_1.50.0 [49] tools_3.3.0 [50] BBmisc_1.9 [51] stringr_1.1.0 [52] sendmailR_1.2-1 [53] munsell_0.4.3 [54] locfit_1.5-9.1 [55] AnnotationDbi_1.34.3 [56] Biostrings_2.40.2 [57] caTools_1.17.1 [58] grid_3.3.0 [59] RCurl_1.95-4.8 [60] TxDb.Celegans.UCSC.ce6.ensGene_3.2.2 [61] rjson_0.2.15 [62] AnnotationForge_1.14.2 [63] labeling_0.3 [64] bitops_1.0-6 [65] base64enc_0.1-3 [66] gtable_0.2.0 [67] DBI_0.4-1 [68] reshape2_1.4.1 [69] R6_2.1.3 [70] GenomicAlignments_1.8.3 [71] Nozzle.R1_1.1-1 [72] dplyr_0.4.3 [73] rtracklayer_1.32.1 [74] KernSmooth_2.23-15 [75] stringi_1.1.1 [76] TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2 [77] BatchJobs_1.6 [78] Rcpp_0.12.7 [79] TxDb.Dmelanogaster.UCSC.dm3.ensGene_3.2.2 [80] TxDb.Mmusculus.UCSC.mm9.knownGene_3.2.2