时间戳

首页 > 解决方案 > 循环数据的重复混合模型或循环数据的混合模型

r - 循环数据的重复混合模型或循环数据的混合模型

问题描述

我想对循环数据进行重复测量 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

解决方案

由于对我来说没有意义的原因（因为我没有读过论文），显然你的“嵌套变量”（Subject）需要特别是数字- 并且作为整数类型不符合条件！

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

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

此外：

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

推荐阅读