Tôi gặp vấn đề khi tìm giải pháp liên quan đến cách chạy bài kiểm tra sau đại học (Tukey HSD) sau khi đo 2 lần (cả hai môn học) ANOVA trong R. Đối với ANOVA, tôi đã sử dụng chức năng aov:
summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1))
Sau khi đọc câu trả lời cho các câu hỏi khác, tôi đã tập hợp rằng trước tiên tôi sẽ phải chạy lại ANOVA bằng cách sử dụng một số tính năng khác (ví dụ: lme). Đây là những gì tôi đã đưa ra.
Lme.mod <- lme(dv ~ x1*x2, random=list(subject=pdBlocked(list(~1, pdIdent(~x1-1), pdIdent(~x2-1)))), data=df1)
anova(Lme.mod)
Cả hai hiệu ứng chính đều có ý nghĩa, nhưng không có hiệu ứng tương tác. Sau đó, tôi đã sử dụng các hàm này để so sánh sau hoc:
summary(glht(Lme.mod, linfct=mcp(x1="Tukey")))
summary(glht(Lme.mod, linfct=mcp(x2="Tukey")))
Tuy nhiên, có một số vấn đề:
Trước hết, tệp Trợ giúp R nói rằng "Hàm mcp phải được sử dụng cẩn thận khi xác định các tham số quan tâm trong các mô hình đa chiều ANOVA hoặc ANCOVA hai chiều (...) phiên bản 1.0-0 và cao hơn tạo ra các so sánh cho các hiệu ứng chính chỉ, bỏ qua các đồng biến và tương tác (các phiên bản cũ hơn tự động tính trung bình theo các điều khoản tương tác). Một cảnh báo được đưa ra. " Và chắc chắn, tôi đã nhận được thông báo cảnh báo sau:
Warning message:
In mcp2matrix(model, linfct = linfct) :
covariate interactions found -- default contrast might be inappropriate
Một điều khó hiểu nữa là mặc dù cả hai tác dụng chính đều có ý nghĩa, nhưng không có sự khác biệt đáng kể nào trong các so sánh sau hoc đối với một trong các yếu tố (x1). Tôi chưa bao giờ gặp phải điều này trước đây. Là các kịch bản / phân tích chính xác / phù hợp, hoặc có điều gì tôi đang thiếu? Bất kỳ trợ giúp sẽ được đánh giá cao nhất!