[BioC] Limma error

Amin Momin gtg269c at mail.gatech.edu
Fri Nov 21 16:00:29 MET 2003


Hi ,

    I have been trying to use limma to construct linear model of my
normalized affymetrix data, but the program gives me an error at the
last step of the calculation by empirical bayes. Can someone using the
package give a comment on my error.


> design <- model.matrix(~ -1+factor(c(1,2,3,4,5,1,2,3,4,5)))
> colnames(design) <- c("group1", "group2", "group3", "group4",
"group5")
> fit <- lmFit(esetsarco, design)
> contrast.matrix <- makeContrasts(group1-group2, group1-group3,
group1-group4, group1-group5, group2-group3, group2-group4,
group2-group5, group3-group4, group3-group5, levels=design)
> fit2 <- contrasts.fit(fit, contrast.matrix)
> fit2 <- eBayes(fit2)
Error in ebayes(fit = fit, proportion = proportion, std.coef = std.coef)
: 
        No residual degrees of freedom in linear model fits



   Thanks 

   Amin Momin



More information about the Bioconductor mailing list