Entering edit mode
12 months ago
mavy
▴
10
Hello I want to check the quality of peaks using ChipQC , but when I am trying to create the object its throwing an error
code I am using
chipObj <- ChIPQC(samples, annotation="hg38")
and the error is CR h3k27me3 epc1 NA bed
2 CS h3k27me3 ep NA bed
3 AR h3k27me3 EP_Aza NA bed
Error in ChIPQC(samples, annotation = "hg38") :
Unable to process. Each bam file must be associated with at most one peakset.
Whereas I have given one bed file for each sample.
my sample file holds the following information. Sample_ID Tissue Factor Condition bam reads ContriID bamControl Peaks Peakcaller
Any help will be appreciated.