Gibt es eine einfache Möglichkeit, zwei glm-Modelle in R zu kombinieren?

Ich habe zwei logistische Regressionsmodelle in R erstellt mit glm(). Beide verwenden dieselben Variablen, wurden jedoch mit unterschiedlichen Teilmengen einer Matrix erstellt. Gibt es eine einfache Möglichkeit, ein Durchschnittsmodell zu erhalten, das die Mittelwerte der Koeffizienten angibt, und dieses dann mit der Funktion Predict () zu verwenden?

[Entschuldigung, wenn diese Art von Frage auf einer Programmierseite veröffentlicht werden soll, lass es mich wissen und ich werde sie dort posten]

Vielen Dank

Möchten Sie den Durchschnitt der vorhergesagten Wahrscheinlichkeiten oder den Durchschnitt der Koeffizienten nehmen? Sie liefern unterschiedliche Ergebnisse, da eine logistische Regression eine nichtlineare Transformation des linearen Prädiktors beinhaltet.

Eine Funktion, die beides ausführen könnte, wäre ungefähr so. Stellen Sie ein avg, "prob"um das erstere oder etwas anderes für das letztere zu erhalten.

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