Print a fitted brsmm model
Examples
# \donttest{
dat <- data.frame(
y = c(
0, 5, 20, 50, 75, 90, 100, 30, 60, 45,
10, 40, 55, 70, 85, 25, 35, 65, 80, 15
),
x1 = rep(c(1, 2), 10),
id = factor(rep(1:4, each = 5))
)
prep <- brs_prep(dat, ncuts = 100)
#> brs_prep: n = 20 | exact = 0, left = 1, right = 1, interval = 18
fit <- brsmm(y ~ x1, random = ~ 1 | id, data = prep)
print(fit)
#>
#> Call:
#> brsmm(formula = y ~ x1, random = ~1 | id, data = prep)
#>
#> Coefficients (mean model with logit link):
#> (Intercept) x1
#> 0.4211 -0.3373
#>
#> Phi coefficients (precision model with logit link):
#> (phi)_(Intercept)
#> -0.5805
#>
#> Random-effects parameters:
#> (re_chol_logsd)_(Intercept)|id
#> -0.6277
#>
#> Random SD: 0.5338
#> ---
#> Mixed beta interval model (Laplace)
#> Observations: 20 | Groups: 4
#> Log-likelihood: -92.1831
#> Convergence code: 0
# }
