10  brms

library(brms)

Amat <- as.matrix(nadiv::makeA(gryphonped))
gryphonRM$animal_pe <- gryphonRM$animal



model_simple1.1 <- brm(
  laydate ~ 1 + (1 | gr(animal, cov = Amat)) + (1 | animal_pe),
  data = gryphonRM,
  family = gaussian(),
  data2 = list(Amat = Amat),
  chains = 2, cores = 2, iter = 1000
)

summary(model_simple1.1)
plot(model_simple1.1)

Wilson, A. J., D. Réale, M. N. Clements, M. M. Morrissey, E. Postma, C. A. Walling, L. E. B. Kruuk, and D. H. Nussey. 2010. An ecologist’s guide to the animal model. Journal of Animal Ecology 79:13–26.