bread.ols: Getting the bread for the vcovHC

Description Usage Arguments Value Examples

View source: R/rms_SandwichAddon.R

Description

The original bread.lm uses the summary.lm function it seems like a quick fix and I've therefore created the original bread definition: $(X'X)^-1$

Usage

1
2
## S3 method for class 'ols'
bread(x, ...)

Arguments

x

The ols model fit

...

arguments passed to methods.

Value

matrix The bread for the sandwich vcovHC function

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
# Generate some data
n <- 500
x1 <- runif(n)*2
x2 <- runif(n)
y <- x1^3 + x2 + rnorm(n)

library(rms)
library(sandwich)
dd <- datadist(x1, x2, y)
org.op <- options(datadist = "dd")

# Main function
f    <- ols(y ~ rcs(x1, 3) + x2)

# Check the bread
bread(f)
# Check the HC-matrix
vcovHC(f, type="HC4m")
# Adjust the model so that it uses the HC4m variance
f_rob <- robcov_alt(f, type="HC4m")
# Get the new HC4m-matrix
# - this function just returns the f_rob$var matrix
vcov(f_rob)
# Now check the confidence interval for the function
confint(f_rob)

options(org.op)

gforge/Greg documentation built on Nov. 16, 2019, 3:43 p.m.