Hi,
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 )
Thanks,
Eszter
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 )
Thanks,
Eszter
Comment