predFCm: Predictive log fold change for microarrays
In limma: Linear Models for Microarray Data

Description Usage Arguments Details Value Author(s) References See Also Examples

Calculate the predictive log fold change for a particular coefficient from a fit object.

1	predFCm(fit, coef=2, var.indep.of.fc=TRUE, all.de=TRUE, prop.true.null.method="lfdr")

`fit`	an `MArrayLM` fitted model object produced by `lmFit` and `eBayes`
`coef`	integer vector indicating which columns in the fit object are to be shrunk
`var.indep.of.fc`	assume the genewise variances are independent of genewise fold changes?
`all.de`	assume all genes are have a non-zero true fold change (`TRUE`)? If `FALSE`, then the proportion of truly non-differentially (non-DE) genes expressed will be estimated.
`prop.true.null.method`	method used to estimate proportion of truly non-DE genes. See `propTrueNull` for possible values.

The predictive log fold changes are calculated as the posterior mean log fold changes in the empirical Bayes hierarchical model. We call them predictive log fold changes because they are the best prediction of what the log fold change will be for each gene in a comparable future experiment.

The log fold changes are shrunk towards zero depending on how variable they are. The var.indep.of.fc argument specifies whether the prior belief is that the log fold changes are independent of the variability of the genes or whether the log fold changes increase with increasing variability of the genes.

If all.de=TRUE, then all genes are assumed to have a non-zero log fold change, even if quite small. If all.de=FALSE, then some genes are assumed to have log fold changes exactly zero. The proportion of non-DE genes is estimated and taken into account in the calculation.

numeric vector of predictive (shrunk) log fold changes

Belinda Phipson and Gordon Smyth

Phipson, B. (2013). Empirical Bayes modelling of expression profiles and their associations. PhD Thesis. University of Melbourne, Australia. http://repository.unimelb.edu.au/10187/17614

lmFit, eBayes, contrasts.fit

#  Simulate gene expression data,
#  6 microarrays with 1000 genes on each array 
set.seed(2004)
y <- matrix(rnorm(6000),ncol=4)

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

# fit a linear model
fit <- lmFit(y,design)
fit <- eBayes(fit)

# output predictive log fold changes for first 5 genes
pfc <- predFCm(fit,coef=2)