R/parcomparisons.R

In danilinares/quickpsy: Fits Psychometric Functions for Multiple Groups

#' Pair comparisons of the parameters using bootstrap
#' \code{parcomparisons} Calculates the bootstrap confidence intervals for the
#' difference in the parameters for two groups for all possible pairs
#' of groups
#' @keywords internal
#' @importFrom rlang .data
parcomparisons <- function(par_dif, par_difbootstrap, ci) {

  ci <- par_difbootstrap %>%
    group_by(!!!syms(setdiff(names(par_dif), c("par", "par2", "dif")))) %>%
    summarise(difinf = quantile(.data$dif, .5*(1 - ci))[[1]],
              difsup = quantile(.data$dif, 1 - .5*(1 - ci))[[1]],
              signif = ifelse(.data$difinf * .data$difsup < 0, "", "*"), .groups = "keep")

  par_dif %>% left_join(ci, by = setdiff(names(par_dif), c("par", "par2", "dif")))
}