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

Feature request: subset does not handle exon and intron expression #147

Open
ghost opened this issue Sep 11, 2018 · 2 comments
Open

Feature request: subset does not handle exon and intron expression #147

ghost opened this issue Sep 11, 2018 · 2 comments

Comments

@ghost
Copy link

ghost commented Sep 11, 2018

I have been getting an error when performing stattest:
Error in solve.default(t(mod) %*% mod) :
system is computationally singular: reciprocal condition number =
This happens in particular when performing a stattest for differential exon expression:

bg_data_filt_e <- ballgown::subset(bg_data, "rowVars(eexpr(bg_data, meas='mcov')) > 1", genomesubset=TRUE)
results_exons <-  stattest(bg_data_filt_e, feature='exon', covariate='Condition' , `meas='mcov')

I have looked at results_exons and it appears as if my filter using eexpr is not actually working in the context of the condition argument. Similarly it does not appear as if you can filter intron expression either (using rowVars(iexpr... ).

If I want to filter exon expression measurements to remove low variance/low expression exons is there another way to do it?
I am trying to export the exon expression data (using eexpr) and filter on low variance (not in a ballgown object). I was thinking I could then perform stattest using the gowntable argument. Since I am very new to R I am finding this difficult, not least because during filtering I am losing rownames which contain the e_id (exon ids).
I've also posted on bioconductor forums: https://support.bioconductor.org/p/112824/#112917

@alyssafrazee alyssafrazee changed the title Filtering exon and intron expression Feature request: subset does not handle exon and intron expression Sep 11, 2018
@alyssafrazee
Copy link
Owner

Hi there, thanks for the note! I responded on Bioconductor with a workaround that should unblock you.

Currently ballgown is maintenance mode only (bugfixes only, no new features) but if developers do eventually pick things back up, I've marked this as a feature request so someone can take it on (and for anyone reading, pull requests are welcome!).

@ghost
Copy link
Author

ghost commented Sep 13, 2018

Hi Alyssa. Thanks for the comment. I have added a reply on the bioconductor forums as I am still experiencing difficulties with my stattest. If you have time to take a look that would be much appreciated.

@ghost ghost closed this as completed Sep 25, 2018
@ghost ghost reopened this Sep 25, 2018
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

1 participant