r - 将 MCMCglmm 对象放入 Anova() 或 Manova()

Question

我正在使用MCMCglmm包来执行混合模型分析。如果我的数据如下：

df_test <- data.frame(ID = 1:10, a = 171:180, b = 71:80 + rnorm(10),
  age = factor(c(rep("young", 3), rep("mid", 4), rep("old", 3)),
    levels = c("young", "mid", "old")))

对于线性模型，我可以轻松summary(Manova(lm(cbind(a, b) ~ age + 0, data = df_test)))获取表格并查看对年龄的明显影响（请忽略此处的共线性问题）。

但是，如果我用MCMCglmmID来考虑混合效应呢？假设该MCMCglmm对象称为“mc_model”，我如何使用Manova或类似的方法来查看（多个）分类（具有多个级别）变量的影响？

谢谢！