tests/na_action.R

## from Terry Therneau
library(survey)
load("naa.rda")

fit1e <- svyglm( pseudo ~ age34 + ccr5 + factor(times),  design= adata.s,na.action=na.exclude)
fit1o <- svyglm( pseudo ~ age34 + ccr5 + factor(times),  design= adata.s)
all.equal(coef(fit1e),coef(fit1o))
all.equal(vcov(fit1e),vcov(fit1o))

Try the survey package in your browser

Any scripts or data that you put into this service are public.

survey documentation built on May 3, 2023, 9:12 a.m.