C'è un modo semplice per combinare due modelli glm in R?

Ho due modelli di regressione logistica in R realizzati con glm(). Entrambi usano le stesse variabili, ma sono stati realizzati usando diversi sottoinsiemi di una matrice. Esiste un modo semplice per ottenere un modello medio che fornisce la media dei coefficienti e quindi utilizzarlo con la funzione predict ()?

[scusate se questo tipo di domanda dovrebbe essere pubblicato su un sito di programmazione fatemi sapere e lo pubblicherò lì]

Grazie

Vuoi prendere la media delle probabilità previste o la media dei coefficienti? Daranno risultati diversi, perché una regressione logistica comporta una trasformazione non lineare del predittore lineare.

Una funzione da svolgere sarebbe qualcosa del genere. Imposta avgper "prob"ottenere il primo o qualcos'altro per il secondo.

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