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

Question

Tengo dos modelos de regresión logística en R hechos con glm() . Ambos utilizan las mismas variables, pero se hicieron utilizando diferentes subconjuntos de una matriz. ¿Existe una manera fácil de obtener un modelo promedio que dé las medias de los coeficientes y luego utilizarlo con la función predict()?

[Lo siento, si este tipo de pregunta debe ser publicada en un sitio de programación, hágamelo saber y lo publicaré allí].

Gracias

Answer 1

¿Quiere tomar la media de las probabilidades previstas o la media 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 cualquiera de las dos cosas sería algo así Establecer avg a "prob" para conseguir lo primero, o algo más para lo 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)
}

Answer 2

En R, el caret tiene algunas funciones interesantes para combinar modelos.