0

我正在使用 stats() 3.5.2 运行 manova：

  • 参与者 1:20
  • 主体因素之间的性别
  • 分组为主体因素
  • 焦虑作为依赖量度
  • BAC 作为从属测量

数据集如下：

treat4 = data.frame (
   participant = rep(1:20,3),
   gender = factor (rep(c(rep("male", 10), rep ("female", 10)),3)),
   group = factor (c(rep("control",20), rep("run",20), rep("party",20))),
   anxiety = round(c(rnorm(20, mean=55, sd=5),rnorm(20, mean=20, sd=5),rnorm(20, mean=75, sd=5))),
   BAC = round(c(rep(0.01,20), rep(0.01,20), rnorm(20, mean= 0.09, sd=0.01)),2))

我应用 manova() 函数并总结如下：

mod = manova(cbind(anxiety,BAC) ~ gender + Error(group),data=treat4)
summary (mod)

这就是我得到的：

Error: group
          Df Pillai approx F num Df den Df Pr(>F)
Residuals  2                                     

Error: Within
          Df   Pillai approx F num Df den Df Pr(>F)
gender     1 0.013447  0.37482      2     55 0.6892
Residuals 56

有几个问题：

1) 性别似乎被视为受试者内因素

2）我没有得到任何关于组因素的统计数据

有什么帮助吗？

4

3 回答 3

0

如果焦虑和 BAC 是您的因变量，您可以使用 cbind 将它们放在 tilda (~) 的左侧，以表示多变量响应，并用于Error()指定组内效应（或随机效应）。tilda (~) 右侧的其余部分是您的组间效果（或固定效果）：

manova(cbind(anxiety,BAC) ~ gender + Error(group),data=treat4)

Call:
manova(cbind(anxiety, BAC) ~ gender + Error(group), data = treat4)

Grand Means:
    anxiety         BAC 
49.96666667  0.03766667 

Stratum 1: group

Terms:
                 Residuals
anxiety           33156.63
BAC             0.09185333
Deg. of Freedom          2

Residual standard errors: 128.7568 0.2143051

Stratum 2: Within

Terms:
                   gender Residuals
anxiety            8.0667 1527.2333
BAC                0.0000    0.0034
Deg. of Freedom         1        56

Residual standard errors: 5.222262 0.007807201
Estimated effects are balanced
于 2020-03-04T22:44:23.247 回答
0

我知道这有点晚了，但我遇到了同样的问题，我认为您可以通过以下方式简单地解决它：

summary(manova(cbind(anxiety,BAC) ~ gender + group + Error(factor(participant)),data=treat4))

基本上，您需要将组添加为 IV（通过执行+ group）。然后您使用Error()来指示它需要如何识别唯一主题，它需要通过参与者编号而不是组来执行此操作。

不要忘记让参与者成为一个因素，否则会出问题！

于 2021-11-29T10:12:56.833 回答
0

感谢@StupidWolf 的回答。

但是，当我将 summary () 应用于模型时：

summary(manova(cbind(anxiety,BAC) ~ gender + Error(group),data=treat4))

我得到以下信息：

Error: group
          Df Pillai approx F num Df den Df Pr(>F)
Residuals  2                                     

Error: Within
          Df   Pillai approx F num Df den Df Pr(>F)
gender     1 0.039097   1.1189      2     55  0.334
Residuals 56

有几个问题：

1) 性别似乎被视为受试者内因素

2）我没有得到任何关于组因素的统计数据

于 2020-03-05T16:37:06.350 回答