我想从 WRS2 包 ( https://rdrr.io/cran/WRS2/man/ancova.html ) 运行一个强大的 ANCOVA。“比率”是我的 DV，“诱饵”是预测变量，“会话”是协变量。
我使用以下功能：
ancova_model = ancova(ratio ~ bait + session, data = training_data)
这是我的数据集：
但是，我收到以下错误消息：
if (nuhat < 2) stop("自由度必须大于或等于 2") 中的错误：需要 TRUE/FALSE 的缺失值另外：警告消息：1：In min(sub[vecn >= 12 ]) : min 没有非缺失参数；返回 Inf 2: In max(sub[vecn >= 12]) : max 没有非缺失参数；返回-Inf
任何想法我做错了什么？
问题是稳健的 ANCOVA 确定协变量的 5 个点，然后在每个点上对独立组执行 yuen t 检验。因此，我怀疑在这么小的样本中，协变量对两组都具有完全相同的值，最终可能会在每个点对每个组进行一次观察（因此 df < 2）。对于每组只有一次观察，无法计算 t 检验，因此无法计算误差。
解决此问题的一种方法是，您可以尝试更改测量组比较的协变量点的数量和类型，如下所示ancova(ratio ~ bait + session, data = training_data, pts = c(5,10))：这可能会给你一个结果，但也可能会抛出警告消息，表明测试已经进行了太少的观察，这意味着你不能真正相信结果。
在我看来，协变量无论如何都没有多大意义。如果我没记错的话，这只是表明两组（水果和兔子）是在同一会话中记录的。因此，这些群体实际上可能是相互依赖的。您可以使用函数来代替依赖 yuen 的 t 检验yuend。