manuscript/comparemm.R

In mac-theobio/effects: isolated confidence intervals and bias-corrected predictor effects

library(shellpipes)
library(varpred)
loadEnvironments()

mm <- model.matrix(fit)
bb <- coef(fit)
v1 <- mm %*% bb |> as.vector()
v2 <- predict(fit)
dd <- data.frame(v1=v1, v2=v2)
mean(v2)
hhsize_mean
mean(df$hhsize)
age_mean <- mean(df$age)
v3 <- varpred(fit, "age", at=list(age=age_mean))
v3
v4 <- varpred(fit, "age", at=list(age=df$age))
v4$preds$fit |> mean()