Question

我是ez包新手，我一直试图使用ezMixed运行混合效果anova，因为我的重复测量数据是不平衡的，并且对于一些重复观察存在一些缺失的数据。我的响应变量'count'是计数数据的平均值，不遵循泊松分布。首先，我不确定我是否使用ezMixed正确指定了我的模型。使用lme（），我的模型看起来像这样：

lmefitmix2 = lme(count ~ MPA*timepointID,random=~1|siteID/timepointID ,na.action=na.omit,data = notolabrusMav2,method="ML")

我的重复测量（在因子内）是嵌套在随机因子siteID中的timepointID。我感兴趣的是测试分组变量“MPA”（效果之间）的主要影响以及重复测量“时间点ID”的影响以及它们的相互作用。我明白如果我的设计已经平衡，我可以使用ezAnova并指定wid =。siteID和= .timepoint ID，但我不知道如何在ezMixed中翻译它。我的尝试如下，但它使得看起来像siteID和timepointID：

er = ezMixed(
    data = notolabrusmav2
    , dv = .(count)
    , random = .(siteID,timepointID)
    , fixed = .(MPA,timepointID)
    , family = 'gaussian'
)
print(er$summary)

我该如何解决这个问题？

另外，print（er $ summary）命令给出了以下结果：

er$summary

           effect errors warnings       bits

1             MPA  FALSE    FALSE   1.710731

2     timepointID  FALSE    FALSE  26.514379

3 MPA:timepointID  FALSE    FALSE -24.895950

所以我理解这里最有力量证据的主要影响是'timepointID'，而两者的相互作用显示出很少的证据强度。

我试图使用ezPredict和ezPlot2生成和绘制预测，如下所示：

myFit = er$models$timepointID$unrestricted
myPreds = ezPredict(myFit)
myPlot = ezPlot2(preds=myPreds,x=timepointID)
print(myPlot$plot)

这给了我以下错误消息：

Error in order(boots$iteration) : argument 1 is not a vector
 print(myPlot$plot)
Error in print(myPlot$plot) : 
  error in evaluating the argument 'x' in selecting a method for function 'print': Error: object 'myPlot' not found

我一直在谷歌搜索它，但我似乎无法找到帮助我修复它的提示......

谢谢你提前帮助; 0

使用ezMixed指定混合模型，并使用ezPlot2（ez包）获取预测和绘图

0 个答案: