Error in ChIPQC(samples) : Unable to process. Each bam file must be associated with at most one peakset.

Question

[bioinfo89]

Hi All,

I get this error when I run the ChIPQC package on my data (single sample). The commands are mentioned below:

## Load sample data
samples <- read.csv('meta/samplesheet.csv')
#View(samples)

## Create ChIPQC object
chipObj <- ChIPQC(samples, annotation="hg19")

The samplesheet.csv:

SampleID,Tissue,Factor,Replicate,bamReads,Peaks
SRR,PCa,AR,1,data/bams/SRR.sambamba.sorted.uniquemapped.bam,data/peakcalls/LNCaP_AR_macs2_peaks.narrowPeak

I am not sure what exactly am I doing wrong.

Any help will be much appreciated. Thanks!

Answer 1

[wing123]

Add another sample could solve this problem, which works for me.