[R] (no subject)
Noha Osman
nmo_138 at usc.edu
Sat Jan 3 22:25:34 CET 2015
Hello everyone,
I try to use edge R package to analysis my data as following
This is my data frame called subdata
A B C D E
1 13707 13866 12193 12671 10178
2 0 0 0 0 1
3 7165 5002 1256 1341 2087
6 8537 16679 9042 9620 19168
10 19438 25234 15563 16419 16582
16 3 3 11 3 5
genotype=factor(c("LMS1","LRS1","MS3","MS4","RS5"))
y=DGEList(counts=data.matrix(subdata),group=genotype)
design=model.matrix(~0+group,data=subdata)
y=estimateGLMCommonDisp(y,design, verbose=TRUE)
I try to calculate the common dispersion to estimate up and down regulation genes but I get this error message as following
Warning message:
In estimateGLMCommonDisp.default(y = y$counts, design = design, :
No residual df: setting dispersion to NA
please could any one help me to dissolve this problem , I really appreciate for that
[[alternative HTML version deleted]]
More information about the R-help
mailing list