search for: stat_m_dep

Dear all: Y=a+bX1+cX2 In the Gee (Generalized Estimation Equation) arguments: The arument Corstr has sveral choices: "independence" "fixed" "stat_M_dep" "non_stat_M_dep" "exchangeable" "AR-M" "unstructured" What does each term mean? How do I choose among them? How do I know the correlation structure of my repeated measures data? Can non-normally distributed X1, X2 be used? Thank you very m...

...and (package gee) doesn't work when the id groups are large, as in my dataset (observations from several summer seasons, for which I imagine an AR correlation structure within each season). The command geeglm (package geepack) seems to work, but provides only few corstr choices (for example not stat_M_dep, which can be useful to investigate models with different correlation structures). Any suggestions? Thanks so much for your time

...a marginal model via GEE but I'm getting strange results and few problems. If I set the working correlation as exchangeable I'm getting the same fitting when I set as independent. Comparing to SAS results it shouldn't happen. If I try to use another working correlation (like AR-M or stat_M_dep), R just exits without giving any error message. Another doubt is how R estimates the scale parameter for Poisson? I get 1.36 in R, 1.25 in SAS (estimation by the square root of DEVIANCE/DOF) and 1.17 (estimation by the square root of Pearson's Chi-Square/DOF). I'm using R 1.4.1 but all...

...log") fams <- c("gaussian", "poisson", "binomial", "Gamma", "quasi") varfuns <- c("constant", "mu", "mu(1-mu)", "mu^2") corstrs <- c("independence", "fixed", "stat_M_dep", "non_stat_M_dep", "exchangeable", "AR-M", "unstructured") linkv <- as.integer(match(c(family$link), links, -1)) famv <- match(family$family, fams, -1) if (famv < 1) stop("unknown family") if (famv <...