循环数据的重复混合模型或循环数据的混合模型答案

【问题标题】：Repeated Mixed model for circular data or mixed model of circular data循环数据的重复混合模型或循环数据的混合模型
【发布时间】：2020-09-26 21:11:58
【问题描述】：

我想对循环数据进行重复测量 ANOVA，但除了 Watson-Williams 检验（ANOVA 等效项）外，我还没有找到一个很好的例子，但我认为它不能解释重复测量。然后我尝试使用bpnreg 包 (https://www.frontiersin.org/articles/10.3389/fpsyg.2018.02040/full) 拟合混合模型，该包使用认知数据对混合模型进行建模。我不知道为什么我不能运行它；它一直在我的数据上出错。我的数据是 2000-2017 年不同年份 10 个不同城市的三个变量的平均值。我想看看一个变量的平均值是否依赖于其他变量，并且变量之间以及城市和年份之间是否存在相互作用。该模型不采用城市名称等分类数据，因此我将城市编码为数字（不确定这是否正确），如果有人知道，请告诉我。以下是我的数据和代码。任何帮助将不胜感激。

以下是我针对 10 个不同城市的数据格式。

Cities   Year   MeanP     MeanE      MeanG      Subject
CityA    2000   158.1235  158.1415   158.1512   1
CityA    2005   158.1735  158.1215   158.1912   1
CityA    2010   158.1035  158.1615   158.1112   1
CityA    2015   158.1935  158.1215   158.1012   1
CityA    2017   158.1335  158.1715   158.1312   1
CityB    2000   154.2347  154.2567   154.2744   2

我的代码如下

library(bpnreg)
data<- read.table("meanscities.csv", header = T, sep=",")
modelfit <- bpnme(MeanG ~ MeanP + MeanE + (1|Subject),
                  data = data,
                  its = 10000, burn = 1000, n.lag = 3, seed = 101)

以下是错误：

mmme(pred.I, data, pred.II) 中的错误：并非所有嵌套变量都是类数字。

但数据是列表格式，除城市名称外都是数字。

【问题讨论】：

标签： r

【解决方案1】：

出于对我来说没有意义的原因（因为我还没有阅读论文），显然您的“嵌套变量”（Subject）需要专门为 numeric - 并且是整数类型不符合条件！

data$Subject <- as.numeric(dd$Subject)

让我克服了这个错误。不幸的是，我遇到了另一个更难弄清楚的错误，但希望这只是因为我使用了您在上面发布的微不足道的数据集...

另外：

您的响应变量必须以弧度为单位，否则您的答案可能毫无意义...
您确定bpnme 不允许使用分类预测变量吗？将分类预测器编码为整数几乎没有意义。如果bpnme 不允许，而您确实需要它们，您可以创建自己的虚拟变量...

【讨论】：