R: 3M Markov mixture model for clustering pathways
pathCluster
R Documentation
3M Markov mixture model for clustering pathways
Description
3M Markov mixture model for clustering pathways
Usage
pathCluster(ybinpaths, M, iter = 1000)
Arguments
ybinpaths
The training paths computed by pathsToBinary.
M
The number of clusters.
iter
The maximum number of EM iterations.
Value
A list with the following items:
h
The posterior probabilities that each path belongs to each cluster.
labels
The cluster membership labels.
theta
The probabilities of each gene for each cluster.
proportions
The mixing proportions of each path.
likelihood
The likelihood convergence history.
params
The specific parameters used.
Author(s)
Ichigaku Takigawa
Timothy Hancock
References
Mamitsuka, H., Okuno, Y., and Yamaguchi, A. 2003. Mining biologically active patterns in
metabolic pathways using microarray expression profiles. SIGKDD Explor. News l. 5, 2 (Dec. 2003), 113-121.
## Prepare a weighted reaction network.
## Conver a metabolic network to a reaction network.
data(ex_sbml) # bipartite metabolic network of Carbohydrate metabolism.
rgraph <- makeReactionNetwork(ex_sbml, simplify=TRUE)
## Assign edge weights based on Affymetrix attributes and microarray dataset.
# Calculate Pearson's correlation.
data(ex_microarray) # Part of ALL dataset.
rgraph <- assignEdgeWeights(microarray = ex_microarray, graph = rgraph,
weight.method = "cor", use.attr="miriam.uniprot", bootstrap = FALSE)
## Get ranked paths using probabilistic shortest paths.
ranked.p <- pathRanker(rgraph, method="prob.shortest.path",
K=20, minPathSize=8)
## Convert paths to binary matrix.
ybinpaths <- pathsToBinary(ranked.p)
p.cluster <- pathCluster(ybinpaths, M=2)
plotClusters(ybinpaths, p.cluster)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(NetPathMiner)
Loading required package: igraph
Attaching package: 'igraph'
The following objects are masked from 'package:stats':
decompose, spectrum
The following object is masked from 'package:base':
union
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/NetPathMiner/pathCluster.Rd_%03d_medium.png", width=480, height=480)
> ### Name: pathCluster
> ### Title: 3M Markov mixture model for clustering pathways
> ### Aliases: pathCluster
>
> ### ** Examples
>
> ## Prepare a weighted reaction network.
> ## Conver a metabolic network to a reaction network.
> data(ex_sbml) # bipartite metabolic network of Carbohydrate metabolism.
> rgraph <- makeReactionNetwork(ex_sbml, simplify=TRUE)
>
> ## Assign edge weights based on Affymetrix attributes and microarray dataset.
> # Calculate Pearson's correlation.
> data(ex_microarray) # Part of ALL dataset.
> rgraph <- assignEdgeWeights(microarray = ex_microarray, graph = rgraph,
+ weight.method = "cor", use.attr="miriam.uniprot", bootstrap = FALSE)
100 genes were present in the microarray, but not represented in the network.
55 genes were couldn't be found in microarray.
Assigning edge weights.
>
> ## Get ranked paths using probabilistic shortest paths.
> ranked.p <- pathRanker(rgraph, method="prob.shortest.path",
+ K=20, minPathSize=8)
Extracting the 20 most probable paths.
>
> ## Convert paths to binary matrix.
> ybinpaths <- pathsToBinary(ranked.p)
> p.cluster <- pathCluster(ybinpaths, M=2)
> plotClusters(ybinpaths, p.cluster)
>
>
>
>
>
> dev.off()
null device
1
>