Learn R Programming

directPA (version 1.5.1)

pathwayStats: Pathway Level Statistics

Description

Takes a vector of statistics with each element corresponds to a gene or phosphorylation site, and calculates a combined statistics for those that belong to the same pathway or kinase.

Usage

pathwayStats(PGs, T, minSize=5, method="Stouffer")

Value

a doublet corresponding to the enrichment after integration across all genes or substrates that belong to the same pathway or kinase, and the size of the mapped genes or substrates to that pathway or kinase.

Arguments

PGs

an array of names indicating genes or substrates that belong to a given pathway or kinase.

T

a vector of statistics (z-scores converted) with each element correspond to a gene or phosphorylation site that belong to the same pathway or kinase.

minSize

the size of annotation groups to be considered for calculating enrichment. Groups that are smaller than the minSize will be removed from the analysis.

method

the p-value integration method for combining accross multiple treatments. Available methods are Stouffer, OSP, Fisher, and maxP. The default method is Stouffer.

Examples

Run this code
# load the example data
data(PM)

# load pathway annotations
data(Pathways)

# convert statistics into z-scores
PM.zscores <- apply(PM, 2, function(x){qnorm(rank(x)/(nrow(PM)+1))})

# Rotate the matrix by contrast 1, -1, -1 (i.e. up-regulation, down-regulation, dow-regulation).
PM.rotated <- rotate3d(PM.zscores, contrast = c(1, -1, -1))

# combine rotated statistics across treatments
gene.pvalues <- apply(PM.rotated, 1, geneStats)

# compute statistics for all reactome pathways
gene.zscores <- qnorm(gene.pvalues, lower.tail = FALSE)
gst <- t(sapply(Pathways.reactome, pathwayStats, gene.zscores))

Run the code above in your browser using DataLab