ChIPQC failed on chrM
1
0
Entering edit mode
@seasky002002-13307
Last seen 4.1 years ago

Hi,

I encountered an error while using ChIPQC. It reported an error after calculation shift for chrM. I pasted the error message below. It looks like happened on chrM, so i change the parameter chromosomes to exclude chrM, and the ChIPQC program run successfully. 

But i still curious why the error happened on chrM. Would someone can give me an answer?

Thank you!

 

 

The error message:

Calculating shift for chrM

 33 / 300
 65 / 300
 97 / 300
 126 / 300
 155 / 300
 184 / 300
 213 / 300
 243 / 300
 272 / 300
 300 / 300
Error in names(res) <- nms :
  'names' attribute [18] must be the same length as the vector [1]
Calls: ChIPQC -> bplapply -> bplapply -> bplapply -> bplapply

 

ChIPQC • 1.0k views
ADD COMMENT
0
Entering edit mode

Hello

I get this same error at the start of my run. Was there a solution to this? My sessionInfo is attached below but also I get a message when I load the library, not sure if this has anything to do with it: No methods found in "RSQLite" for requests: dbGetQuery

Bam file has 93 contigs
Error in names(res) <- nms : 
  'names' attribute [28] must be the same length as the vector [4]
Calls: ChIPQC -> bplapply -> bplapply -> bplapply -> bplapply

 

> sessionInfo()
R version 3.3.1 (2016-06-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS release 6.8 (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] parallel  stats4    stats     graphics  grDevices utils     datasets 
[8] methods   base     

other attached packages:
 [1] ChIPQC_1.10.3              DiffBind_2.2.12           
 [3] SummarizedExperiment_1.4.0 Biobase_2.34.0            
 [5] GenomicRanges_1.26.4       GenomeInfoDb_1.10.3       
 [7] IRanges_2.8.2              S4Vectors_0.12.2          
 [9] BiocGenerics_0.20.0        ggplot2_2.2.1             

ADD REPLY
0
Entering edit mode

 

 

 

hi,

The original error looks to do with the way ChIPQC calculates cross coverage. The ChrM may cause some problems here because it is circular.

I believe there is a way to check for circular chromosomes so I will look to remove any circular chromosomes while providing a warning for the user.

For now the solution is to remove the mitochondrial chromosome.

I usually provide ChIPQC with a vector of the main chromosomes omitting ChrM and the contigs which haven't been fully assembled into main chromosomes.

Something like

> chromosomesToTest <- c(paste0("chr",c(1:22,"X","Y")))

> chromosomesToTest
 [1] "chr1"  "chr2"  "chr3"  "chr4"  "chr5"  "chr6"  "chr7"  "chr8"  "chr9"  "chr10" "chr11" "chr12"
[13] "chr13" "chr14" "chr15" "chr16" "chr17" "chr18" "chr19" "chr20" "chr21" "chr22" "chrX"  "chrY"

> my_experiment = ChIPQC(experiment="CellLine_Experiment_Sheet.csv", annotation="hg19",chromosomes=chromosomesToTest)

I will look to update ChIPQC to remove the ChrM internally in the next version.

best,

tom

ADD REPLY
0
Entering edit mode
@thomas-carroll-7019
Last seen 13 months ago
United States/New York/The Rockefeller …

hi,

Could you let me know what version of ChIPQC you are using? The result from a call to sessionInfo() should be good.

best,

tom

 

 

ADD COMMENT

Login before adding your answer.

Traffic: 551 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