regTermTest {survey} | R Documentation |

## Wald test for a term in a regression model

### Description

Provides Wald test and working likelihood ratio (Rao-Scott) test of the
hypothesis that all coefficients associated with a particular
regression term are zero (or have some other specified
values). Particularly useful as a substitute for `anova`

when not fitting by maximum likelihood. The Wald tests use a chisquared or F distribution, the LRT uses a linear combination of chisquared distributions as in `pchisqsum`

.

### Usage

regTermTest(model, test.terms, null=NULL,df=Inf, method=c("Wald","LRT"))

### Arguments

`model` |
A model object with `coef` and `vcov` methods |

`test.terms` |
Character string or one-sided formula giving name of
term or terms to test |

`null` |
Null hypothesis values for parameters. Default is zeros |

`df` |
Denominator degrees of freedom for an F test. If
`NULL` these are estimated from the model |

`method` |
If `"Wald"` , the Wald-type test; if `"LRT"`
the Rao-Scott test based on the estimated log likelihood ratio |

### Value

An object of class `regTermTest`

or `regTermTestLRT`

.

### References

Rao, JNK, Scott, AJ (1984) "On Chi-squared Tests For Multiway Contingency Tables with Proportions Estimated From Survey Data" Annals of Statistics 12:46-60.

### See Also

`anova`

, `vcov`

, `contrasts`

,`pchisqsum`

### Examples

data(esoph)
model1 <- glm(cbind(ncases, ncontrols) ~ agegp + tobgp *
alcgp, data = esoph, family = binomial())
anova(model1)
regTermTest(model1,"tobgp")
regTermTest(model1,"tobgp:alcgp")
regTermTest(model1, ~alcgp+tobgp:alcgp)
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)
regTermTest(model2, ~ell,df=NULL)
regTermTest(model2, ~ell, method="LRT")
regTermTest(model2, ~ell+meals, method="LRT")

[Package

*survey* version 3.18

Index]