Các lmerTest
gói cung cấp một anova()
chức năng cho các mô hình hỗn hợp tuyến tính với tùy chọn xấp xỉ của bậc tự do (df) Satterthwaite của (mặc định) hoặc Kenward-Roger. Sự khác biệt giữa hai phương pháp này là gì? Khi nào nên chọn cái nào?
Các lmerTest
gói cung cấp một anova()
chức năng cho các mô hình hỗn hợp tuyến tính với tùy chọn xấp xỉ của bậc tự do (df) Satterthwaite của (mặc định) hoặc Kenward-Roger. Sự khác biệt giữa hai phương pháp này là gì? Khi nào nên chọn cái nào?
Câu trả lời:
Tôi cũng quan tâm đến việc tìm ra sự khác biệt có thể là gì. Điều tốt nhất tôi có thể cung cấp cho bạn, bây giờ, là bài đăng trên blog này cho thấy rằng xấp xỉ Kenward-Roger là hơi, nhưng có lẽ không đáng kể, bảo thủ hơn so với xấp xỉ Satterthwaite. Tác giả cũng lưu ý rằng cả hai đều bảo thủ hơn so với xấp xỉ bình thường, nhưng một lần nữa, không nhiều nếu kích thước mẫu đủ cao. Tôi không chắc liệu đây có phải là một kết luận chung chung của tác giả hay không.
Chỉnh sửa: Tôi sẽ thêm rằng bài viết "So sánh mức độ của mẫu số của các phương pháp xấp xỉ tự do trong mô hình hỗn hợp hai chiều không cân bằng" của KB Gregory dường như chỉ ra rằng không có phương pháp nào thường là phương pháp tốt hơn, mặc dù rõ ràng có những trường hợp Xấp xỉ Kenward-Roger mất một số mức độ bảo thủ.
Một sự khác biệt khác giữa hai phương pháp được mô tả trong Luke (2017):
Cả hai cách tiếp cận Kenward-Roger (Kenward & Roger, 1997) và Satterthwaite (1941) đều được sử dụng để ước tính mức độ tự do của mẫu số cho thống kê F hoặc mức độ tự do cho thống kê t. SAS PROC MIXED sử dụng phép tính gần đúng Satterthwaite (Học viện SAS, 2008). Mặc dù xấp xỉ Satterthwaite có thể được áp dụng cho các mô hình ML hoặc REML, nhưng xấp xỉ Kenward-Roger chỉ được áp dụng cho các mô hình REML.
(đậm của tôi)