get_elbow_limma: Calculates fold values from an MArrayLM object.

Description Usage Arguments Details Value Examples

View source: R/elbowlib.R

Description

Performs the ELBOW fold change test on an MArrayLM object. This is a wrapper class to help integrate the ELBOW method into Bioconductor. followed tutorial from: http://bioinformatics.knowledgeblog.org/2011/06/20/analysing-microarray-data-in-bioconductor/

Usage

1
  get_elbow_limma(marraylm, columns = NULL)

Arguments

marraylm

is the MArrayLM object to analyze.

columns

is the list of sample columns to obtain the elbow fold cut-off values for. This can be specified as a vector or a single value.

Details

Some of the code in this method is based on how toptable accesses the MArrayLM object to read fold values. Therefore, any MArrayLM object which works with the toptable method should also work with this method.

Value

a matrix specified as follows

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
########
		# LOAD DATA INTO LIMMA
		########
		library("limma")

		# load a filtered expression set into R
		# NOTE: see the vignette for instructions on preparing
		#       a filtered dataset with your own data.
		data(GSE20986_eset_exprs, package="ELBOW")
		data(GSE20986_design, package="ELBOW")

		# fit the linear model to the filtered expression set
		fit <- lmFit(GSE20986_eset_exprs, GSE20986_design)

		# set up a contrast matrix to compare tissues v cell line
		contrast.matrix <- makeContrasts(huvec_choroid = huvec - choroid, huvec_retina = huvec - retina, huvec_iris <- huvec - iris, levels=GSE20986_design)

		# Now the contrast matrix is combined with the per-probeset linear model fit.
		huvec_fits <- contrasts.fit(fit, contrast.matrix)
		huvec_ebFit <- eBayes(huvec_fits)

		########
		# GET THE ELBOW LIMIT (this function)
		########
		get_elbow_limma(huvec_ebFit)

ELBOW documentation built on May 31, 2017, 3:48 p.m.

Search within the ELBOW package
Search all R packages, documentation and source code