Displaying 4 results from an estimated 4 matches for "gee1".
Did you mean:
gee
2010 Feb 10
1
using step() with package geepack
...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:
library(geepack)
data(dietox) # the example from ?geeglm
# GEE model, using gaussian instead of poisson:
gee1 <- geeglm(Weight~1, data=dietox, id=Pig, family=gaussian, corstr="ar1")
gee2 <- update(gee1, ".~.+Time")
anova(gee1, gee2) # gives a p-value for adding the Time variable
step(gee2, direction="forward") # fails
cheers,
Steve
[[alternative HTML vers...
2004 Mar 19
0
yags, GEEs, and GLMMs
...on mixed models with binary response variables?......).
I have used the glmmPQL function (about which I also have a whole list of questions....) as mentioned in MASS, but now I was told that GEE models might be better suited for my data. So I got the yags package and tried a model like this:
> gee1<-yags(alive~sex*weight, family=binomial,id=id,corstr='exchangeable', alphainit=0.)
But then the summary only gives me this:
> summary(gee1)
Length Class Mode
0 yagsResult list
How do I get R to give me the parameter estimates, std errors, t- and p-va...
2004 Mar 19
0
yags, GEEs and GLMMs
...on mixed models with binary response variables?......).
I have used the glmmPQL function (about which I also have a whole list of questions....) as mentioned in MASS, but now I was told that GEE models might be better suited for my data. So I got the yags package and tried a model like this:
> gee1<-yags(alive~sex*weight, family=binomial,id=id,corstr='exchangeable', alphainit=0.)
But then the summary only gives me this:
> summary(gee1)
Length Class Mode
0 yagsResult list
How do I get R to give me the parameter estimates, std errors, t- and p-v...
2008 Mar 05
1
problem with geepack
...n
Individual = a code indicating the id of the specimen.
I have up to 12 replicate observations for individual. Some observation
are missing and I organized the data frame to have exactly 12 rows for
each specimen, with NAs where there is a missing observation.
The following model worked fine:
gee1<-geese(Activity~Water, id=Individual, data=dataF, family=gaussian)
but when I use other variables e.g
gee2<-geese(Activity~Hgt+Size+Water, id=Individual, data=dataF,
family=gaussian)
I get the error message
"Error in geese.fit(x, y, id, offset, soffset, w, waves, zsca, zcor,
corp,...