[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