R海拾遗——批量回归

glm_re <- function(gs, data) {
    re<-glm(as.formula(gs),data=data,family=binomial())
    or <- exp(coef(re))[2]
    zhixin<- exp(confint(re))
    or_low <- zhixin[2,1]
    or_high <- zhixin[2,2]
    result <- data.frame(or,or_low,or_high)
    return(result)
}

# 其中sufi为公式的后半部分
make_gs <- function(sufi) {
    gs_ <- c(
        paste0("y1",sufi),
        paste0("y2",sufi),
        paste0("y3",sufi),
        paste0("y4",sufi),
        paste0("y5",sufi),
        paste0("y6",sufi)
    )
return(gs_)
}

log_re <- function(gs_,data) {
    r1 <- glm_re(gs_[1], data)
    r2 <- glm_re(gs_[2], data)
    r3 <- glm_re(gs_[3], data)
    r4 <- glm_re(gs_[4], data)
    r5 <- glm_re(gs_[5], data)
    r6 <- glm_re(gs_[6], data)
    result3 <- rbind(r1,r2,r3,r4,r5,r6)
return(result3)
}