Hello,

I have vehemently searched the Bioconductor support site for answers with no success, and hoping someone might be able to help me. I am running ChIPQC on fission yeast samples and while it gets through most of the samples when generating the object, I can't get past a certain point.

The error is listed below (out of bound indices and out of bound ranges), and I have no idea how to troubleshoot but I have some thoughts on why this might be happening.

1. I have tried running it with consensus = T and consensus =F, but still get the same error.

2. I do not have an annotation file and am wondering if this is the root of the problem? 

3. I have tried with and without columns for the bamControl and ControlID in my samples.csv file (as shown in the help ?ChIPQC), although the tutorial example does not have these columns. Is it recommended to have input file information?

3. When I load these same samples into DifffBind I am able to conduct the analyses just fine. Using DiffBind I see that the samples that are throwing an error, are those that have low FrIP (5-10%). Perhaps this is the problem? 

Any help with this would be great and much appreciated!

Many thanks,

Meeta

chipQC <- ChIPQC(experiment=samples, consensus = TRUE,
                 chromosomes=c('AB325691ASM294v2', 'IASM294v2','IIASM294v2',               'IIIASM294v2','MTASM294v2', 'MTRASM294v2'))

Error in sample WT-0h-replicate2: subscript contains NAs or out-of-bounds indices
Error in sample WT-0h-replicate3: subscript contains NAs or out-of-bounds indices
Error in sample WT-0h-replicate4: subscript contains NAs or out-of-bounds indices
Error in sample WT-24h-replicate4: subscript contains NAs or out-of-bounds indices
Error in sample sample43_S43_combined: subscript contains out-of-bounds ranges
Error in sample chipseq41_S41_combined: subscript contains out-of-bounds ranges
Error in ChIPQC(experiment = samples, consensus = TRUE, chromosomes = c("AB325691ASM294v2",  : 

Errors in generating sample data.

In addition: Warning messages:
1: In if (class(sample) != "ChIPQCsample") { :
  the condition has length > 1 and only the first element will be used
2: In if (class(sample) != "ChIPQCsample") { :
  the condition has length > 1 and only the first element will be used
3: In if (class(sample) != "ChIPQCsample") { :
  the condition has length > 1 and only the first element will be used
4: In if (class(sample) != "ChIPQCsample") { :
  the condition has length > 1 and only the first element will be used

> sessionInfo()
R version 3.2.5 (2016-04-14)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu precise (12.04.5 LTS)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               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    LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats4    parallel  stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] ChIPQC_1.6.1               DiffBind_1.16.3            locfit_1.5-9.1             GenomicAlignments_1.6.3   
 [5] Rsamtools_1.22.0           Biostrings_2.38.4          XVector_0.10.0             limma_3.26.8              
 [9] SummarizedExperiment_1.0.2 Biobase_2.30.0             GenomicRanges_1.22.4       GenomeInfoDb_1.6.3        
[13] IRanges_2.4.8              S4Vectors_0.8.11           BiocGenerics_0.16.1        ggplot2_2.1.0             
[17] RSQLite_1.0.0              DBI_0.3.1 

Hi Meeta, 

We need some more information to figure out what is happening. Is "samples" a samplesheet, a DiffBind (DBA) object, or a SampleList? Could you share more of your script and/or data (peak files etc.), perhaps send me a link at my email:

rory.stark@cruk.cam.ac.uk

Regards-

Rory