[BioC] limma lmFit error
J.delasHeras at ed.ac.uk
J.delasHeras at ed.ac.uk
Wed Feb 14 15:20:50 CET 2007
I think it means you stated that "D2R" was your reference, as in a
common reference design, where every slide uses that same reference,
but your design is not like that: your D4 and D7 samples are
completely unlinked to D2R, therefore limma can only estimate D2.
You may want to look at Chapter 9 in the limma users guide (single
channel analysis of two-color data), p50.
Jose
Quoting Ben Nota <ben.nota at ecology.falw.vu.nl>:
> Dear BioC,
>
> I'm getting an error using the lmFit function within limma. I'm using R
> 2.4.0 with limma version 2.9.1.
>
> A little background about the experiment:
>
> I use two color agilent arrays containing duplicate spots. I have metal
> exposure compared to non exposed (= reference) at 3 different time
> points (2, 4 and 7 days). Four bio-replicates per time point.
>
> Explanation of RNA sample names: Day 2 reference (D2R) vs. Day 2 Metal
> (D2M), Day 4 reference (D4R) vs. Day 4 Metal (D4M), etc.
>
> I normalized my data using loess
>
> My targets look like this:
>
> > targets
> Array FileName Cy3 Cy5 Date
> 1 Array15A Array15A.txt D4R D4M 31/01/07
> 2 Array15B Array15B.txt D4M D4R 31/01/07
> 3 Array17A Array17A.txt D4R D4M 31/01/07
> 4 Array17B Array17B.txt D4M D4R 31/01/07
> 5 Array46A Array46A.txt D7R D7M 31/01/07
> 6 Array46B Array46B.txt D7M D7R 31/01/07
> 7 Array47A Array47A.txt D7R D7M 31/01/07
> 8 Array47B Array47B.txt D7M D7R 31/01/07
> 9 Array10A Array10A.txt D2M D2R 17/10/06
> 10 Array11A Array11A.txt D2R D2M 17/10/06
> 11 Array12A Array12A.txt D2M D2R 16/11/06
> 12 Array12B Array12B.txt D2R D2M 16/11/06
>
> My design:
>
> > design<-modelMatrix(targets,ref="D2R")
> Found unique target names:
> D2M D2R D4M D4R D7M D7R
> > design
> D2M D4M D4R D7M D7R
> 1 0 1 -1 0 0
> 2 0 -1 1 0 0
> 3 0 1 -1 0 0
> 4 0 -1 1 0 0
> 5 0 0 0 1 -1
> 6 0 0 0 -1 1
> 7 0 0 0 1 -1
> 8 0 0 0 -1 1
> 9 -1 0 0 0 0
> 10 1 0 0 0 0
> 11 -1 0 0 0 0
> 12 1 0 0 0 0
>
> And my correlation is calculated:
>
> > cor<-duplicateCorrelation(MA,design,ndups=2)
> > cor$consensus.correlation
> [1] 0.8890682
>
> But if I want to fit my linear model I'll get the following error message:
>
> > fit<-lmFit(MA, design, ndups=2,cor=cor$consensus.correlation)
> Coefficients not estimable: D4R D7R
> Error in dimnames(cov.coef) <- list(coef.names, coef.names) :
> length of 'dimnames' [1] not equal to array extent
>
> Does anybody knows what this error means and what I'm doing wrong? I
> hope so!
>
> Thanks in advance!
>
> Best regards, Ben.
>
> --
> Benjamin Nota
> Vrije Universiteit
> Department of Animal Ecology
> De Boelelaan 1085
> 1081 HV AMSTERDAM, The Netherlands
> Tel: +31 (0)20-5987217
> Fax: +31 (0)20-5987123
> ben.nota at falw.vu.nl
>
> http://www.bio.vu.nl/do/
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
--
Dr. Jose I. de las Heras Email: J.delasHeras at ed.ac.uk
The Wellcome Trust Centre for Cell Biology Phone: +44 (0)131 6513374
Institute for Cell & Molecular Biology Fax: +44 (0)131 6507360
Swann Building, Mayfield Road
University of Edinburgh
Edinburgh EH9 3JR
UK
More information about the Bioconductor
mailing list