similar to: using step() with package geepack

Hi all I am analyzing a data set containing information about the behaviour of marine molluscs on a vertical wall. Since I have replicate observations on the same individuals I was thinking to use the geepack library. The data are organised in a dataframe with the following variables Date = date of sampling, Size = dimensions (mm) Activity duration of activity (min) Water = duration of

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

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

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

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>

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

Hi, Let's say one has an object with multiple classes, and a generic function to apply to it has associated S3 methods for more than one of those classes. Further, the method it chooses (I do not know how; some order in the class vector?) is not the suitable one and it produces an error. Would there be some way to make the generic function to choice the correct method, or in case that for

Dear R users I am analysing data from a group of twins and their siblings. The measures that we are interested in are all correlated within families, with the correlations being stronger between twins than between non-twin siblings. The measures are all calculated from survey answers and by definition have limited ranges (e.g. -5 to +5), though within the range they are approximately normally

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

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 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:

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

Hi there, I just installed the package under windows, and it works very well. However, as I tried to install it on a cluster which uses linux OS, "Linux lionxc.aset.psu.edu 2.6.18-194.el5 #1 SMP Tue Mar 16 21:52:39 EDT 2010 x86_64 x86_64 x86_64 GNU/Linux" I always failed. The message during the installation process is as follows: > install.packages("geepack") Warning

*H*i, I am trying to fit a GEE model with *geeglm* function. The model is the following: Modelo<-geeglm(sqrt ~Tra+ Mes, id=Lugar , data=datos, family=gaussian(identity), corstr="independence") *Tra( is a experimental treatment, 2 levels)*, *Mes* (is the month of take data, 4 levels) and *Lugar* (is the site of study, 3 levels) are categorical variables and *sqrt* (sqrt of Total

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

I am trying to create a loop that will perform a series of analyses. I am using geeglm from geepack, which fails if there are any null values. Creating a subset solves this, but do not seem to be able to set the subset dynamically based on a changing variable. while (j <= y.num) { strSubset = as.character(df.IV$IV[j]) #Gives column name in quotes df.data.sub = subset(df.data,

I'm trying to fit a generalized mixed effects model to a data set where each subject has paired categorical responses y (so I'm trying to use a binomial logit link). There are about 183 observations and one explanatory factor x. I'm trying to fit something like: (lmer(y~x+(1|subject))) I also tried fitting the same type of model using glmmPQL from MASS. In both cases, I get a

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 (

? 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 included 2 internal functions in the package 'dlnm', called 'mkbasis' and 'mklagbasis'. Despite they are not meant to be called by the users, I included them in the namespace in order to make them available, keeping the process more transparent and giving the opportunity to change or improve them. I included an help page 'dlnm-internal.Rd' to