มีวิธีง่ายๆในการรวมสองรุ่น glm ใน R หรือไม่?

glm()ฉันมีสองรูปแบบการถดถอยโลจิสติกในการวิจัยที่ทำด้วย พวกเขาทั้งสองใช้ตัวแปรเดียวกัน แต่ทำโดยใช้ชุดย่อยที่แตกต่างกันของเมทริกซ์ มีวิธีง่ายๆในการรับแบบจำลองเฉลี่ยซึ่งให้ค่าสัมประสิทธิ์แล้วใช้กับฟังก์ชันทำนาย () หรือไม่

[ขออภัยถ้าคำถามประเภทนี้ควรโพสต์บนเว็บไซต์การเขียนโปรแกรมแจ้งให้เราทราบและฉันจะโพสต์ที่นั่น]

ขอบคุณ

คุณต้องการที่จะใช้ค่าเฉลี่ยของความน่าจะเป็นที่คาดการณ์หรือค่าสัมประสิทธิ์เฉลี่ยหรือไม่? พวกเขาจะให้ผลลัพธ์ที่แตกต่างกันเนื่องจากการถดถอยโลจิสติกเกี่ยวข้องกับการแปลงแบบไม่เชิงเส้นของตัวทำนายเชิงเส้น

ฟังก์ชั่นที่จะทำอย่างใดอย่างหนึ่งจะเป็นเช่นนี้ กำหนด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)
}