Question: rowVars function, used in ballgown analysis not available for R 3.5.2, alternatives?
1
gravatar for RNAseqer
6 weeks ago by
RNAseqer 40
RNAseqer 40 wrote:

Hello,

I've hit a bit of a snag here: I had to update to the latest version of R to get ballgown working, but one of the functions I see included in filtering out low expression genes, rowVars, is not available for this version of R.

Using this tutorial: https://rpubs.com/kapeelc12/Ballgown

The rowVars function is employed in filtering:

bg_filt = subset(bg,"rowVars(texpr(bg)) >1",genomesubset=TRUE)

I have R studio, would it be possible to download an older version of R and switch to it temporarily to use rowVars, then switch back to continue the ballgown output analysis?

If so does anyone know what the last version of R rowVars worked on?

Or should I try to find an older version of ballgown, and downgrade my R version? I know there are many functions on bioconductor that do not necessarily work on the latest version of R, so is there a benifit to running older versions of R generally when working with these packages?

rowvars rna-seq ballgown R • 188 views
ADD COMMENTlink modified 6 weeks ago by kristoffer.vittingseerup1.7k • written 6 weeks ago by RNAseqer 40
5
gravatar for Bastien Hervé
6 weeks ago by
Bastien Hervé4.0k
Limoges, CBRS, France
Bastien Hervé4.0k wrote:

Is it the rowVars function from metaMA ?

install.packages('metaMA')
library(metaMA)
?rowVars
ADD COMMENTlink written 6 weeks ago by Bastien Hervé4.0k

Huh, must be! I tried to get it as part of matrixstats but that didnt work, this download did. Thanks a bunch! This was a real help.

ADD REPLYlink written 6 weeks ago by RNAseqer 40

This comment has been moved to an answer so it can get accepted and mark this question as solved.

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

ADD REPLYlink written 6 weeks ago by WouterDeCoster38k
4
gravatar for kristoffer.vittingseerup
6 weeks ago by
European Union
kristoffer.vittingseerup1.7k wrote:

I would not recommend using ballgown in the first place - it performs VERY poorly compared to all the other DE tools out there (and have not been updated in years). Take a look at this recent benchmark (which agrees with most other benchmarks but is one of the few that actually test ballgown). Instead go for DESeq2 or edgeR - if you have StringTie data you can directly get them into R via tximport and use use DESeq2 as described here.

ADD COMMENTlink modified 25 days ago • written 6 weeks ago by kristoffer.vittingseerup1.7k
1

This is going to save me some time and heartache. Thanks for the heads-up!

ADD REPLYlink written 4 weeks ago by RNAseqer 40
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 657 users visited in the last hour