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
ADD COMMENT
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)))

ADD REPLY
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!

ADD COMMENT

Login before adding your answer.

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