Error when using ChIPQC package
1
1
Entering edit mode
2.6 years ago
owen ▴ 10

Hi, I am using ChIPQC package for assessing the quality of ChIP-Seq data. However, when I was trying to constructing a ChIPQCexperiment object, there was always an error:

chipObj = ChIPQC(samples, annotation = "hg19", blacklist = NULL, chromosome = "chr1")

FA_1 G3 FA endogenous basal 1 narrow
FB_1 G3 FB endogenous basal 1 narrow
FC_1 G29 FC endogenous basal 1 narrow
Compiling annotation...
Error in names(samplelist) <- unique(c(rownames(meta), names(controlist))) :
'names' attribute [6] must be the same length as the vector [5]


Did anyone know how to deal with it? Thank you.

ChIP-Seq qc R • 1.7k views
1
Entering edit mode

You seem to have either more samples than sample names or vice versa. What is the output of length(samples) and unique(c(rownames(meta), names(controlist)))

1
Entering edit mode
19 months ago

I also spent almost 2 whole days trying to figure out this problem. After remapping the bams, the problem went away! In my case, I got the bam files from one of my colleagues, which was made quite some time ago. I thought, it was fine! But, I learned my lesson! Never trust bam files coming from another person!