Rで2つのglmモデルを組み合わせる簡単な方法はありますか？

で作成されたRのロジスティック回帰モデルが2つありglm()ます。どちらも同じ変数を使用しますが、マトリックスの異なるサブセットを使用して作成されました。係数の平均を与える平均モデルを取得し、これをpredict（）関数で使用する簡単な方法はありますか？

[この種の質問をプログラミングサイトに投稿する必要がある場合は、お知らせください。投稿します]

ありがとう

予測された確率の平均、または係数の平均を取得しますか？ロジスティック回帰には線形予測子の非線形変換が含まれるため、結果は異なります。

どちらかを行う関数はこのようなものでしょう。前者、または後者の何かを取得するように設定avgし"prob"ます。

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