ggm_compare: Compare Gaussian Graphical Models

In donaldRwilliams/GGMnonreg: Non-Regularized Gaussian Graphical Models

Description

Establish whether each of the corresponding edges are significantly different in two groups

Usage

ggm_compare(Yg1, Yg2, method = "spearman", alpha = 0.05)

Arguments

Yg1	The data matrix of dimensions n (observations) by p (nodes) for group one.
Yg2	The data matrix of dimensions n (observations) by p (nodes) for group two.
method	Character string. Which type of correlation coefficients to be computed. Options include "pearson" (default), "kendall", "spearman", and "polychoric".
alpha	The desired significance level (defaults to 0.05). Note that 1 - alpha corresponds to specificity.

Value

An object of class ggm_compare including:

• adj: Adjacency matrix, where a 1 indicates a difference.

• wadj: Weighted adjacency matrix (partial correlation differences that were significantly different)

• cis: Confidence intervals for the partial correlation differences.

Examples

# data

Yg1 <- na.omit(subset(bfi, gender == 1)[,1:10])
Yg2 <- na.omit(subset(bfi, gender == 2)[,1:10])

# compare relations
fit <- ggm_compare(Yg1, Yg2)

donaldRwilliams/GGMnonreg documentation built on Nov. 13, 2021, 9:57 a.m.