Skip to contents

Variance-covariance matrix of estimated coefficients

Usage

# S3 method for class 'brs'
vcov(object, model = c("full", "mean", "precision"), ...)

Arguments

object

A fitted "betaregscale" object.

model

Character: which component ("full", "mean", or "precision").

...

Ignored.

Value

A square numeric matrix.

See also

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
# }