Có cách nào dễ dàng để kết hợp hai mô hình glm trong R không?

Tôi có hai mô hình hồi quy logistic trong R được thực hiện với glm(). Cả hai đều sử dụng cùng một biến, nhưng được tạo bằng các tập con khác nhau của ma trận. Có một cách dễ dàng để có được một mô hình trung bình cung cấp các phương tiện của các hệ số và sau đó sử dụng mô hình này với hàm dự đoán () không?

[xin lỗi nếu loại câu hỏi này nên được đăng trên một trang web lập trình cho tôi biết và tôi sẽ đăng nó ở đó]

Cảm ơn

Bạn có muốn lấy trung bình của xác suất dự đoán, hoặc trung bình của các hệ số? Họ sẽ đưa ra các kết quả khác nhau, bởi vì hồi quy logistic liên quan đến một biến đổi phi tuyến của bộ dự báo tuyến tính.

Một chức năng để làm hoặc là một cái gì đó như thế này. Đặt avgđể "prob"lấy cái trước, hoặc cái gì khác cho cái sau.

pred_comb <- function(mod1, mod2, dat, avg="prob", ...)
{
    xb1 <- predict(mod1, dat, type="link", ...)
    xb2 <- predict(mod2, dat, type="link", ...)
    if(avg == "prob")
        (plogis(xb1) + plogis(xb2))/2
    else plogis((xb1 + xb2)/2)
}