Kiểm tra sự tương đương của các mô hình không lồng nhau


12

Giả sử là hàm tuyến tính của x và một hình nộm d . Giả thuyết của tôi là d chính nó là giống như một chỉ số chủ nghĩa khoái lạc của một vector của các biến khác, Z . Tôi có hỗ trợ cho điều này trong M A N O V A của Z (tức là z 1 , z 2 , ..., z n ) trên d . Có cách nào để kiểm tra sự tương đương của hai mô hình này không:yxddZMMộtNÔiVMộtZz1z2znd

Mô hình 1: y= =b0+b1x+b2d+e1

Mô hình 2: y= =g0+ZG+e2

Trong đó là vectơ cột của các tham số.G

Câu trả lời:


8

Để bắt đầu với bạn phải xác định khái niệm tương đương . Người ta có thể nghĩ rằng hai mô hình tương đương nhau khi chúng tạo ra độ chính xác dự báo gần như nhau (mô hình này có liên quan đến chuỗi thời gian và dữ liệu bảng điều khiển), một mô hình khác có thể quan tâm nếu phù hợp với mô hình gần . Cái trước là đối tượng để xác thực chéo khác nhau (thường dùng dao hoặc một số thử nghiệm ngoài mẫu, Rob accuracy()làm điều này một cách độc đáo), cái sau dùng để giảm thiểu một số tiêu chí thông tin.

BTôiCMộtTôiC

Một cuộc thảo luận thú vị được đưa ra trong cuốn sách bắt buộc của Cameron và Trivingi (Chương 8.5 cung cấp đánh giá xuất sắc về các phương pháp), chi tiết lý thuyết cụ thể hơn được tìm thấy ở Hong và Preston tại đây .

Nói một cách đơn giản, lựa chọn từ hai mô hình càng kỹ lưỡng hơn (có ít tham số để ước tính, do đó, nhiều mức độ tự do hơn) sẽ được đề xuất là thích hợp hơn. Một tiêu chí thông tin giới thiệu một hàm hình phạt đặc biệt hạn chế việc đưa các biến giải thích bổ sung vào mô hình tuyến tính về mặt khái niệm tương tự như các hạn chế được đưa ra bởi điều chỉnhR2.

Tuy nhiên, bạn có thể không chỉ quan tâm đến việc chọn mô hình giảm thiểu tiêu chí thông tin đã chọn. Khái niệm tương đương ngụ ý rằng một số thống kê kiểm tra nên được xây dựng. Do đó, bạn có thể đi kiểm tra tỷ lệ khả năng hoặc Cox hoặc VoungLR thử nghiệm, Davidson-MacKinnon J kiểm tra.

Cuối cùng, theo các thẻ, bạn có thể chỉ quan tâm đến các Rchức năng:

library(lmtest)
coxtest(fit1, fit2)
jtest(fit1, fit2)

Trường hợp fit1fit2là hai mô hình hồi quy tuyến tính được trang bị không lồng nhau, coxtestlà CoxLRthử nghiệm và jtestDavidson-MacKinnonJ kiểm tra.


Cảm ơn, Dmitrij. Nếu tôi hiểu chính xác, cả coxtest và jtest về cơ bản đều là các thử nghiệm lồng nhau được sửa đổi. Bước 1: Chạy mô hình với nhóm hồi quy kết hợp từ model1 và model2. Bước 2: Kiểm tra riêng từng model1 và model2 dưới dạng các tập hợp con của "siêu mẫu". Tôi có đúng không Ngoài ra, trên lưu ý về các biện pháp IC, có cách nào để so sánh thống kê sự khác biệt của AIC / BIC giữa các mô hình 1 và 2 không? Lưu ý: Tôi KHÔNG cố chọn mô hình "tốt nhất", nhưng bạn nói đúng rằng về cơ bản tôi đang thử kiểm tra xem hai mô hình có giống nhau không.
dùng3671

@user, bạn không cần thêm vào siêu mẫu bất cứ điều gì, chỉ cần cung cấp jtesthoặc coxtestvới sự phù hợp không lồng nhau từ Bước 1. Tiêu chí thông tin cho không lồng nhau sẽ là một hướng dẫn tốt cho mô hình nào phù hợp hơn về mặt thống kê (phân tích), nhưng đối với thử nghiệm giả thuyết tôi sẽ chỉ đi cho bất kỳLR(thực sự khả năng đăng nhập là một phần của bất kỳ tiêu chí thông tin nào) tự kiểm tra. Kết luận sẽ có phần gần gũi, nhưng vì có hai chức năng hình phạt được đưa ra một cách xác định nên sẽ hơi khó để so sánh chúng theo thống kê.
Dmitrij Celov
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.