Variance-covariance matrix of estimated coefficients
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)
)
prep <- brs_prep(dat, ncuts = 100)
#> brs_prep: n = 20 | exact = 0, left = 1, right = 1, interval = 18
fit <- brs(y ~ x1, data = prep)
vcov(fit)
#> (Intercept) x1 (phi)
#> (Intercept) 0.74717129 -0.444368049 -0.010878707
#> x1 -0.44436805 0.292891282 0.009448033
#> (phi) -0.01087871 0.009448033 0.076320255
vcov(fit, model = "mean")
#> (Intercept) x1
#> (Intercept) 0.7471713 -0.4443680
#> x1 -0.4443680 0.2928913
# }
