R에 두 개의 glm 모델을 결합하는 쉬운 방법이 있습니까?

R로 만든 두 가지 로지스틱 회귀 모델이 있습니다 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)
}