select.estimate: Graph Selection for 'estimate' Objects

In BGGM: Bayesian Gaussian Graphical Models

Graph Selection for estimate Objects

Description

Provides the selected graph based on credible intervals for the partial correlations that did not contain zero \insertCiteWilliams2019BGGM.

Usage

## S3 method for class 'estimate'
select(object, cred = 0.95, alternative = "two.sided", ...)

Arguments

object	An object of class estimate.default.
cred	Numeric. The credible interval width for selecting the graph (defaults to 0.95; must be between 0 and 1).
alternative	A character string specifying the alternative hypothesis. It must be one of "two.sided" (default), "greater" or "less". See note for futher details.
...	Currently ignored.

Details

This package was built for the social-behavioral sciences in particular. In these applications, there is strong theory that expects all effects to be positive. This is known as a "positive manifold" and this notion has a rich tradition in psychometrics. Hence, this can be incorporated into the graph with alternative = "greater". This results in the estimated structure including only positive edges.

Value

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

• pcor_adj Selected partial correlation matrix (weighted adjacency).

• adj Adjacency matrix for the selected edges

• object An object of class estimate (the fitted model).

References

\insertAllCited

See Also

estimate and ggm_compare_estimate for several examples.

Examples

# note: iter = 250 for demonstrative purposes

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

# estimate
fit <- estimate(Y, iter = 250,
                progress = FALSE)


# select edge set
E <- select(fit)

BGGM documentation built on Sept. 11, 2024, 5:19 p.m.