Displaying 1 result from an estimated 1 matches for "mate1".
Did you mean:
  mate
  
2005 Jan 05
0
lme, glmmPQL, multiple random effects
...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.   :0.0000   1:120
  RW:30   2      : 6   2      : 6   1st Qu.:0.0000
  WR:30   3      : 6   3      : 6   Median :1.0000
  WW:30   4      : 6   4      : 6   Mean...