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.? 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 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