[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