Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Filter average 0 junction read on the plot #125

Open
Qingru44004 opened this issue Mar 26, 2024 · 1 comment
Open

Filter average 0 junction read on the plot #125

Qingru44004 opened this issue Mar 26, 2024 · 1 comment

Comments

@Qingru44004
Copy link

Hi,

As discussed in another issue, sometimes the junction read may round to 0 if its average is less than 1. Now, I want to add a filter of the events, like for every JC shown in the plot, at least XX% of the samples should have XX reads. How should I get those raw records to get rid of the 0 in the plot? Meanwhile, do you have a suggestion for the cutoff on my filteration?

@EricKutschera
Copy link
Contributor

You could try modifying the rmats2sashimiplot code. This post has links to the lines of code for calculating and ploting the junction counts: #112 (comment)

rMATS v4.3.0 includes --individual-counts which outputs counts for each junction. You could filter based on those columns

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants