Entering edit mode

Ari Eszter
▴
40

@ari-eszter-5784
Last seen 7.7 years ago

Dear edgeR Users,
I have RNS-Seq count data by genes, 4 treatment groups (each with 2
biological replicates) and two factors (E: C or H, and A: 13 or 25).
And I
applied the GLM model of edgeR with these designs:
E, A, E+A and E*A
I would like to decide which model and which variant describes my data
better.
1. Am I right that the $coefficients of a "DGELRT" object are log
likelihood values?
2. Does anyone have a suggestion how can I apply AIC (Akaike
Information
Criterion) test or Likelihood ratio (LR) test on edgeR GLM results to
choose the best model?
3. Is there a possibility to plot the multiple GLM regression for a
gene?
(I would like plot something like this:
http://www.jerrydallal.com/LHSP/pix/regpix4.gif )
Bests,
Eszter
--
Eszter Ari
Institut für Populationsgenetik
Vetmeduni Vienna
Veterinärplatz 1
1210 Wien, Austria
[[alternative HTML version deleted]]