Computes Wald confidence intervals for model parameters using the normal approximation.
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)
confint(fit)
#> 2.5 % 97.5 %
#> (Intercept) -1.4390767 1.9492767
#> x1 -1.2809286 0.8405165
#> (phi) -0.9343775 0.1485488
confint(fit, model = "mean")
#> 2.5 % 97.5 %
#> (Intercept) -1.439077 1.9492767
#> x1 -1.280929 0.8405165
# }
