Compute Bayesian information criterion for a mixedLSR model
Arguments
- a
A list of coefficient matrices.
- n
The sample size.
- llik
The log-likelihood of the model.
Examples
n <- 50
simulate <- simulate_lsr(n)
model <- mixed_lsr(simulate$x, simulate$y, k = 2, init_lambda = c(1,1), alt_iter = 0)
#> mixedLSR Start: 1
#> Computing Final Model...
#> Done!
bic_lsr(model$A, n = n, model$llik)
#> [1] 5686.224