[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