Внутри-предметный контраст в R от смешанного ANOVA - PullRequest
Новая
       15

Внутри-предметный контраст в R от смешанного ANOVA

0 голосов
/ 22 апреля 2019

Я провел смешанный анализ ANOVA 3х4 и должен следить за контрастом предметов, однако я не знаю, как это сделать.

Я повторяю анализ, сделанный в статье на R, чтобы сделать его воспроизводимым, и авторы оригинала провели ANOVA, а затем внутриконтрастные контрасты. Они сообщили F (df1, df2) = .... p = ..., что заставляет меня предположить, что они сделали какой-то дополнительный анализ ANOVA, но я не могу понять, что и с какими переменными.

Пример сбалансированных данных с примерами переменных (исходные данные не сбалансированы) 


library(tidyverse)
library(ez)

dat <- tibble(subj = rep(c(1:24), each = 4),
              trial = rep(c(1:4), 24),
              score = rnorm(96, 0.8, 0.2),
              cond = rep(c(1:4), each = 24))

dat$trial <- factor(dat$trial, levels=
                                c("1",
                                  "2",
                                  "3",
                                  "4"))

dat$cond <- factor(dat$cond, levels=
                      c("1",
                        "2",
                        "3",
                        "4"))

Я провел ANOVA, используя ezANOVA 

aov <- ezANOVA(data=dat,
        dv=score,
        wid=subj,
        within=trial, 
        between=cond,
        type=3, # Type III as the data is unbalanced
        detailed=T 
)

Кто-нибудь знает, как я могу проводить внутриконтрастные контрасты в отношении того, что в этом случае будет пробным? Либо «вручную», либо с помощью пакета? В исходных данных важны как основные эффекты, так и взаимодействие.

Я уже смотрел на Использование R для вычисления между-внутри-subj контрастами взаимодействия ANOVA с использованием автомобиля или lme , но это мне не помогло.

Большое спасибо!

...