r - manova () in R with between and within subject 因子

Question

私は 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) グループ因子の統計が得られない

何か助けはありますか？