Drop All Possible Single Terms to a geem Model Using Wald or Score Test — drop1.geem • geeasy

Compute all the single terms in the scope argument that can dropped from the model, and compute a table of the corresponding Wald test statistics.

# S3 method for class 'geem'
drop1(
  object,
  scope,
  test = c("Wald", "none", "score", "sasscore"),
  method = c("robust", "naive", "sandwich"),
  ...
)

Arguments

object: a fitted object of class geese.
scope: a formula giving the terms to be considered for adding or dropping.
test: the type of test to include.
method: Indicates which method is used for computing the standard error. robust is the default and corresponds to the modified sandwich estimator. naive is the classical naive variance estimate. sandwich is an alias for robust.
...: other arguments. Not currently used

Value

An object of class anova summarizing the differences in fit between the models.

See also

drop1, geem

Author

Claus Ekstrom claus@ekstroem.dk

Examples


library(geeM)
#> Loading required package: Matrix
library(geepack)
data(ohio)
fit <- geem(resp ~ age + smoke + age:smoke, id=id, data=ohio,
            family="binomial", corstr="exch", scale.fix=TRUE)
drop1(fit)
#> Single term deletions
#> 
#> Model:
#> resp ~ age + smoke + age:smoke
#>           DF    Wald Pr(>Chi)
#> age:smoke  1 0.64379   0.4223