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]]

