Манова () в R с между и внутри субъектных факторов - PullRequest
Новая
       66

Манова () в R с между и внутри субъектных факторов

0 голосов
/ 04 марта 2020

Я использую stats () 3.5.2 для запуска мановы с:

  • участником 1: 20
  • полом между субъектным фактором
  • в группе как в рамках предметного фактора
  • тревога как зависимая мера
  • BA C как зависимая мера

Набор данных следующий: 

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) Я не получаю статистику по групповому фактору

Любая помощь?

Ответы [ 2 ]

0 голосов
/ 05 марта 2020

Спасибо @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) Я не получаю статистику по групповому фактору

0 голосов
/ 05 марта 2020

, если тревога и BA C являются вашими зависимыми переменными, вы помещаете их в левую часть тильды (~) с помощью cbind, чтобы указать многомерный ответ, и используйте 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
...