powered by
simex
mcsimex
extract.covmat(model)
lm
glm
gam
nls
lme
nlme
x <- rnorm(100,0,100) u <- rnorm(100,0,25) w <- x+u y <- x +rnorm(100,0,9) model <- lm(y~w, x = TRUE, y =TRUE) extract.covmat(model) summary(model)$cov.unscaled*summary(model)$sigma^2
Run the code above in your browser using DataLab