daDoctoR/R/rep_glm.R

#' A repeated logistic regression function
#'
#' @description For bivariate analyses. The confint() function is rather slow, causing the whole function to hang when including many predictors and calculating the ORs with CI.
#' @param y Effect meassure.
#' @param vars variables in model. Input as c() of columnnames, use dput().
#' @param string variables to test. Input as c() of columnnames, use dput().
#' @keywords logistic regression
#' @export
#' @examples
#' rep_glm()

rep_glm<-function(y,vars,string,ci=FALSE,data){
  ## x is data.frame of predictors, y is vector of an aoutcome as a factor
  ## output is returned as coefficient, or if or=TRUE as OR with 95 % CI.
  ##

  require(dplyr)

  d<-data
  x<-select(d,one_of(c(string)))
  v<-select(d,one_of(c(vars)))
  dt<-cbind(y,v)
  m1<-length(coef(glm(y~.,family = binomial(),data = dt)))

  if (!is.factor(y)){stop("Some kind of error message would be nice, but y should be a factor!")}

  if (ci==TRUE){

    df<-data.frame(matrix(ncol = 3))
    names(df)<-c("pred","or_ci","pv")

    for(i in 1:ncol(x)){

      dat<-cbind(dt,x[,i])

      m<-glm(y~.,family = binomial(),data=dat)

      l<-suppressMessages(round(exp(confint(m))[-c(1:m1),1],2))
      u<-suppressMessages(round(exp(confint(m))[-c(1:m1),2],2))
      or<-round(exp(coef(m))[-c(1:m1)],2)

      or_ci<-paste0(or," (",l," to ",u,")")

      pv<-round(tidy(m)$p.value[-c(1:m1)],3)

      x1<-x[,i]

      if (is.factor(x1)){
        pred<-paste(names(x)[i],levels(x1)[-1],sep = "_")
      }

      else {pred<-names(x)[i]}

      df<-rbind(df,cbind(pred,or_ci,pv))

    }}

  if (ci==FALSE){

    df<-data.frame(matrix(ncol = 3))
    names(df)<-c("pred","b","pv")

    for(i in 1:ncol(x)){
      dat<-cbind(dt,x[,i])

      m<-glm(y~.,family = binomial(),data=dat)

      b<-round(coef(m)[-c(1:m1)],3)

      pv<-round(tidy(m)$p.value[-c(1:m1)],3)

      x1<-x[,i]

      if (is.factor(x1)){
        pred<-paste(names(x1)[i],levels(x1)[-1],sep = "_")
      }

      else {pred<-names(x)[i]}

      df<-rbind(df,cbind(pred,b,pv))

    }}

  pa<-as.numeric(df[,3])
  pa<-ifelse(pa<0.001,"<0.001",pa)

  t <- ifelse(pa<=0.1|pa=="<0.001","include","drop")

  pa <- ifelse(pa<=0.05|pa=="<0.001",paste0("*",pa),
               ifelse(pa>0.05&pa<=0.1,paste0(".",pa),pa))

  r<-data.frame(df[,1:2],pa,t)[-1,]

  return(r)
}