¿Hay una manera fácil de combinar dos modelos glm en R?

Tengo dos modelos de regresión logística en R hechos con glm(). Ambos usan las mismas variables, pero se hicieron usando diferentes subconjuntos de una matriz. ¿Hay una manera fácil de obtener un modelo promedio que proporcione las medias de los coeficientes y luego usar esto con la función predic ()?

[perdón si este tipo de pregunta debe publicarse en un sitio de programación, hágamelo saber y lo publicaré allí]

Gracias

¿Desea tomar el promedio de las probabilidades pronosticadas o el promedio de los coeficientes? Darán resultados diferentes, porque una regresión logística implica una transformación no lineal del predictor lineal.

Una función para hacer sería algo como esto. Establezca avgpara "prob"obtener el primero, o algo más para el segundo.

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