BGGM (version 2.1.5)

select.explore: Graph selection for explore Objects

Description

Provides the selected graph based on the Bayes factor Williams2019_bfBGGM.

Usage

# S3 method for explore
select(object, BF_cut = 3, alternative = "two.sided", ...)

Value

The returned object of class select.explore contains a lot of information that is used for printing and plotting the results. For users of BGGM, the following are the useful objects:

alternative = "two.sided"

• pcor_mat_zero Selected partial correlation matrix (weighted adjacency).

• pcor_mat Partial correlation matrix (posterior mean).

• Adj_10 Adjacency matrix for the selected edges.

• Adj_01 Adjacency matrix for which there was evidence for the null hypothesis.

alternative = "greater" and "less"

• pcor_mat_zero Selected partial correlation matrix (weighted adjacency).

• pcor_mat Partial correlation matrix (posterior mean).

• Adj_20 Adjacency matrix for the selected edges.

• Adj_02 Adjacency matrix for which there was evidence for the null hypothesis (see note).

alternative = "exhaustive"

• post_prob A data frame that included the posterior hypothesis probabilities.

• neg_mat Adjacency matrix for which there was evidence for negative edges.

• pos_mat Adjacency matrix for which there was evidence for positive edges.

• neg_mat Adjacency matrix for which there was evidence for the null hypothesis (see note).

• pcor_mat Partial correlation matrix (posterior mean). The weighted adjacency matrices can be computed by multiplying pcor_mat with an adjacency matrix.

Arguments

object

An object of class explore.default

BF_cut

Numeric. Threshold for including an edge (defaults to 3).

alternative

A character string specifying the alternative hypothesis. It must be one of "two.sided" (default), "greater", "less", or "exhaustive". See note for further details.

...

Currently ignored.

Details

Exhaustive provides the posterior hypothesis probabilities for a positive, negative, or null relation @see Table 3 in @Williams2019_bfBGGM.

References

See Also

explore and ggm_compare_explore for several examples.

Examples

# \donttest{
#################
### example 1 ###
#################

#  data
Y <- bfi[,1:10]

# fit model
fit <- explore(Y, progress = FALSE)

# edge set
E <- select(fit,
            alternative = "exhaustive")

# }