Hello,
My name is Javier Simón Seanchez and I'm a post-doc at the VUmc in
Amsterdam.
The reason of this e-mail is that I'm trying to run edgeR in an expression
dataset and getting the following error when calculating the GLM common
dispersion:

*Error in beta[k, ] <- betaj[decr, ] :
  NAs are not allowed in subscripted assignments
*

Im running cases versus controls and I want to modulate for different
tissues. How can I overcome this error?

My design matrix looks like this:
    (Intercept)    Tissue1    Tissue2    Tissue3    Tissue4    case_control
Sample1    1    0    0    0    0    1
Sample2    1    1    0    0    0    1
Sample3    1    0    0    1    0    1
Sample4    1    0    0    0    1    1
Sample5    1    0    0    0    0    1
Sample6    1    1    0    0    0    1
Sample7    1    0    1    0    0    1
Sample8    1    0    0    1    0    1
Sample9    1    0    0    0    1    1
Sample10    1    0    0    0    0    1
Sample11    1    1    0    0    0    1
Sample12    1    0    1    0    0    1
Sample13    1    0    0    1    0    1
Sample14    1    0    0    0    1    1
Sample15    1    0    0    0    0    1
Sample16    1    1    0    0    0    1
Sample17    1    0    1    0    0    1
Sample18    1    0    0    1    0    1
Sample19    1    0    0    0    1    1
Sample20    1    0    0    0    0    1
Sample21    1    1    0    0    0    1
Sample22    1    0    1    0    0    1
Sample23    1    0    0    1    0    1
Sample24    1    0    0    0    1    1
Sample25    1    0    0    0    0    1
Sample26    1    1    0    0    0    1
Sample27    1    0    1    0    0    1
Sample28    1    0    0    1    0    1
Sample29    1    0    0    0    1    1
Sample30    1    0    0    0    0    1
Sample31    1    1    0    0    0    1
Sample32    1    0    1    0    0    1
Sample33    1    0    0    1    0    1
Sample34    1    0    0    0    1    1
Sample35    1    0    0    0    0    1
Sample36    1    1    0    0    0    1
Sample37    1    0    0    1    0    1
Sample38    1    0    0    0    1    1
Sample39    1    0    0    0    0    1
Sample40    1    1    0    0    0    1
Sample41    1    0    0    1    0    1
Sample42    1    0    0    0    1    1
Sample43    1    0    1    0    0    1
Sample44    1    0    0    0    0    1
Sample45    1    1    0    0    0    1
Sample46    1    0    1    0    0    1
Sample47    1    0    0    1    0    1
Sample48    1    0    0    0    1    1
Sample49    1    0    0    0    0    1
Sample50    1    1    0    0    0    1
Sample51    1    0    1    0    0    1
Sample52    1    0    0    1    0    1
Sample53    1    0    0    0    1    1
Sample54    1    0    0    0    0    1
Sample55    1    1    0    0    0    1
Sample56    1    0    1    0    0    1
Sample57    1    0    0    1    0    1
Sample58    1    0    0    0    1    1
Sample59    1    1    0    0    0    1
Sample60    1    0    0    0    1    1
Sample61    1    0    0    0    0    1
Sample62    1    1    0    0    0    1
Sample63    1    0    1    0    0    1
Sample64    1    0    0    1    0    1
Sample65    1    0    0    0    1    1
Sample66    1    1    0    0    0    1
Sample67    1    0    1    0    0    1
Sample68    1    0    0    0    0    1
Sample69    1    0    0    1    0    1
Sample70    1    0    0    0    1    1
Sample71    1    0    0    0    0    0
Sample72    1    1    0    0    0    0
Sample73    1    0    1    0    0    0
Sample74    1    0    0    1    0    0
Sample75    1    0    0    0    1    0
Sample76    1    0    0    0    0    0
Sample77    1    1    0    0    0    0
Sample78    1    0    1    0    0    0
Sample79    1    0    0    1    0    0
Sample80    1    0    0    0    1    0
Sample81    1    0    0    0    0    0
Sample82    1    1    0    0    0    0
Sample83    1    0    1    0    0    0
Sample84    1    0    0    1    0    0
Sample85    1    0    0    0    1    0
Sample86    1    0    0    0    0    0
Sample87    1    1    0    0    0    0
Sample88    1    0    1    0    0    0
Sample89    1    0    0    1    0    0
Sample90    1    0    0    0    1    0
Sample91    1    0    0    0    0    0
Sample92    1    1    0    0    0    0
Sample93    1    0    1    0    0    0
Sample94    1    0    0    1    0    0
Sample95    1    0    0    0    1    0

Thanks a lot for your help.

	[[alternative HTML version deleted]]

