TITLE(gee @@ Function to solve a Generalized Estimation Equation Model (version 3.5) )
Produces an object of class "gee" which is a generalized Estimation Equation fit of the data.
USAGE(
gee(formula, id,
    data, subset, na.action,
    offset=NA, R=NA, b=NA,
    tol=0.001, maxiter=as.integer(25),
    family = gaussian, corstr="independence",
    Mv=1, silent=T, contrasts=NULL,
    scale.fix = F, scale.value = 1)
)
ARGUMENTS(
ARG(formula@@)
a formula expression as for other regression models, of the form response ~ predictors. See the documentation of lm and formula for details.
ARG(id@@)
a vector which identifies the clusters.  The length of id should be the same as the number of observations.  Data are assumed to be sorted
so that observations on a cluster are contiguous rows for all entities
in the formula.
ARG(data@@)
  an optional data frame in which to interpret the variables occurring in the formula, along with the id, offset and n variables.
ARG(subset@@)
  expression saying which subset of the rows of the  data should  be  used in the fit.  This can be a logical vector (which is replicated to have length equal to the number of observations), or a numeric vector indicating which observation numbers are to be included, or a  character  vector of the row names to be included.  All observations are included by default.
ARG(na.action@@)
  a function to filter missing data.  For gee only na.omit should be used here.
ARG(offset@@)
a vector giving the offset.  If the length of offset is less than the length of the response, then offset is replicated to the appropriate length.
ARG(R@@)
a square matrix of dimension maximum cluster size containing the user specified correlation.  This is only appropriate if corstr="fixed".
ARG(b@@)
an initial estimate for the parameters.
ARG(tol@@)
the tolerance used in the fitting algorithm.
ARG(maxiter@@)
the maximum number of iterations.
ARG(family@@)
a family object: a list of functions and expressions for
defining link and variance functions.  Families supported
in gee are gaussian, binomial, poisson, Gamma, and quasi;
see the glm and family documentation.  Some links are
not currently available: cloglog, 1/mu^2 and sqrt have
not been hard-coded in the cgee engine at present.
The inverse guassian variance function is not available.
All combinations of remaining functions can be obtained
either by family selection or by the use of quasi.
Future releases will allow S-coded functions
for link, variance and correlation function specification.
ARG(corstr@@)
a character specifying the Correlation structure. The following are permitted:
     "independence"
     "fixed"
     "stat_M_dep"
     "non_stat_M_dep"
     "exchangeable"
     "AR-M"
     "unstructured"
ARG(Mv@@)
When the corstr is "stat_M_dep", "non_stat_M_dep", or "AR-M" then Mv must be specified.
ARG(silent@@)
a logical variable controlling whether parameter estimates at each iteration are printed.
ARG(contrasts@@)
a list giving contrasts for some or all of the factors appearing
in the model formula.  The elements of the list should have the
same name as the variable and should be either a contrast matrix
(specifically, any full-rank matrix with as many rows as there are
levels in the factor), or else a function to compute such a matrix
given the number of levels.
ARG(scale.fix@@)
a logical variable; if true, the scale parameter is fixed at
the value of scale.value
ARG(scale.value@@)
numeric variable giving the value to which the scale parameter
should be fixed; used only if scale.fix == T.
)
PARA BOLD(VALUE)
An object of class "gee" representing the fit.
PARA BOLD(SIDE EFFECTS)
PARA BOLD(REFERENCES)
Liang, K.Y. and Zeger, S.L. (1986).
Longitudinal data analysis using generalized linear models.
Biometrika 73: 13-22. 
Zeger, S.L. and Liang, K.Y. (1986).
Longitudinal data analysis for discrete and continuous outcomes.
Biometrics 42: 121-130.
PARA BOLD(SEE ALSO)
LANG(glm), LANG(lm), LANG(formula).
BLANK

PARA BOLD( )
BLANK

EXAMPLES(
gee(y~x,id)   # Gaussian model with independent correlation structure
gee(d ~ x1+x2,id,link="logit",family=binomial,corstr="exchangeable")
gee(d~x1+x2,id=id,family=binomial(link=(probit)),corstr="stat_M_dep",M=1)
)
PARA ITALIC(File automatically converted from S(-PLUS) help format)
