Question: CAGEr: consensusCluster doesn't have dominant CTSS info
gravatar for Jingwen
10 months ago by
Jingwen10 wrote:


I am using the latest CAGEr (1.22.1) in R 3.5.0

When following the demo in CAGEr website, I found the consensusClusterGR function didn't return dominant CTSS information when specifying sample name:

> consensusClustersGR( ce, sample = "zf_unfertilized_egg", returnInterquantileWidth = TRUE, qLow = 0.1, qUp = 0.9)

The output is exactly the same as the output in the demo webpage, no dominant CTSS. Furthermore, when exporting consensusCluster to bed file:

> exportToBed(object = ce, what = "consensusClusters", qLow = 0.1, qUp = 0.9,oneFile=FALSE)

I got such error message: Error in consensusClustersQuantileLow(object) : Not supported for ‘CAGEexp’ objects. Use ‘consensusClustersQuantile()’ instead.

Then I tried with the original CAGEr format, CAGEset. There was still no dominant ctss info when calling consensusClusterGR function. Again error messge when exporting consensusCluster to bed files.

> exportToBed(object = ce, what = "consensusClusters", qLow = 0.1, qUp = 0.9, oneFile = FALSE)
Error in dominant_ctss_start[clusters.q$dominant_ctss == clusters.q$start] <- clusters.q$start[clusters.q$dominant_ctss ==  : 
  NAs are not allowed in subscripted assignments

Did I miss some critical step before exporting consensus clusters? Thanks for your help.



cager • 127 views
ADD COMMENTlink written 10 months ago by Jingwen10

Thanks Jingwen, I will correct this. By the way, I am also tracking this issue on GitHub: If you would like to submit a patch you are very welcome !

ADD REPLYlink written 10 months ago by Charles Plessy90
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 333 users visited in the last hour