I'm wondering if anyone has any experience using continuous covariates in a JunctionSeq model for testing differential exon usage (DEU)? From the readJunctionSeqCounts function description it says: "design A data frame containing the condition variable and all desired covariates. All variables should be factors". Also, in the vignette all of the examples only involve factors.
My questions then are:
- Is there some way to coerce the program to include numeric factors for things like age or components outputted from SVA?
- If not, are there other programs that test for DEU that might allow for continuous covariates in the model?
•
written
by
dhibar •
1 answer
Unless you want to start altering the JunctionSeq code, I don't believe they've added the ability analyze continuous covariates.
For exon coverage, you could use DEXSeq. Or, if you use the QoRTs counts for JunctionSeq, you could analyze the exon and junction counts separately (so, junctions get a different dispersion value) with your counts-based differential expression algorithm of your choice, but you wouldn't have the visualizations provided by JunctionSeq.
Log in to answer this question.
-
Paired analysis with sex and age as covariates for SVA
written by Curious_Gene •
-
Function crr from cmprsk package with categorical variable
written by Nona •
-
What exactly should and shouldn't go in Combat's model matrix?
written by n,n
-
Tools for screening the influence/importance of covariates in multidimensional data
written by Papyrus
-
SVA variable of interest in eQTL analysis
written by Colari19
-
Surrogate variables in eQTL studies
written by rodd
-
Covariates for eQTL analysis
written by Colari19
-
How to in R: Include principal components as covariates in analysis
written by RNAseqer
-
SVA analysis output R
written by robertt_bw •
-
How to adjust data for "continuous" covariates ?
written by housebie •
