Comments (8)
Hi Ruggero,
I assume you have found out about the RED/RLD scores? We use condition 1 as the reference. You can also see the "APA regulation scores" section in our paper.
Best,
Vivian
from maaper.
Thank you for letting me know! When I clicked on the Github link in the email notification, it somehow took me to this page.
from maaper.
Hi Ruggero,
Thanks for posting this question on Github.
The RLDi score compares expression of all intronic/internal PASs (weighted by length) vs. all PASs on the 3'-most exon, while the REDi score compares the intronic/internal PAS with the largest proportion change vs. the PAS with the largest proportion change on the 3' most exon. Therefore, it is possible that the two scores have opposite signs, especially when there are >2 PASs on the same gene, as you have observed.
For example, consider a gene with 3 intronic PASs and 1 PAS on the 3' most exon. Their proportions in cell type A are respectively:
40%, 30%, 20%, 10% (3' most PAS)
and their proportions in cell type B are:
34%, 24%, 28%, 14%
Then, with some constraints on the PAS positions, the RLDi and REDi of this gene can have opposite signs.
I hope this answers your questions.
Best,
Vivian
from maaper.
Hi Vivian,
Yes it does, thank you for your reply.
As a follow up question I would ask if there are any plans to report the number of reads per PAS in the PASs table in addition to the relative amount as percentage.
The reason I ask this is to individuate sites that have low coverage, as that may also contribute to ambiguous situations and hard to interpret results.
Thank you in advance :)
Ruggero
from maaper.
Hi Ruggero,
Given the current outputs, you can use the number of reads per gene and the PAS percentages to calculate the number of reads per PAS.
Best,
Vivian
from maaper.
Hi Vivian,
I see, so the number of reads corresponds to the sum of all the PASs represented, got it.
Thank you for the info :)
Best,
Ruggero
from maaper.
Hi Vivian,
Last question was from an user called AidenSb, not me. I did find the info about the RED/RLD scores and thank you for the info on the reference condition, though it wasn't me asking :)
Best,
Ruggero
from maaper.
I suspect he might have just deleted his comment, no harm done ^^
from maaper.
Related Issues (9)
- How do you determine fragment length HOT 4
- Individual Condition Tracks HOT 1
- Comparing celltypes from multiple samples. HOT 2
- running maaper for more than one sample HOT 9
- 1 fragments longer than 600 ... Error in density.default(dist_all, from = 0, to = dist_thre, n = dist_thre + : argument 'x' must be numeric HOT 6
- Error in density.default(dist_all, from = 0, to = dist_thre, n = dist_thre + : argument 'x' must be numeric
- How to study other species? HOT 1
- Effects of parameters HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from maaper.