[BioC] lmfit "No residual degrees of freedom"
David martin
vilanew at gmail.com
Mon Jan 17 14:09:36 CET 2011
Hello,
I'm having problems computing comparisons between different groups (this
is not array data and i'm not using limma). Any on what the error must be ?
> mydata.eset
ExpressionSet (storageMode: lockedEnvironment)
assayData: 2754598 features, 3 samples
element names: exprs
protocolData: none
phenoData
sampleNames: sampleA sampleB sampleC
varLabels: Group Sample
varMetadata: labelDescription
featureData: none
experimentData: use 'experimentData(object)'
Annotation:
fit <- lmFit(mydata.eset, design)
....
contrast.matrix
> contrast.matrix
group1-group2
group1 1
group2 -1
group3 0
group1-group3
group1 1
group2 0
group3 -1
group2-group3
group1 0
group2 1
group3 -1
fit2 <- contrasts.fit(fit, contrast.matrix)
fit2 <- eBayes(fit2)
Error in ebayes(fit = fit, proportion = proportion, stdev.coef.lim =
stdev.coef.lim) :
No residual degrees of freedom in linear model fits
thanks,
david
More information about the Bioconductor
mailing list