View source: R/cglasso_S3methods.R
summary.cglasso | R Documentation |
‘summary
’ produces a summary of the models fitted using cglasso
or cggm
.
## S3 method for class 'cglasso' summary(object, GoF = AIC, print.all = TRUE, digits = 3L, ...)
object |
an R object inheriting class |
GoF |
a valid goodness-of-fit function, such as |
print.all |
logical. Should all summary statistics be printed? |
digits |
the minimum number of significant digits to be used: see ‘print.default’. |
... |
the penalty parameter passed to the goodness-of-fit function (argument ‘ |
The function summary.cglasso
computes the summary statistics needed to evaluate the goodness-of-fit of the models fitted by cglasso
or cggm
.
Model evaluation can be made in two ways. The easiest way is to use a valid goodness-of-fit function, such as AIC.cglasso
or BIC.cglasso
. In this case, further arguments can be passed to these functions by the argument ‘...’. The second way consists in passing the output of a goodness-of-fit function, that is, an R object of class ‘GoF
’. Usually, this approach is preferable when the computation of the chosen goodness-of-fit measure is time-consuming, such as when the sample size is small relative to the number of parameters and the AIC
or BIC
functions are used to evaluate a sequence of fitted models. In these cases, we suggest the computation of several measures of goodness-of-fit in a preliminary step and then the use of summary.cglasso
in a subsequent step to evaluate the sequence of fitted models.
To improve the readability of the results, the output is divided into two sections.
The first section is structured in a sequence of tables showing, for each combination of the two tuning parameters ‘lambda
’ and ‘rho
’, the number of estimated non-zero regression coefficients (‘df.B
’), the number of estimated non-zero partial correlation coefficients (‘df.Tht
’), the degrees-of-freedom (‘df
’) and the value of the goodness-of-fit measure used to evaluate the fitted models. To help the user with the identification of the optimal fitted model, the last column of each table reports the ranking of the models (where the optimal model is marked with the symbol ‘<-
’).
The second section reports the summary statistics of the selected optimal model.
A named list with the following elements is silently returned:
table |
data.frame containing the summary statistics used to evaluate the sequence of fitted models. |
lambda.id |
position of the optimal lambda-value identified by the chosen goodness-of-fit function. |
rho.id |
position of the optimal rho-value identified by the chosen goodness-of-fit function. |
Luigi Augugliaro (luigi.augugliaro@unipa.it)
cglasso
, cggm
, AIC.cglasso
and BIC.cglasso
.
set.seed(123) # Y ~ N(0, Sigma) and probability of left/right censored values equal to 0.05 n <- 100L p <- 3L rho <- 0.3 Sigma <- outer(1L:p, 1L:p, function(i, j) rho^abs(i - j)) Z <- rcggm(n = n, Sigma = Sigma, probl = 0.05, probr = 0.05) out <- cglasso(. ~ ., data = Z) summary(out) # models evaluation by AIC summary(out, GoF = BIC) # models evaluation by BIC summary(out, GoF = BIC, mle = TRUE, g = 0.5) # models evaluation by eBIC # Y ~ N(b0 + XB, Sigma) and probability of left/right censored values equal to 0.05 n <- 100L p <- 3L q <- 2 b0 <- runif(p) B <- matrix(runif(q * p), nrow = q, ncol = p) X <- matrix(rnorm(n * q), nrow = n, ncol = q) rho <- 0.3 Sigma <- outer(1L:p, 1L:p, function(i, j) rho^abs(i - j)) Z <- rcggm(n = n, b0 = b0, X = X, B = B, Sigma = Sigma, probl = 0.05, probr = 0.05) out <- cglasso(. ~ ., data = Z, lambda = 0.01) summary(out) # models evaluation by AIC summary(out, GoF = BIC) # models evaluation by BIC summary(out, GoF = BIC, mle = TRUE, g = 0.5) # models evaluation by eBIC out <- cglasso(. ~ ., data = Z, rho = 0.01) summary(out) # models evaluation by AIC summary(out, GoF = BIC) # models evaluation by BIC summary(out, GoF = BIC, mle = TRUE, g = 0.5) # models evaluation by eBIC out <- cglasso(. ~ ., data = Z) summary(out) # models evaluation by AIC summary(out, GoF = BIC) # models evaluation by BIC summary(out, GoF = BIC, mle = TRUE, g = 0.5) # models evaluation by eBIC
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.