R: Examples for 'survey::regTermTest'

### ** Examples

 data(esoph)
 model1 <- glm(cbind(ncases, ncontrols) ~ agegp + tobgp *
     alcgp, data = esoph, family = binomial())
 anova(model1)

Analysis of Deviance Table

Model: binomial, link: logit

Response: cbind(ncases, ncontrols)

Terms added sequentially (first to last)


            Df Deviance Resid. Df Resid. Dev
NULL                           87     367.95
agegp        5  121.045        82     246.91
tobgp        3   36.639        79     210.27
alcgp        3  127.933        76      82.34
tobgp:alcgp  9    5.451        67      76.89

 regTermTest(model1,"tobgp")

Wald test for tobgp
 in glm(formula = cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp, 
    family = binomial(), data = esoph)
F =  6.080611  on  3  and  67  df: p= 0.0010073

 regTermTest(model1,"tobgp:alcgp")

Wald test for tobgp:alcgp
 in glm(formula = cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp, 
    family = binomial(), data = esoph)
F =  0.5880663  on  9  and  67  df: p= 0.80234

 regTermTest(model1, ~alcgp+tobgp:alcgp)

Wald test for alcgp alcgp:tobgp
 in glm(formula = cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp, 
    family = binomial(), data = esoph)
F =  8.358247  on  12  and  67  df: p= 2.2352e-09

 data(api)
 dclus2<-svydesign(id=~dnum+snum, weights=~pw, data=apiclus2)
 model2<-svyglm(I(sch.wide=="Yes")~ell+meals+mobility, design=dclus2, family=quasibinomial())
 regTermTest(model2, ~ell)

Wald test for ell
 in svyglm(formula = I(sch.wide == "Yes") ~ ell + meals + mobility, 
    design = dclus2, family = quasibinomial())
F =  6.055448  on  1  and  36  df: p= 0.018792

 regTermTest(model2, ~ell,df=NULL)

Wald test for ell
 in svyglm(formula = I(sch.wide == "Yes") ~ ell + meals + mobility, 
    design = dclus2, family = quasibinomial())
F =  6.055448  on  1  and  36  df: p= 0.018792

 regTermTest(model2, ~ell, method="LRT", df=Inf)

Working (Rao-Scott) LRT for ell
 in svyglm(formula = I(sch.wide == "Yes") ~ ell + meals + mobility, 
    design = dclus2, family = quasibinomial())
Working 2logLR =  6.781297 p= 0.0096772 
df=1

 regTermTest(model2, ~ell+meals, method="LRT", df=NULL)

Working (Rao-Scott+F) LRT for ell meals
 in svyglm(formula = I(sch.wide == "Yes") ~ ell + meals + mobility, 
    design = dclus2, family = quasibinomial())
Working 2logLR =  4.692659 p= 0.11957 
(scale factors:  1.6 0.37 );  denominator df= 36

 regTermTest(model2, ~ell+meals, method="WorkingWald", df=NULL)

Working (Rao-Scott+F)  for ell meals
 in svyglm(formula = I(sch.wide == "Yes") ~ ell + meals + mobility, 
    design = dclus2, family = quasibinomial())
Working Wald statistic =  4.196737 p= 0.14263 
(scale factors:  1.6 0.37 );  denominator df= 36

Examples for 'survey::regTermTest'

Wald test for a term in a regression model