r aov（）有2个相关度量

Question

[我正在尝试使用一项独立措施（组，三个级别：控制，跑步，聚会）和两项相关措施（焦虑和BAC）进行重复测量方差分析。

这里是data.frame，

我有以下数据集，我想在其中执行manova（）：

• 1个主题内因素（group）
• 1主体间因素（性别）
• 2个相关指标：焦虑和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（）模型中包括主体间因素gender的帮助。

谢谢！

我知道ezanova（）仅占1个从属度量。我试图使用aov（）函数来说明2个从属度量，但不清楚在Error参数中输入什么。

感谢您的帮助

Answer 1

如果焦虑和BAC是您的因变量，则将它们与cbind一起放在tilda（〜）的左侧，以指示多变量响应，然后使用Error()指定组内效应（或随机效应）。蒂尔达（〜）右侧的其余部分是组间效果（或固定效果）：

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

Answer 2

感谢@StupidWolf的回答。

但是，当我将摘要（）应用于模型时：

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）我没有任何关于组因子的统计信息>