findSigGenes: Find genes significantly differentially expressed between two...

Description Usage Arguments Details Value Author(s) References Examples

View source: R/findSigGenes.r

Description

For a given expression profile of two conditions, find the genes differencially expressed using T-test, fold change or SAM algorithm.

Usage

1
2
findSigGenes(Expr, Label, Method = "tTest", Directed = TRUE, 
	FdrCut = 0.01, FDCut = 1)

Arguments

Expr

A data frame, the expression profile to find differentially expressed genes, the rownames should be the ID of genes.

Label

A vector of 0/1s, indicating the class of samples in the expression profile, 0 represents case, 1 represents control.

Method

A string, specifying the method to calculate the differentially expressed genes, should be one of the "tTest"or"foldChange".

Directed

Logical, if the the up or down regulated set should be distinguished.

FdrCut

Numeric, the fdr cutoff for T test, can be ignored if not using t-test.

FDCut

Numeric, the cutoff for fold change, can be ignored if not using fold change.

Details

For a given expression profile of two conditions, lncPath package provide two method to find differentially expressed genes: t-text and fold change. The row of the expression profile should be gene IDs and the column of the expression profile should be names of samples. Samples should be under two conditions and the label should be given as 0 and 1. For t-test, fold change and SAM, different threshold can be set for significant differentially expressed genes.

Value

A vector of strings, the IDs of differentially expressed genes.

Author(s)

Junwei Han <[email protected]>, Zeguo Sun <[email protected]>

References

Subramanian, A., Tamayo, P., Mootha, V.K., Mukherjee, S., Ebert, B.L., Gillette, M.A., Paulovich, A., Pomeroy, S.L., Golub, T.R., Lander, E.S. et al. (2005) Gene set enrichment analysis: a knowledgebased approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A, 102, 15545-15550.

Examples

1
2
3
4
5
6
7
8
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.
Profile <- getExampleData("Profile")
Labels <- getExampleData("Labels")

SigGenes <- findSigGenes(Profile, Labels)
head(SigGenes)

LncPath documentation built on Sept. 26, 2018, 5:04 p.m.