CEDA (version 1.1.1)

runLimma: Modeling CRISPR screen data by R package limma

Description

The lmFit function in R package limma is employed for group comparisons.

Usage

runLimma(data, design, contrast.matrix)

Value

A data frame with rows corresponding to sgRNAs and columns corresponding to limma results

Arguments

data

A numeric matrix containing log2 expression levels of sgRNAs with rows corresponding to sgRNAs and columns corresponding to samples.

design

A design matrix with rows corresponding to samples and columns corresponding to coefficients to be estimated.

contrast.matrix

A matrix with columns corresponding to contrasts.

Examples

y <- matrix(rnorm(1000*6),1000,6)
condition <- gl(2,3,labels=c("Treatment","Baseline"))
design <- model.matrix(~ 0 + condition)
contrast.matrix <- makeContrasts("conditionTreatment-conditionBaseline",levels=design)
limma.fit <- runLimma(y,design,contrast.matrix)