Как я могу протестировать эффекты в ANOVA с разделенным графиком, используя подходящие сравнения моделей для использования с аргументами X
и в R? Я знаком с Далгаардом и (2007) [1]. К сожалению, это только чистит дизайн Split-Plot. Делать это в полностью рандомизированной схеме с двумя факторами внутри субъекта:M
anova.mlm()
?anova.mlm
N <- 20 # 20 subjects total
P <- 3 # levels within-factor 1
Q <- 3 # levels within-factor 2
DV <- matrix(rnorm(N* P*Q), ncol=P*Q) # random data in wide format
id <- expand.grid(IVw1=gl(P, 1), IVw2=gl(Q, 1)) # intra-subjects layout of data matrix
library(car) # for Anova()
fitA <- lm(DV ~ 1) # between-subjects design: here no between factor
resA <- Anova(fitA, idata=id, idesign=~IVw1*IVw2)
summary(resA, multivariate=FALSE, univariate=TRUE) # all tests ...
Следующие сравнения моделей приводят к тем же результатам. Ограниченная модель не включает рассматриваемый эффект, но все другие эффекты того же порядка или ниже, полная модель добавляет рассматриваемый эффект.
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw2, test="Spherical") # IVw1
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw1, test="Spherical") # IVw2
anova(fitA, idata=id, M=~IVw1 + IVw2 + IVw1:IVw2,
X=~IVw1 + IVw2, test="Spherical") # IVw1:IVw2
Дизайн Split-Splot с одним фактором внутри и одним фактором:
idB <- subset(id, IVw2==1, select="IVw1") # use only first within factor
IVb <- gl(2, 10, labels=c("A", "B")) # between-subjects factor
fitB <- lm(DV[ , 1:P] ~ IVb) # between-subjects design
resB <- Anova(fitB, idata=idB, idesign=~IVw1)
summary(resB, multivariate=FALSE, univariate=TRUE) # all tests ...
Это anova()
команды для копирования тестов, но я не знаю, почему они работают. Почему тесты следующих сравнений моделей дают одинаковые результаты?
anova(fitB, idata=idB, X=~1, test="Spherical") # IVw1, IVw1:IVb
anova(fitB, idata=idB, M=~1, test="Spherical") # IVb
Два фактора внутри субъекта и один фактор между субъектами:
fitC <- lm(DV ~ IVb) # between-subjects design
resC <- Anova(fitC, idata=id, idesign=~IVw1*IVw2)
summary(resC, multivariate=FALSE, univariate=TRUE) # all tests ...
Как мне воспроизвести результаты, приведенные выше, с соответствующими модельными сравнениями для использования с аргументами X
и ? Какова логика этих сравнений моделей?M
anova.mlm()
РЕДАКТИРОВАТЬ: Suncoolsu указал, что для всех практических целей данные из этих проектов должны быть проанализированы с использованием смешанных моделей. Тем не менее, я все еще хотел бы понять, как повторить результаты summary(Anova())
с anova.mlm(..., X=?, M=?)
.
[1]: Далгаард, П. 2007. Новые функции для многомерного анализа. Р Новости, 7 (2), 2-7.
lme4
пакет, чтобы соответствовать модели, а неlm
. Но это может быть очень конкретное книжное представление. Я позволю другим комментировать это. Я могу привести пример, основанный на том, как я его интерпретирую, который отличается от вашего.