I am performing an exploratory analysis of the ENCODE Txn Factor ChiPseq data. For some reason that I cannot understand, the covplot function only works for certain input files. It works perfectly for certain input files but when I change just the input file and maintain the script the same I get the following error. All the files I want to analyze come from the ENCODE databse and are .narrowpeaks or .broadpeaks format. Then I transform them to a GRanges file to analyze them.

> covplot <- covplot(my_granges,weightCol ="score") #coverage plot

Error in data.frame(chr = names(cov[i]), start = start(x), end = end(x),  : 
  arguments imply differing number of rows: 1, 118, 119

For example, when I use the ENCODE ChiPseq data for BCL11A in hESC.narrowpeaks or BCL11A in hESC Replicate 1.broadpeaks it works, but when I use BCL11A in hESC Replicate2.broadpeaks, it does not work.

please post the link to download this encode dataset, I will have a look at it.

Bests,

Guangchuang
 

It works right now if this: install_github("GuangchuangYu/ChIPseeker") Thanks man.

Dear Guangchuang,

Thanks for your reply!

 - This is replicate 1 of BCL11A ChIPseq in hESC from ENCODE (the one that works): 

https://www.genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=424063059_EqT7VUHDesA7tH02ocEyI12cm5kE&g=wgEncodeHaibTfbsH1hescBcl11aPcr1xPkRep1

- This is replicate 2 of BCL11A ChIPseq in hESC from ENCODE (the one that does not work):

https://www.genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=424063059_EqT7VUHDesA7tH02ocEyI12cm5kE&g=wgEncodeHaibTfbsH1hescBcl11aV0416102PkRep2

You have to click on "metadata", and then the link to download the file will appear.

Thank you very much for looking into this!

Best,

Sara

Thanks for your report. It has been fixed in version >=1.4.1

Please wait 2-3 days for Bioconductor to update tarball.
 

Thank you very much!

It works perfectly!