Modified gee function to output extra objects for saws

Share:

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

1
2
3
4
5
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

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