similar to: empty cell when running GEE for binary data

Hi A quick question. Standard errors reported by gee/yags differs from the ones in geeglm (geepack). require(gee) require(geepack) require(yags) mm <- gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable") mm2 <- geeglm(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable", std.err = "san.se") mm3 <- yags(breaks ~

Hello R-users: I am using "yags" for fitting GEE which is giving me the same result as "Proc GENMOD". Now I have couple of questions related to yags output. (By the way, someone told me to run the geeglm for the same analysis and I did run but did not get the same result as of genmod and don't know how to correct the geeglm codes so that all three will be same!)

Dear R users Recently I received advice from this fine group on gee() and sample weights One suggestion was to use geeglm() I hope someone can help me to solve a problem that arises when converting a code from gee to geeglm. *Here is a code that I wrote with the original data, not weighted: * > m1 <- gee( Bin ~ educ+agemean+ residencysize + yearx , id = rad09 , data = Males,

Hi, I'm using 'geepack' to run Generalized Estimating Equations. I'm aware that I can use anova to compare two models, but would it be possible to test QIC on R? It seems that there were similar questions a couple of years ago, but the question has not been answered yet. I'd appreciate if someone could show me the code! Thank you, Sachi [[alternative HTML version

I have included data at the bottom of this email. It can be read in by highlighting the data and then using this command: dat <- read.table("clipboard", header = TRUE,sep="\t") I can obtain solutions with both of these: library(gee) fit.gee<-gee(score ~ chem + time, id=id, family=gaussian,corstr="exchangeable",data=dat) and library(yags) fit.yags <-

Dear R users, I fitted a GEE model using the function 'geese' (or 'geeglm') with user defined correlation matrix. I want to get the var-cov matrix of the regression coefficients. But the output provides only limited information. I would be very much thankful if you could kindly let me know how to get it..since I am struggling lot getting this. Thanks -- View this message in

Hi, having used both the gee pacakge and the geepack package, i am unsure of how to interpret the results. Here are the results from the geeglm function from the geepack package > gee2<-geeglm(data$erythema~data$product, data = data, id=subject, > family=binomial, corstr="independence") Warning message: In model.response(mf, "numeric") : using

Dear all, I need to fit a gee model with an auto-regressive correlation structure and I faced some problems. I attach a simple example: ####################################################### library(gee) library(geepack) # I SIMULATE DATA FROM POISSON DISTRIBUTION, 10 OBS FOR EACH OF 50 GROUPS set.seed(1) y <- rpois(500,50) x <- rnorm(500) id <- rep(1:50,each=10) # EXAMPLES FOR

Hi, I'm running GEE using geepack. I set corstr = "ar1" as below: > m.ar <- geeglm(L ~ O + A, + data = firstgrouptxt, id = id, + family = binomial, corstr = "ar1") > summary(m.ar) Call: geeglm(formula = L ~ O + A, family = binomial, data = firstgrouptxt, id = id, corstr = "ar1") Coefficients:

Hello, I am interested in running generalized estimating equation models in R. Currently there are two main packages for doing so in R, geepack and gee. I understand that even though one can obtain similar to almost identical results using either of the two, that there are differences between the packages. The paper that introduces the geepack package (

I'm writing a paper aimed at motivating the use of GEE within the field of economics. However, after computing using the geeglm function, I noticed there's one intercept in the summary output. I assume this means the function is pooling the data. That means my code is not what I want. I want a "fixed effects" model, meaning I want the intercept to vary by cluster. Here's

Dear List: I found an error when I called the 'gee' function. I cannot solve and explain it. There are no errors when I used the 'geeglm' function. Both functions fit the gee model. The project supervisor recommends me to use the 'gee' function. But I cannot explain to him why this error happens. Would you help me solve this problem? I appreciate your help. In

Dear all, Is it possible to incorporate a nested design in GEE? I have measurements on trees that where measured in two years. The trees are nested in plots. Each plot contains 24 trees. The number of plots is 72. Hence we would expect 2 * 24 * 72 = 3456 data points. A few are missing, so we end up wih 3431 data points. This is what I have tried until now. #assuming independence between trees

Hi, I'm analyzing my data using GEE, which looks like below: > interact <- geeglm(L ~ O + A + O:A, + data = data1, id = id, + family = binomial, corstr = "ar1") > summary(interact) Call: geeglm(formula = lateral ~ ontask + attachment + ontask:attachment, family = binomial, data = firstgroupnowalking, id = id, corstr = "ar1") Coefficients:

Goal: use GEE or GLMM to analyze repeated measures data in R GEE problem: can?t find a way to do GEE with negative binomial family in R GLMM problem: not sure if I?m specifying random effect correctly Study question: Does the interaction of director and recipient group affect rates of a behavior? Data: Animals (n = 38) in one of 3 groups (life stages): B or C. Some individuals (~5)

Hi , I use library(gee),library(geepack),library(yags) perform GEE data analysis , but all of them cannot do variable selection! Both step and stepAIC can do variable selection based on AIC criterion under linear regression and glm, but they cannot work when model is based on GEE. I want to ask whether any variable selection function or package under GEE model avaliable now? Thanks! Best,

Hello, Is it possible to derive Pan's QIC (2001 Biometrics 57:120) from either a fitted gee() object in the gee package or from a geese() fit in the geepack package? If so, would anyone be kind enough to provide me with code to do so? I realize that QIC is part of the output from yags() but I would like to use one of the other functions. Thanks. Richard

Dear all, I am trying to do a GEE on count data and I am having problems with how to order the data. Below is a simplified example of what my data looks like.. Route Time Day Pass Distance 1 30 1 4 0 1 60 1 12 200 1 120 1 25 600 1 30 2 8 0 1 60 2 17 200 1

Hello, R exits every time I want to fit a GEE with AR-1 correlation structure. It just closes itself without any error. My Data has about 3300 observation (using the big data with about 20 000 observation it doesn't work at all). I'm a Windows 2000 User of R 2.0.1 but it was the same problem on a MacOSX. Is it my fault? Thank you for your help! Johanna Brandt

i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/win-library/2.6' trying URL