[BioC] Limma : Problem with lmFit() and block

Laetitia Marisa Laetitia.Marisa at cgm.cnrs-gif.fr
Tue Apr 19 10:46:51 CEST 2005


Hello,

When I use the block variable in the lmFit, data in the MArrayLM have 
different dimension :

 > dim(analysis$fit)
[1] 6936    3

 > dim(analysis$fit$coef)
[1] 6936    3

 > dim(analysis$fit$genes)
[1] 2312    6


Here is the command line I've  passed and a description of the data I've 
used :

analysis$fit <- lmFit( msData, design=analysis$effects$Designmatrix , 
block = msData$block , ndups=msData$printer$ndups, 
spacing=msData$printer$spacing,  weights=NULL )

 > dim(msData)
[1] 6936   15
 > msData$block
 [1] 1 1 2 3 3 4 4 5 5 6 6 7 7 8 8
 > msData$printer$ndups
[1] 3
 > msData$printer$spacing
[1] 1

 > analysis$effects$Designmatrix
         A  B  C
Array1   1  0  0
Array2  -1  0  0
Array3   1  0  0
Array4   1  0  0
Array5  -1  0  0
Array6   0  1  0
Array7   0 -1  0
Array8   0  1  0
Array9   0 -1  0
Array10  0  0  1
Array11  0  0 -1
Array12  0  0  1
Array13  0  0 -1
Array14  0  0  1
Array15  0  0 -1

Am I doing something wrong or it is a problem in the lmFit function?

Thanks a lot.

Laetitia Marisa.



More information about the Bioconductor mailing list