genas {limma}R Documentation

Genuine Association of Gene Expression Profiles

Description

Calculates biological correlation between two gene expression profiles.

Usage

genas(fit, coef=c(1,2))

Arguments

fit

an MArrayLM fitted model object produced by lmFit or contrasts.fit and followed by eBayes

coef

numeric vector of length 2 to indicate which contrasts/columns in the fit object are to be used

Details

The biological correlation between the true log fold changes of pairs of comparisons is computed. This method is to be applied when multiple groups (such as treatment groups, mutants or knock-outs) are compared back to the same control group.

This method is an extension of the empirical Bayes method of limma. It aims to separate the technical correlation, which comes from comparing multiple treatment/mutant/knock-out groups to the same control group, from biological correlation, which is the true correlation of the gene expression profiles between two treatment/mutant/knock-out groups.

The fit object should include only differentially expressed genes. One approach is to calculate the true proportion of differentially expressed genes using convest on the F p-value produced by lmFit. Any reasonable set of genes displaying some degree of differential expression should be adequate.

Value

genas produces a list with the following components.

technical.correlation

estimate of the technical correlation

covariance.matrix

estimate of the covariance matrix from which the biological correlation is obtained

biological.correlation

estimate of the biological correlation

deviance

the likelihood ratio test statistic used to test whether the biological correlation is equal to 0

p.value

the p.value associated with deviance

Author(s)

Belinda Phipson and Gordon Smyth

See Also

lmFit, eBayes, contrasts.fit

Examples

library(limma)
#  Simulate gene expression data,
#  6 microarrays with 100 genes on each array 
set.seed(2004)
y<-matrix(rnorm(600),ncol=6)

# two experimental groups and one control group with two replicates each
group<-factor(c("A","A","B","B","control","control"))
design<-model.matrix(~0+group)
colnames(design)<-c("A","B","control")

# fit a linear model
fit<-lmFit(y,design)
contrasts<-makeContrasts(A-control,B-control,levels=design)
fit2<-contrasts.fit(fit,contrasts)
fit2<-eBayes(fit2)

# calculate biological correlation between the gene expression profiles of (A vs control) and (B vs control)
genas(fit2)

[Package limma version 3.10.2 Index]