MirSEA: Identify dysregulated pathways based on microRNA (miRNA) set...

View source: R/MirSEA.R

MirSEAR Documentation

Identify dysregulated pathways based on microRNA (miRNA) set enrichment analysis

Description

This function propose a novel method of miRNA set enrichment analysis(MiRSEA)to identify the dysregulated pathways by calculating the enrichment score of miRNA set which co-regulate a biological pathway(or prior gene set)

Usage

MirSEA(input.ds, input.cls, p_value,p2miR,
reshuffling.type = "miR.labels", nperm = 1000, 
weighted.score.type = 1, ms.size.threshold.min = 10, 
ms.size.threshold.max = 500)

Arguments

input.ds

Input miRNA expression Affymetrix dataset file in GCT format

input.cls

Input class vector (phenotype) file in CLS format

p_value

A weighting matrix of p value of the hypergeometric. (rows are pathway ,cols are microRNAs(miRNAs))

p2miR

pathway-miRNA correlation(pmSET) profile

reshuffling.type

Type of permutation reshuffling: "sample.labels" or "miR.labels" (default: "miR.labels")

nperm

Number of random permutations (default: 1000)

weighted.score.type

Enrichment correlation based weighting:When weight= 0, ES reduces to the standard Kolmogorov-Smirnov statistic,when weight=1, we are weighting the miRNAs by their dw-score normalized by the sum of the dw-scores over all of the miRNAs in the miRNA set,when weight=2,it represent over weight (default: 1)

ms.size.threshold.min

Minimum size (in miRNAs) for database miRNA sets to be considered (default: 10)

ms.size.threshold.max

Maximum size (in miRNAs) for database miRNA sets to be considered (default: 500)

Details

MiRSEA integrates pathway (e.g.the strength of the pathway regulated by miRNAs.) and differential expression among miRNAs in identifying dysregulated pathways.MiRSEA can order pathway by the enrichment score of miRNA set, which is co-regulated by a miRNA set.

Value

report.phen1

It is the summary of the result of the up regulated pathway

report.phen2

It is the summary of the result of the down regulated pathway.Each rows of the dataframe represents a pathway. Its columns include "Pathway Name", "SIZE", "Pathway Source", "Pathway Enrichment Score", "NOM p-val", "FDR q-val", "Tag percentage"(Percent of miRNA set before running enrichment peak),"MiR percentage"(Percent of miRNA list before running enrichment peak),"Signal strength" (enrichment signal strength).

Author(s)

Junwei Hanhanjunwei1981@163.com,Siyao Liu liusiyao29@163.com

References

Subramanian A, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America 102(43):15545-15550.

Lu M, Shi B, Wang J, Cao Q, & Cui Q (2010) TAM: a method for enrichment and depletion analysis of a microRNA category in a list of microRNAs. BMC bioinformatics 11:419.

See Also

EnrichmentScore, EnrichmentScore2,S2N,Corrp2miRfile

Examples

## Not run: 
#get example of expression data
#input.ds <- readLines("F:/lsy/xin data/GSE36915.gct")
#input.cls <- readLines("F:/lsy/xin data/GSE36915.cls")
input.ds <- GetExampleData("dataset")
input.cls <- GetExampleData("class.labels")

#get example of p value matrix
p_value <- GetExampleData("p_value")
#get example of correlation profile
p2miR <- GetExampleData("p2miR")

#identify dysregulated pathways by using the function MirSEA
MirSEAresult <- MirSEA(input.ds,input.cls,p_value,p2miR,
reshuffling.type = "miR.labels", nperm = 1000, 
weighted.score.type = 1, ms.size.threshold.min = 10, 
ms.size.threshold.max = 500)
#print the summary results of pathways to screen
 
summaryResult1 <- MirSEAresult$report.phen1
summaryResult1[1:5,]
summaryResult2 <- MirSEAresult$report.phen2
summaryResult2[1:5,]
 
#write the summary results of pathways to tab delimited file.
write.table(summaryResult1,file="summaryResult1.txt",sep="\t",row.names=FALSE)

write.table(summaryResult2,file="summaryResult2.txt",sep="\t",row.names=FALSE)

## End(Not run)

MiRSEA documentation built on May 29, 2024, 4:36 a.m.