getMaxStats {BioTIP} | R Documentation

## Get the cluster index and network nodes of biomodule

### Description

This function retrives the cluster index and network-node ids for the identified biomodule (that shows the maximum MCI score) at each state in the study.

### Usage

    getMaxStats(membersL, idx)

### Arguments

`membersL` | A two-layer nested list of character or numeric values, any one out of the five elements outputed by the function `getMCI`.
`idx` | A vector of integers that are cluster ids of the biomodule (the module with the highest MCI score) per state. This is the first element of the result from `getMaxMCImember`

### Value

A list describing the biomodule of each state, corresponding to one of the five elements (members, MCI, Sd, PCC, and PCCo) outputted by the function `getMCI`. The calss of the vector depends on the class of the input parameter `membersL`.

### Author(s)

Zhezhen Wang zhezhen@uchicago.edu

### See Also

`getMCI`

### Examples

    test = list('state1' = matrix(sample(1:10,6),4,3),'state2' = matrix(sample(1:10,6),4,3),'state3' = matrix(sample(1:10,6),4,3))
    # assign colnames and rownames to the matrix
    for(i in names(test)){
     colnames(test[[i]]) = 1:3
     row.names(test[[i]]) = c('g1','g2','g3','g4')
    }

    cluster = list(c(1,2,2,1),c(1,2,3,1),c(2,2,1,1))
    names(cluster) = names(test)
    for(i in names(cluster)){
     names(cluster[[i]]) = c('g1','g2','g3','g4')
    }
    membersL_noweight <- getMCI(cluster,test)
    idx = c(1,2,1)
    names(idx) = names(membersL_noweight[['sd']])
    selectedSD = getMaxStats(membersL_noweight[['sd']],idx)

[Package BioTIP version 1.0.0 Index]