Question: contrasts in limma and missing values
gravatar for Gordon Smyth
11.9 years ago by
Gordon Smyth38k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth38k wrote:
Dear Hans-Ulrich, This behaviour follows directly from arithmetic in R, for example c(1,0) %*% c(1,NA) produces NA rather than 0. The reason for this is that the unknown value corresponding to the NA could in principle have been Inf, meaning that 0*NA can't be inferred to be 0. You could argue that limma should assume that infinite expression values are impossible, therefore it should treat 0*NA as 0, but I haven't so far been convinced that this finesse is worth the effort. In your example, you could get the contrast you want without NAs simply by telling limma explicitly that the 3rd coefficient doesn't appear in your contrast, which you can do by fit2 <-[,-3],c(1,-1)) Best wishes Gordon >Date: Mon, 29 Oct 2007 19:42:27 +0100 >From: Hans-Ulrich Klein <h.klein at=""""> >Subject: [BioC] contrasts in limma and missing values >To: BioC Mailing List <bioconductor at=""""> >Message-ID: <47262993.8030201 at> >Content-Type: text/plain; charset=ISO-8859-1; format=flowed > >Dear List, > >I have a question concerning missing values in limma. Consider the >following short example: > >M = rbind(c(10,9,5,4,5,12,13,7,8),c(10,9,5,4,5,12,13,7,8)) >w = rbind(c(1,1,1,1,1,1,1,1,1),c(1,1,1,1,1,0,0,0,0)) >class = factor(c(1,1,2,2,2,1,1,2,2)) >lab = factor(c(1,1,1,1,1,2,2,2,2)) > >design = model.matrix(~ 0 + class + lab) >contrasts = makeContrasts(contrasts="class1-class2", levels=design) > >fit = lmFit(M, design=design, weights=w) >fit2 =, contrasts) > >Obviously, the effect of "lab" can not be estimated for gene 2. All >measurements of gene 2 in lab2 have weights = 0. However, the contrast >of interest (here class1 - class2) can be estimated by using the >observations from lab1, but in the limma-object (fit2 in the example >above) the contrast is NA. Can someone explain me whether/why this is >sensible? > >Thanks for your help, >Hans-Ulrich
limma • 562 views
ADD COMMENTlink written 11.9 years ago by Gordon Smyth38k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 135 users visited in the last hour