Làm cách nào tôi có thể kiểm tra các hiệu ứng trong ANOVA phân chia lô bằng cách sử dụng các so sánh mô hình phù hợp để sử dụng với các đối số X
và trong R? Tôi quen thuộc với và Dalgaard (2007) [1]. Thật không may, nó chỉ chải thiết kế Split-Plot. Làm điều này trong một thiết kế hoàn toàn ngẫu nhiên với hai yếu tố bên trong chủ đề: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 ...
Các so sánh mô hình sau đây dẫn đến kết quả tương tự. Mô hình bị hạn chế không bao gồm hiệu ứng trong câu hỏi nhưng tất cả các hiệu ứng khác có cùng thứ tự hoặc thấp hơn, mô hình đầy đủ sẽ thêm hiệu ứng trong câu hỏi.
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
Thiết kế Split-Splot với một yếu tố bên trong và một giữa các chủ thể:
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 ...
Đây là các anova()
lệnh để sao chép các bài kiểm tra, nhưng tôi không biết tại sao chúng lại hoạt động. Tại sao các thử nghiệm so sánh mô hình sau đây dẫn đến kết quả tương tự?
anova(fitB, idata=idB, X=~1, test="Spherical") # IVw1, IVw1:IVb
anova(fitB, idata=idB, M=~1, test="Spherical") # IVb
Hai yếu tố bên trong chủ thể và một yếu tố giữa các chủ thể:
fitC <- lm(DV ~ IVb) # between-subjects design
resC <- Anova(fitC, idata=id, idesign=~IVw1*IVw2)
summary(resC, multivariate=FALSE, univariate=TRUE) # all tests ...
Làm cách nào để sao chép các kết quả đã cho ở trên với các so sánh mô hình tương ứng để sử dụng với X
và các M
đối số của anova.mlm()
? Logic đằng sau những so sánh mô hình này là gì?
EDIT: suncoolsu chỉ ra rằng cho tất cả các mục đích thực tế, dữ liệu từ các thiết kế này nên được phân tích bằng các mô hình hỗn hợp. Tuy nhiên, tôi vẫn muốn hiểu làm thế nào để sao chép kết quả summary(Anova())
với anova.mlm(..., X=?, M=?)
.
[1]: Dalgaard, P. 2007. Các chức năng mới để phân tích đa biến. Tin tức R, 7 (2), 2-7.
lme4
gói để phù hợp với mô hình VÀ KHÔNGlm
. Nhưng đây có thể là một quan điểm dựa trên cuốn sách rất cụ thể. Tôi sẽ để người khác nhận xét về nó. Tôi có thể đưa ra một ví dụ dựa trên cách tôi diễn giải nó khác với bạn.