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)