Dear R-ers,
I am just a simple 'end-user' of R and am trying to analyse data with a
binary response variable (dead or alive) in relation to weight and sex (of young
birds). As some of the birds have the same biological mother, I am using mixed
models with the identity of the mother as a random factor. (please, Mick
Crawley, when are you going to write a chapter 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-values
etc of the model?
And if someone is willing to answer questions about the glmmPQL function, please
let me know and I'll send you the most urgent questions, but I don't
want to bother everyone with it.
Many thanks,
Ellen
P.S. Sorry, i think I sent this same message to the wrong address earlier,
don't know where it went to.
Ellen Kalmbach - Animal Ecology Group - University of Groningen - Kerklaan 30 -
NL-9751NN Haren - Netherlands - +31-503633408 - e.kalmbach@biol.rug.nl
[[alternative HTML version deleted]]
