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


