Displaying 1 result from an estimated 1 matches for "rbinarydata".
Did you mean:
binarydata
2005 Jan 05
0
lme, glmmPQL, multiple random effects
.... (data available as an *.rda for anyone who wants it...)
Following the advice from various postings from R-help and from Pinhero
and Bates, I can successfully (I think) code an lme() model, using
pdBlocked, pdIdent and a dummy grouping variable in a grouped data
object.
> load("~/RbinaryData.rda")
> dd<-data.frame(RbinaryData,dum=factor(rep(1,120)))
> names(dd)
[1] "Cross" "Female" "Male" "Mate1" "dum"
> summary(dd)
Cross Female Male Mate1 dum
RR:30 1 : 6 1 : 6 Min....