tests/regTermTest-missing.R
In bschneidr/fastsurvey: Analysis of Complex Survey Samples
library(survey)
data(nhanes)
design <- svydesign(id=~SDMVPSU, strata=~SDMVSTRA, weights=~WTMEC2YR, nest=TRUE,data=nhanes)
a<-svyglm(formula = I(race == 1) ~ HI_CHOL + agecat + RIAGENDR, design = subset(design,!is.na(HI_CHOL)), family=quasibinomial)
b<-svyglm(formula = I(race == 1) ~ HI_CHOL + agecat + RIAGENDR, design =design , family=quasibinomial)
ta<-regTermTest(a, ~HI_CHOL)
tb<-regTermTest(b, ~HI_CHOL)
stopifnot(isTRUE(all.equal(ta$chisq, tb$chisq)))
stopifnot(isTRUE(all.equal(ta$lambda, tb$lambda)))
bschneidr/fastsurvey documentation built on March 13, 2024, 11:12 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
library(survey)
data(nhanes)
design <- svydesign(id=~SDMVPSU, strata=~SDMVSTRA, weights=~WTMEC2YR, nest=TRUE,data=nhanes)
a<-svyglm(formula = I(race == 1) ~ HI_CHOL + agecat + RIAGENDR, design = subset(design,!is.na(HI_CHOL)), family=quasibinomial)
b<-svyglm(formula = I(race == 1) ~ HI_CHOL + agecat + RIAGENDR, design =design , family=quasibinomial)
ta<-regTermTest(a, ~HI_CHOL)
tb<-regTermTest(b, ~HI_CHOL)
stopifnot(isTRUE(all.equal(ta$chisq, tb$chisq)))
stopifnot(isTRUE(all.equal(ta$lambda, tb$lambda)))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.