CAGEr: Quantile Values (q_0.1 and q_0.9) in CAGEset / CAGEexp objects
1
1
Entering edit mode
mirkocelii ▴ 40
@mirkocelii-23498
Last seen 4.4 years ago

Hello, Are q0.1 and q0.9 values the same in CAGEexp and CAGEset objects? To cut a long story short, I'm trying to run a code in which they used CAGEset objects, but taking advantage of CAGEr vignette I created a CAGEexp object. the CAGEseq object value os a position of $q0.1 and $q0.9 from base n=1 of the chromosome (so big numbers close to those in $start $end columns) the CAGEexp object I created has very small values, like 1, 10, 25, they look like the distance from the base n=1 of the interval defined by start and end columns. is it like that? to get the InterquantileWidth tagClusters() should I just edit those columns by doing

 x$"q_0.1" = x$"start" + x$"q_0.1" - 1
 x$"q_0.9" = x$"start" + x$"q_0.9" - 1

many thanks !

CAGEr CAGEset quantile CAGEexp • 1.1k views
ADD COMMENT
0
Entering edit mode
@charles-plessy-7857
Last seen 15 months ago
Japan

Hello again,

CAGEset and CAGEexp objects differ in the way they internally store quantile information, by both of them should return the same values with commands such as CAGEr:::tagClustersQuantileLow(exampleCAGEexp, 1, "0.1"). (Please note that exampleCAGEexp and exampleCAGEset are not the same data).

CAGEexp objects store quantile positions relatively to the start of the cluster. This allows for efficient Rle compression, especially for small quantiles or with unfiltered clusters where many of them have a size of 1.

ADD COMMENT
0
Entering edit mode

Dear Charles, thanks again for your support!

"CAGEexp objects store quantile positions relatively to the start of the cluster."

Cool this is what i have.

"CAGEset and CAGEexp objects differ in the way they internally store quantile information"

So, you confirm that CAGEsetq0.1 = CAGEexpstart + CAGEexpq_0.1 -1, right?

ADD REPLY
0
Entering edit mode

I think so. By the way, in that development round (Bioc 3.12), I definitely want to sanity check the quantile information code of CAGEexp objeccts, and then remove the CAGEset code. While doing so, we might change how quantile information is stored internally. If you have interest to help, please join us!

ADD REPLY

Login before adding your answer.

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