Comments (5)
Hi @jjyotikataria,
did you already have a look at the explanations of these metrics in our vignette and in our preprint on bioRxiv?
from fraser.
Yes, @c-mertes @ischeller . I went through the vignette and the paper. But I am not so much confident about what the splicing metrics sigh5, sigh3, theta3, theta5 are basically denoting. Also, can you give me a basic explanation regarding the p values and z scores.
from fraser.
@c-mertes @ischeller I went through almost all the papers and references. But I am not able to be understanding it completely. Can you please tell me what the p value, z score basically represents in context of the genes of the samples we are working upon. Also, what does psi of 0 to range 1 represents practically?
from fraser.
I still not really sure what part you did not understand regarding the p values and z scores, so I will give a brief basic explanation and hope that it helps:
Basically, psi5 quantifies acceptor site usage for a given donor site, and psi3 quantifies donor site usage for a given acceptor site. The interpretation of them is therefore on the intron-level, and not at the gene-level.
For a given intron with donor site D and acceptor site A, the metric psi5 indicates the ratio of the number of reads that map to this intron over the number of all reads that map to this intron or any other intron that also uses donor site D (see Figure 1 in our preprint for an illustration). For example, a psi5=0.8 for this intron means that 80% of split reads that contain the donor site D use the acceptor site A, while 20% of those split reads use another acceptor site. As the psi5 values are ratios, they can take values between 0 and 1 and the psi5 values of all introns that share a donor site sum up to 1.
Therefore, an aberrant psi5 value in one sample indicates aberrant acceptor site usage (and analogously aberrant donor site usage for aberrant psi3 values).
In FRASER, we model the expected value of psi3/psi5/theta for every intron, and then compute p values based on the beta-binomial distribution by comparing the observed read count of the intron with this expected psi3/psi5/theta value to assess the statistical significance of the observed count. We also provide delta psi3/psi5/theta values (i.e. the difference between observed and expected values). We don't recommend to use the z scores to find aberrant events, we mainly include them for comparison with other methods. Rather use the delta psi3/psi5/theta value (which are easy to interpret) together with the p value to find aberrant splicing events that are of interest to you.
For convenience, we also provide aggregated p values per gene. A significant gene-level p value indicates that this gene contained at least one intron with an aberrant psi3/psi5/theta values. But the psi3/psi5/theta values themselves can only be interpreted on the intron-level.
I hope this helps to answer your questions.
from fraser.
@ischeller Very elaborately explained. Thank you. I clearly understood now what psi5 and psi3 ratios indicate. I have another question dear. Like you said, "Therefore, an aberrant psi5 value in one sample indicates aberrant acceptor site usage (and analogously aberrant donor site usage for aberrant psi3 values." How do one determine and tag that particular psi value be it psi5 or psi3 as aberrant?
Also, observed or the experimental psi values are the values which we get from the real time data by giving the required bam files. What are exactly the expected psi values and what is the criteria for this expectation
It will make clear what delta psi is i.e. difference between actual observed and the expected values.
from fraser.
Related Issues (20)
- Error in creating dir HOT 3
- Question regarding the sample size of FRASER HOT 2
- Question about using the mergeCounts function
- Ppadjust is always 1 HOT 5
- Paired-end BAM files recognized as single end
- FRASER maxing out memory on sge cluster HOT 3
- Remove group and condition from the code/vignette as it is not used within FRASER
- filter out reads with MAPQ=0 & duplicated reads
- psiValue calculation
- annotateRangesWithTxDb function error HOT 1
- "Please provide hgnc symbols to compute gene p values!" error HOT 1
- Error in H5Fcreate(file) : HDF5. File accessibilty. Unable to open file.
- Error in working directory
- ERROR: Paired-end reads were detected in single-end read library
- A question about resultsByGenes
- FRASER Error in check_returned_array
- FRASER Resources and Run Time
- Offline use of countSplitReads HOT 2
- FRASER with non model organisms
- Error in predictYCpp(as.matrix(H), D, b) : Mat::init(): requested size is too large; suggest to enable ARMA_64BIT_WORD
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 fraser.