我使用各种函数拟合我的离散计数数据进行比较。我使用 geepack 拟合了一个GEE模型，log(count)使用lme（ nlme ）的线性混合效果模型，GLMM使用{{1 （ lme4 ），以及在R中使用glmer（ gamm4 ）的GAMM。
我有兴趣比较这些模型，并希望绘制一组新数据（预测变量）的预期（预测）值。我的目标是比较特定条件下每个模型的预测效果（x变量）。特别感兴趣的是边际（GEE）和条件估计之间的比较。
我认为我的主要问题可能是使用正确的标签和属性等以正确的形式获取新数据。我仍然是一个R新手，并且很难与这些东西斗争（不幸的是，在我的大学没有这个课程。）
我目前已经安装了模型
gamm4并且可以毫无问题地提取其固定效果系数和标准误差。我也没有问题从对数比例转换它们或估计考虑随机效应的置信区间。
我还有我的新数据框gee1 lme1 lmer1 gamm1
，它有23个变量的365个观测值（一年中每一天的平均环境数据）。
我坚持如何预测新的计数估计值。我玩了model.matrix函数，但无法使它工作。例如，我试过：
newdat任何建议或好的参考将不胜感激。任何人都可以帮助解决上面的错误吗？