Returns the consensus network pathways from a multiple alignment.
Usage
consensusPathways(object)
# S4 method for class 'ConsortiumMetabolismAlignment'
consensusPathways(object)Arguments
- object
A ConsortiumMetabolismAlignment object of type
"multiple".
Methods (by class)
consensusPathways(ConsortiumMetabolismAlignment): Consensus pathways from a multiple ConsortiumMetabolismAlignment
Examples
# \donttest{
cm1 <- synCM("comm_1", n_species = 3, max_met = 5)
cm2 <- synCM("comm_2", n_species = 4, max_met = 6)
cms <- ConsortiumMetabolismSet(cm1, cm2, name = "test")
#>
#> ── Creating CMS "test" ─────────────────────────────────────────────────────────
#> ℹ Validating 2 <ConsortiumMetabolism> objects
#> ✔ Validating 2 <ConsortiumMetabolism> objects [19ms]
#>
#> ℹ Collecting metabolites from 2 consortia
#> ✔ Collecting metabolites from 2 consortia [30ms]
#>
#> ℹ Re-indexing 6 unique metabolites
#> ✔ Re-indexing 6 unique metabolites [26ms]
#>
#> ℹ Expanding 2 binary matrices to 6-dimensional space
#> ✔ Expanding 2 binary matrices to 6-dimensional space [24ms]
#>
#> ℹ Computing 6 x 6 levels matrix
#> ✔ Computing 6 x 6 levels matrix [25ms]
#>
#> ℹ Computing pairwise overlap (1 pairs via crossprod)
#> ✔ Computing pairwise overlap (1 pairs via crossprod) [23ms]
#>
#> ℹ Assembling pathway data from 2 consortia
#> ✔ Assembling pathway data from 2 consortia [30ms]
#>
#> ℹ Building dendrogram from 2 x 2 dissimilarity matrix
#> ✔ Building dendrogram from 2 x 2 dissimilarity matrix [22ms]
#>
#> ℹ Extracting dendrogram node positions
#> ✔ Extracting dendrogram node positions [24ms]
#>
#> ℹ Collecting 2 consortium graphs
#> ✔ CMS "test" created: 2 consortia, 6 metabolites (0.2s)
#> ℹ Collecting 2 consortium graphs
#> ✔ Collecting 2 consortium graphs [84ms]
#>
cma <- align(cms)
#> Computing multiple alignment for 2 consortia using "FOS".
consensusPathways(cma)
#> consumed produced nConsortia proportion
#> 1 met2 met1 1 0.5
#> 2 met3 met1 1 0.5
#> 3 met4 met1 1 0.5
#> 4 met5 met1 1 0.5
#> 5 met1 met2 1 0.5
#> 6 met3 met2 1 0.5
#> 7 met5 met2 1 0.5
#> 8 met6 met2 1 0.5
#> 9 met2 met3 1 0.5
#> 10 met4 met3 1 0.5
#> 11 met1 met4 1 0.5
#> 12 met5 met4 1 0.5
#> 13 met6 met4 1 0.5
#> 14 met1 met5 1 0.5
#> 15 met2 met5 1 0.5
#> 16 met3 met5 1 0.5
#> 17 met4 met5 1 0.5
#> 18 met6 met5 1 0.5
# }