Я пытаюсь провести анализ анова в R для набора данных с одним внутренним фактором и одним между факторами. Данные взяты из эксперимента по проверке сходства двух методов тестирования. Каждый субъект был протестирован в Методе 1 и Методе 2 (внутренний фактор), а также входил в одну из 4 разных групп (междуфактор). Я пробовал использовать функции aov, Anova (в автомобильном пакете) и ezAnova. Я получаю неправильные значения для каждого метода, который я пробую. Я не уверен, в чем моя ошибка, если это непонимание R или самой Anova. Я включил код, который использовал, и считаю, что он должен работать. Я пробовал массу вариантов этого, надеясь наткнуться на ответ. Этот набор данных сбалансирован, но у меня много подобных наборов данных, и многие из них не сбалансированы. Спасибо за любую помощь, которую вы можете предоставить.
library(car)
library(ez)
#set up data
sample_data <- data.frame(Subject=rep(1:20,2),Method=rep(c('Method1','Method2'),each=20),Level=rep(rep(c('Level1','Level2','Level3','Level4'),each=5),2))
sample_data$Result <- c(4.76,5.03,4.97,4.70,5.03,6.43,6.44,6.43,6.39,6.40,5.31,4.54,5.07,4.99,4.79,4.93,5.36,4.81,4.71,5.06,4.72,5.10,4.99,4.61,5.10,6.45,6.62,6.37,6.42,6.43,5.22,4.72,5.03,4.98,4.59,5.06,5.29,4.87,4.81,5.07)
sample_data[, 'Subject'] <- as.factor(sample_data[, 'Subject'])
#Set the contrats if needed to run type 3 sums of square for unblanaced data
#options(contrats=c("contr.sum","contr.poly"))
#With aov method as I understand it 'should' work
anova_aov <- aov(Result ~ Method*Level + Error(Subject/Method),data=test_data)
print(summary(anova_aov))
#ezAnova method,
anova_ez = ezANOVA(data=sample_data, wid=Subject, dv = Result, within = Method, between=Level, detailed = TRUE, type=3)
print(anova_ez)
Кроме того, значения, которые я должен получать на выходе SAS