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 - [EMAIL PROTECTED]
[[alternative HTML version deleted]]
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
