Czy istnieje prosty sposób na połączenie dwóch modeli GLM w R?

Mam dwa modele regresji logistycznej w R wykonane glm(). Oba używają tych samych zmiennych, ale zostały wykonane przy użyciu różnych podzbiorów macierzy. Czy istnieje prosty sposób na uzyskanie średniego modelu, który podaje średnie współczynników, a następnie użycie go z funkcją predykcji ()?

[przepraszam, jeśli tego typu pytanie powinno zostać opublikowane na stronie programistycznej, daj mi znać, a tam go opublikuję]

Dzięki

Czy chcesz wziąć średnią z przewidywanych prawdopodobieństw, czy średnią współczynników? Dadzą różne wyniki, ponieważ regresja logistyczna obejmuje nieliniową transformację predyktora liniowego.

Każda z tych funkcji byłaby czymś takim. Ustaw, avgaby "prob"uzyskać pierwsze, lub coś innego dla drugiego.

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