Off topic:ChIPQC without peak data fails to generate a report (due to invocation of plotCorHeatmap)
2
0
Entering edit mode
Coby Viner ▴ 50
@coby-viner-10939
Last seen 3.0 years ago
University of Toronto, Canada

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

 

ChIPQC chipqc • 766 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 908 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6