RNA-seq time series data analysis
1
0
Entering edit mode
@nataliafghf-11494
Last seen 8.2 years ago

Hi everybody,

I'm going to apply multiDE, an R package for detection of DEG in multiple treatment conditions, on some RNA-seq time series data... (I wanna assume each time point as a treatment condition)

Let's assume Yidg denotes the normalized read counts for sample "i", in condition "d" for gene "g".

We also assume that Yidg marginally follows the negative binomial distribution with expectation "μdg" and dispersion parameter "ϕg" (i.e., the variance of Yidg is μdg+ϕgμ2dg).

The statistical methodology behind this package is a two factorial log linear model : logμdg = μ + αd + βg + γdg = μ + αd + βg + UdVg,

where μ is the grand mean, αd is the main effect for condition d, βg is the main effect for gene g, and γdg:=UdVg is the interaction effect between gene g and condition d.

My professor has asked me to estimate the main effect for condition (α), the main effect for gene (β) and the effect of interaction between gene and condition (γ). While the package can only show "Ud" in its output...

I'm in grave need of someone to help me please find out how I can estimate those effects...

My main problem is I don't know how I can calculate μdg. Maybe if I can calculate it, then applying a regression strategy would be helpful to estimate those effects...

here it is the link to the full paper: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4917940/

Thanks in advance

RNA-seq Log-linear Main effects • 1.6k views
ADD COMMENT
1
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 16 hours ago
The city by the bay

multiDE does not seem to be a Bioconductor package. If you're having problems with using it, then I would suggest that you contact the package authors directly, as this is not the appropriate forum.

More generally, several Bioconductor packages can support identification of DE genes from RNA-seq time course data. Section 9.6 of the limma user's guide describes how to parametrize a design matrix for time course samples, and that knowledge can be applied to many well-tested pipelines for RNA-seq data analysis such as voom or edgeR. I would suggest, at the very least, trying out some of these methods.

ADD COMMENT
0
Entering edit mode

Thank you so much for your useful info... well, in fact, I have to run multiDE because my thesis has been defined on it... But I will try those mentioned packages too... Thanks

ADD REPLY

Login before adding your answer.

Traffic: 451 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6