Error in code

[guha22]

Dear Dr. Bolker,

Hope this email finds you well.

I have been trying to run the code on this page. The phylo_glmm code runs fine. However, the phylo_lmm code gives me the following error (though it still seems to fit a model).

In case you could let me know what I have done wrong I would be grateful.

Regards,
Guha

---

CODE:::

phylo_lmm_fit <- phylo_lmm(phen~cofactor+(1|phylo) +(1|obs),

•                        data=datG,phylo=phylo,
  

•                        control=lmerControl(check.nobs.vs.nlev="ignore",check.nobs.vs.nRE="ignore"),
  

•                        phyloZ=phyloZ)
  

Warning messages:
1: In split.default(x, g) :
data length is not a multiple of split variable
2: In split.default(seq_along(x), f, drop = drop, ...) :
data length is not a multiple of split variable

lme4_resLMM <- rbind(coef(summary(phylo_lmm_fit))[,1:2],

•                  cbind(matrix(unlist(VarCorr(phylo_lmm_fit))),NA))
  

phylo_lmm_fit
Linear mixed model fit by REML ['lmerMod']
REML criterion at convergence: 1563.8
Random effects:
Groups Name Std.Dev.
phylo (Intercept) 41.406
obs (Intercept) 0.443
Residual 9.690
Number of obs: 200, groups: phylo, 398; obs, 200
Fixed Effects:
(Intercept) cofactor
37.348 5.207