Chisquare on variables from a dataframe and give gt_summary output with all chi-square parameters including standardized residuals

chisq <- function(data, i, j, chisq_flag=1, verbose=1){
  if(verbose == 1){
    print(paste("Chisquare test between", i, "and", j))
  }
  
  i_factor <- as.factor(data[[i]])
  j_factor <- as.factor(data[[j]])
  
  if (chisq_flag == 1){
    chisq_results <- chisq.test(j_factor, i_factor)
  
    if(verbose == 1){
      print(chisq_results$p.value)
    }
  } else {
    chisq_results <- chisq.test(j_factor, i_factor, correct=FALSE)
    
    if(verbose == 1){
      print(chisq_results$p.value)
    }
  }
  
  gt_results <- gt_summary(j_factor, i_factor)
  
  print(gt_results)