getCluster {HGMND}

Get the cluster structure of the HGMND estimate

Description

After estimating the conditional dependence matrices of the multiple datasets using the HGMND method, the cluster structure can be revealed by comparison of these matrices.

Usage

  getCluster(est.HGMND, method = "F", tol = 1e-5)

Arguments

est.HGMND

a list, the result of the function HGMND with class "est.HGMND".

method

the method of evaluating the difference of two conditional dependence matrices. The function norm from the base package is used to calculate the matrix norm of the element-wise difference of two matrices. It must be chosen from "O", "I", "F", "M", "2", corresponding to the same settings in the function norm. Default to "F", the Frobenius norm.

tol

tolerance in evaluating the difference of two conditional dependence matrices. If the calculated difference is no larger than tol, they are regarded as in one cluster. Default to 1e-5.

Value

the function getCluster returns the clustering structure of the multiple conditional dependence matrices.

mat.comapre

a matrix of 0 or 1. If the element on the ith row and jth column of the matrix is 1, the ith and the jth conditional dependence matrices are in the same cluster, 0 otherwise.

est.cluster

a vector with length same as the number of conditional dependence matrices indicating the cluster label of each matrix.

Examples

# This is an example of HGMND with simulated data
data(HGMND_SimuData)
h              <- genscore::get_h_hp("mcp", 1, 5)
HGMND_SimuData <- lapply(HGMND_SimuData, function(x) scale(x, center = FALSE))
mat.chain      <- diag(length(HGMND_SimuData))
diag(mat.chain[-nrow(mat.chain), -1]) <- 1

result <- HGMND(x        = HGMND_SimuData,
                setting  = "gaussian",
                h        = h,
                centered = FALSE,
                mat.adj  = mat.chain,
                lambda1  = 0.086,
                lambda2  = 3.6,
                gamma    = 1,
                tol      = 1e-3,
                silent  = TRUE)
Theta       <- result[["Theta"]]
res.cluster <- getCluster(result)