mgee: Modified gee function to output extra objects for saws

In saws: Small-Sample Adjustments for Wald Tests Using Sandwich Estimators

Modified gee function to output extra objects for saws

Description

This function calls the gee function from the gee package, then applies the geeUOmega function to it to create a score matrix (i.e., estimating equation) and information array (i.e., minus the derivative of the estimating equation). Since the mgee function just calls the gee function all help for gee applies to mgee.

Usage

mgee(formula = formula(data), id = id, data = parent.frame(),
    subset, na.action, R = NULL, b = NULL, tol = 0.001, 
    maxiter = 25, family = gaussian, corstr = "independence", 
    Mv = 1, silent = TRUE, contrasts = NULL, scale.fix = FALSE,
    scale.value = 1, v4.4compat = FALSE)

Arguments

formula	see gee help
id	see gee help
data	see gee help
subset	see gee help
na.action	see gee help
R	see gee help
b	see gee help
tol	see gee help
maxiter	see gee help
family	see gee help
corstr	see gee help
Mv	see gee help
silent	see gee help
contrasts	see gee help
scale.fix	see gee help
scale.value	see gee help
v4.4compat	see gee help

Value

A gee object with two extra elements to the list, u and omega (see saws).

Note

You can alternatively take the output from gee and apply the geeUOmega function. But see the warning for that function.

Author(s)

last few lines by M.P. Fay, for the rest see gee package DESCRIPTION

See Also

gee, geeUOmega

Examples

##  example from gee help
data(warpbreaks)
mout<-mgee(breaks~tension,id=wool,data=warpbreaks,corstr="exchangeable")
saws(mout)

saws documentation built on June 24, 2022, 1:07 a.m.