Có ai có đề xuất hoặc gói sẽ tính hệ số xác định một phần?
Hệ số xác định từng phần có thể được định nghĩa là phần trăm biến thể không thể giải thích được trong mô hình rút gọn, nhưng có thể được giải thích bằng các yếu tố dự đoán được chỉ định trong mô hình (er) đầy đủ. Hệ số này được sử dụng để cung cấp cái nhìn sâu sắc về việc một hoặc nhiều dự đoán bổ sung có thể hữu ích trong mô hình hồi quy được chỉ định đầy đủ hơn hay không.
Tính toán cho một phần r ^ 2 tương đối đơn giản sau khi ước tính hai mô hình của bạn và tạo các bảng ANOVA cho chúng. Tính toán cho một phần r ^ 2 là:
(SSEreduces - SSEfull) / SSEreduces
Tôi đã viết hàm tương đối đơn giản này sẽ tính toán điều này cho mô hình hồi quy tuyến tính đa biến. Tôi không quen thuộc với các cấu trúc mô hình khác trong R nơi chức năng này có thể không hoạt động tốt:
partialR2 <- function(model.full, model.reduced){
anova.full <- anova(model.full)
anova.reduced <- anova(model.reduced)
sse.full <- tail(anova.full$"Sum Sq", 1)
sse.reduced <- tail(anova.reduced$"Sum Sq", 1)
pR2 <- (sse.reduced - sse.full) / sse.reduced
return(pR2)
}
Bất kỳ đề xuất hoặc mẹo về các chức năng mạnh mẽ hơn để thực hiện nhiệm vụ này và / hoặc triển khai hiệu quả hơn của mã trên sẽ được đánh giá cao.