fit.limma: Wrapper script for lmFit, contrast.fit and ebayes

View source: R/limma.R

fit.limmaR Documentation

Wrapper script for lmFit, contrast.fit and ebayes

Description

fit.limma fits a limma model, computes contrasts and applies an emperical bayes procedure to shrink variance estimates

Usage

fit.limma(data, MM, CM)

Arguments

data

An expression set

MM

A model matrix

CM

A contrast matrix

Author(s)

Florian Klinglmueller float@lefant.net

Examples


data <- matrix(rnorm(600),ncol=6)
design <- data.frame('Treatment' = factor(rep(0:1,each=3)))
mm <- model.matrix(~0+Treatment,design)
colnames(mm) <- c('control','treatment')
cm <- makeContrasts('control-treatment',levels=mm)
out <- fit.limma(data,mm,cm)

topTable(out)


floatofmath/bt88.03.704 documentation built on March 25, 2022, 12:28 p.m.