BGGM (version 1.0.0)

select.ggm_compare_estimate: Edge Differences and (Practical) Equivalence Between GGMs

Description

Edge Differences and (Practical) Equivalence Between GGMs

Usage

# S3 method for ggm_compare_estimate
select(object, cred = 0.95, rope = NULL, prob = 0.95, ...)

Arguments

object

object of class ggm_compare_estimate

cred

credible interval width used for the decision rule

rope

region of practical equivalence

prob

posterior probability (see notes)

...

not currently used

Value

A list of class select.ggm_compare_estimate: rope = NULL:

• mat_adj adjacency matrix (one for each contrast)

• mat_pcor selected partial correlations (one for each contrast)

• call match.call()

• object object of class ggm_compare_estimate

• rope region of practical equivalence

• cred credible interval

• prob posterior probability

Examples

# NOT RUN {
# data
Y1 <- BGGM::bfi[1:500,1:5]
Y2 <- BGGM::bfi[501:1000, 1:5]

# fit model
fit <- ggm_compare_estimate(Y1, Y2)

# posterior summary of differences
summary(fit)

# select (threshold) with credible intervals
sel <- select(fit)

# summary
summary(sel)

# selected differences
sel$mat_pcor

# adjacency matrix
sel$mat_adj
# }