[R-sig-ME] Complete separation example

[Jarrod Hadfield]

Hi Vince,

Sex effect on horn type in Soay sheep is completely separated (the  
trait is partly sex-limited):

data(SShorns)
m1<-glm(I(horn=="polled")~sex, data=SShorns, family=binomial)

Cheers,

Jarrod


Quoting Vincent Dorie <vjd4 at nyu.edu> on Tue, 4 Mar 2014 12:20:41 -0500:

> Hi list, can anyone point me to a data set that exhibits complete or  
> quasi-complete separation when fit as a glmm? It's for a vignette,  
> so something public would be ideal.
>
> With thanks,
> Vince
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>



-- 
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.