| pggm | R Documentation |
Profile boosting for Gaussian graphical model.
pggm(
S,
nObs,
maxK = floor(min(nObs - 1, NROW(S) - 1, 50)),
digits = 8,
verbose = FALSE
)
S |
Covariance matrix. |
nObs |
Number of observations. |
maxK |
Maximum number of identified edges. |
digits |
Integer indicating the number of decimal places or significant digits to be used. |
verbose |
Print the procedure path? |
Index set of identified features.
library(MASS)
library(Matrix)
set.seed(2025)
n <- 1000
p <- 10
Omega <- Diagonal(p)
diag(Omega[1:4, 2:5]) <- diag(Omega[2:5, 1:4]) <- 0.5
Sigma <- chol2inv(chol(Omega))
X <- mvrnorm(n, rep(0, p), Sigma, empirical=TRUE)
S <- cov(X)
system.time( egg <- pggm(S, n) )
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.