Apakah ada cara mudah untuk menggabungkan dua model glm dalam R?

Saya memiliki dua model regresi logistik dalam R yang dibuat dengan glm(). Keduanya menggunakan variabel yang sama, tetapi dibuat menggunakan subset matriks yang berbeda. Apakah ada cara mudah untuk mendapatkan model rata-rata yang memberikan nilai koefisien dan kemudian menggunakannya dengan fungsi predict ()?

[maaf jika pertanyaan jenis ini harus diposting di situs pemrograman, beri tahu saya dan saya akan mengirimkannya di sana]

Terima kasih

Apakah Anda ingin mengambil rata-rata probabilitas yang diprediksi, atau rata-rata koefisien? Mereka akan memberikan hasil yang berbeda, karena regresi logistik melibatkan transformasi nonlinier dari prediktor linier.

Fungsi untuk melakukan keduanya adalah seperti ini. Tetapkan avguntuk "prob"mendapatkan yang pertama, atau yang lain untuk yang terakhir.

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)
}