similar to: p and wald values intra-groups geeglm: geepack

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 (

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091126/7134fc17/attachment-0001.pl>

Hi, I have a follow analysis Trat1 = quantitative variable Trat2 = qualitative variable with 3 levels (A, B, C) Trat3 = qualitative variable with 3 levels (D, E, F) Resp = Response I try to get the parameters to compare with zero, so I make this model: glm(Resp~Trat1*Trat2+Trat1*Trat3-Trat1-1) The -Trat1 is to make comparison of slope with zero. The -1 is to make comparison os intercept with

Dear all I am struggling with how to deal with missing values using geeglm. I know that geeglm only works with complete datasets, but I cannot seem to get the na.omit function to work. For example assuming DataMiss contains 3 columns, each of which has missing observations, and an id column with no missing info then identifies the clusters. Outcome: 2 level integer Predictor: numeric variable

Hello R-helpers, I would like to calculate the standard error for the predicted value from geeglm. As an example, I would like to calculate the GEE mean of treatments and their standard error. I first specified the model as mod <- geeglm(resp ~ trt, data=dat,id=id,family=Gaussian,corstr="ar1",weights=weight) Then I predicted the GEE mean and se using the following code

event.nab.2 is 0/1 and I dichotomized va to get va.2 to see if I could get geeglm to work. glm has no problem with the data but geeglm chokes. Each subject (patient.id) has at most 2 observations and more than 3/4 of the subjects have 2 observations. I have even worse problems trying to use glmmPQL from MASS and worse still trying to use lmer from lme4. But I figured a marginal model would work.

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,

Dear R users, could you please help me figure out why I am getting an error? Initially my data looks like this: > attributes(compl)$names [1] "UserID" "compl_bin" "Sex.x" "PHQ_base" "PHQ_Surv1" "PHQ_Surv2" "PHQ_Surv3" [8] "PHQ_Surv4" "EFE"

Dear all, I am trying to use the GEE methodology to fit a trend for the number of butterflies observed at several sites. In total, there are 66 sites, and 19 years for which observations might be available. However, only 326 observations are available (instead of 1254). For the time being, I ignore the large number of missing values, and the fact that GEE is only valid under MCAR. When I run the

Dear all, We have got data (response and predictor variables) for each country of the world; I started by fitting standard GLM and tested for spatial correlation using variogram models (geoR) fitted to the residuals of the GLM. Spatial autocorrelation is significant. Therefore, I think about using general estimation equations (geeglm or geese in geepack) allowing for residual spatial

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 ~

Hi, I'm using geeglm function to account for the repeated measure. fit1<- geeglm( binary.outcome ~ age + race + gender + fever.yes.no, data=mydata, id=ID, family=binomial, corstr="exchangeable") summary(fit1)$coef gives too large estimates and standard deviation: Estimate Std.err Wald Pr(>|W|) (Intercept) 3.07e+16

Hi, I managed to make R core dump (linux and Mac OSX versions), but I think I've figured out why. First, here's the message I get on core dump (on linux - no message on Mac): R: ../inst/include/tnt/fmat.h:529: TNT::Vector<T> TNT::matmult(const TNT::Fortran_Matrix<T>&, const TNT::Vector<T>&) [with T = double]: Assertion `A.num_cols() == x.dim()'

Hi all, I am trying to do a generalized estimating equation (GEE) with the "geepack" package and I am not 100% sure what exactly the "id" argument means. It seems to be an important argument because results differ considerably defining different clusters. I have a data set of counts (poisson distribution): numbers of butterfly species counted every month during a period of

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

I'm using the package geepack to fit GEE models. Does anyone know of methods for add1 and drop1 for a 'geeglm' model object, or perhaps a method for extractAIC based on the QIC of Pan 2001? I see there has been some mention of this on R-help a few years ago (RSiteSearch("QIC")). The package does provide an anova method for its model objects, and update() seems to work:

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070905/6d1002c1/attachment.pl

Dear R-users, I was doing some Poisson regression using the geepack package and I found a "strange" behaviour of the geeglm command when defining an offset. Maybe it's my limited knowledge of R syntax, or maybe it's something else. Let's make an example. After loading the geepack library, you may write

Hello R Forum members. I have installed for my statistician user, apparently without error, both the concord and geepack packages. The target system is R 2.10.1 on a 64-bit RedHat Enterprise Linux platform. However when she attempts to invoke a function in geepack, for example... geeglm((abuse_total ~ case),id=mother, family="poisson") the resultant error is Error: could not find

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