R/jamovi_mdiff_functions.R

In esci: Estimation Statistics with Confidence Intervals

jamovi_mdiff_initialize <- function(self, grouping_variable = TRUE) {

  # Set some variables for convenience -----------------------
  #   Is analysis from summary data or raw?
  #   Are we evaluating a hypothesis?
  #   Is this a contrast?
  from_raw <- (self$options$switch == "from_raw")

  # Get a handle for each table
  tbl_overview <- NULL
  tbl_es_mean_difference <- NULL
  tbl_es_mean_ratio <- NULL
  tbl_es_smd <- NULL
  tbl_es_median_difference <- NULL
  tbl_es_median_ratio <- NULL
  tbl_es_odds_ratio <- NULL
  tbl_hypothesis_evaluations <- NULL
  # tbl_htest_summary <- NULL
  assume_equal_variance <- NULL
  try(tbl_overview <- self$results$overview, silent = TRUE)
  try(tbl_es_mean_difference <- self$results$es_mean_difference, silent = TRUE)
  try(tbl_es_mean_ratio <- self$results$es_mean_ratio, silent = TRUE)
  try(tbl_es_smd <- self$results$es_smd, silent = TRUE)
  try(tbl_es_median_difference <- self$results$es_median_difference, silent = TRUE)
  try(tbl_es_median_ratio <- self$results$es_median_ratio, silent = TRUE)
  try(tbl_es_odds_ratio <- self$results$es_odds_ratio, silent = TRUE)
  try(assume_equal_variance <- self$options$assume_equal_variance, silent = TRUE)
  try(tbl_point_null <- self$results$point_null, silent = TRUE)
  try(tbl_interval_null <- self$results$interval_null, silent = TRUE)
  try(effect_size <- self$options$effect_size, silent = TRUE)


  if (effect_size == "mean_difference") {

    if (!is.null(tbl_overview) & !is.null(assume_equal_variance)) {
      if (assume_equal_variance) {
        tbl_overview$setNote(
          key = "overview_table",
          note = "Variances are assumed equal, so <i>s</i><sub>p</sub> was used to calculate each CI.",
          init = TRUE
        )
      } else {
        tbl_overview$setNote(
          key = "overview_table",
          note = "Variances are not assumed equal, and so the CI was calculated separately for each mean.",
          init = TRUE
        )

      }
    }
  } else {
    tbl_overview$setNote(
      key = "overview_table",
      note = NULL,
      init = TRUE
    )
  }

  if (!is.null(tbl_es_mean_difference) & !is.null(assume_equal_variance)) {
    if (assume_equal_variance) {
      tbl_es_mean_difference$setNote(
        key = "overview_table",
        note = "Variances are assumed equal, so <i>s</i><sub>p</sub> was used to calculate each CI.",
        init = FALSE
      )
    } else {
      tbl_es_mean_difference$setNote(
        key = "overview_table",
        note = "Variances are not assumed equal, so the Welch method was used to calculate each CI on a difference.",
        init = FALSE
      )

    }
  }

  # Prep output -------------------------------------------
  # Set CI and MoE columns to reflect confidence level
  conf_level <- jamovi_sanitize(
    my_value = self$options$conf_level,
    return_value = 95,
    na_ok = FALSE,
    convert_to_number = TRUE,
    lower = 75,
    lower_inclusive = FALSE,
    upper = 100,
    upper_inclusive = FALSE,
    my_value_name = "Confidence level"
  )

  jamovi_set_confidence(tbl_overview, conf_level)
  jamovi_set_confidence(tbl_es_mean_difference, conf_level)
  jamovi_set_confidence(tbl_es_smd, conf_level)
  jamovi_set_confidence(tbl_es_mean_ratio, conf_level)
  jamovi_set_confidence(tbl_es_median_difference, conf_level)
  jamovi_set_confidence(tbl_es_median_ratio, conf_level)


  # Outcomes: 1 if from summary, length of outcome_variables if raw
  outcome_count <- if(from_raw) {
    length(self$options$outcome_variable)
  } else {
    1
  }

  # For now, only 1 contrast can be specified
  contrast_count <- 1

  # How many levels?
  #  For raw, check grouping_variable
  #  For summary, check group_labels
  if (grouping_variable) {
    if (from_raw) {
      if (!is.null(self$options$grouping_variable)) {
        level_source <- self$options$grouping_variable
        level_count <- length(levels(as.factor(self$data[, level_source])))
      } else {
        level_count = 2
      }
    } else {
      level_source <- try(level_source <- self$options$grouping_variable_levels, silent = TRUE)
      if (is(level_source, "try-error")) {
        level_count = 2
      } else {
        level_count <- length(levels(as.factor(self$data[, level_source])))
      }
    }
  } else {
    level_count <- 1
  }



  # Rows needed for each table -------------------------------
  overview_rows <- level_count * outcome_count
  mdiff_rows <- contrast_count * outcome_count * 3
  smd_rows <- contrast_count * outcome_count

  jamovi_init_table(tbl_overview, overview_rows)
  jamovi_init_table(tbl_es_mean_difference, mdiff_rows, breaks = 3)
  jamovi_init_table(tbl_es_smd, smd_rows)
  jamovi_init_table(tbl_es_mean_ratio, smd_rows)
  jamovi_init_table(tbl_es_median_difference, mdiff_rows, breaks = 3)
  jamovi_init_table(tbl_es_median_ratio, smd_rows, breaks = 3)

  if (!is.null(tbl_point_null)) {
    jamovi_init_table(
      tbl_point_null,
      outcome_count
    )
  }
  if (!is.null(tbl_interval_null)) {
    jamovi_init_table(
      tbl_interval_null,
      outcome_count
    )
  }


  width <- jamovi_sanitize(
    my_value = self$options$es_plot_width,
    return_value = 600,
    convert_to_number = TRUE,
    lower = 10,
    lower_inclusive = TRUE,
    upper = 3000,
    upper_inclusive = TRUE
  )
  height <- jamovi_sanitize(
    my_value = self$options$es_plot_height,
    return_value = 400,
    convert_to_number = TRUE,
    lower = 10,
    lower_inclusive = TRUE,
    upper = 4000,
    upper_inclusive = TRUE
  )

  keys <- if (from_raw)
    self$options$outcome_variable
  else
    jamovi_sanitize(
      self$options$outcome_variable_name,
      "My outcome variable",
      na_ok = FALSE
    )

  for (my_key in keys) {
    self$results$estimation_plots$addItem(key = my_key)
    image <- self$results$estimation_plots$get(my_key)
    image$setSize(width , height)
  }


}


jamovi_add_htest_mdiff <- function(self, estimate) {
  evaluate_h <- self$options$evaluate_hypotheses

  if(evaluate_h) {
    # Test results
    effect_size = self$options$effect_size
    if (effect_size == "mean_difference") effect_size <- "mean"
    if (effect_size == "median_difference") effect_size <- "median"

    rope_upper <- jamovi_sanitize(
      self$options$null_boundary,
      na_ok = FALSE,
      return_value = 0,
      lower = 0,
      lower_inclusive = TRUE,
      convert_to_number = TRUE
    )

    rope_units <- "raw"
    try(rope_units <- self$options$rope_units, silent = TRUE)

    test_results <- test_mdiff(
      estimate,
      effect_size = effect_size,
      rope = c(rope_upper * -1, rope_upper),
      rope_units = rope_units,
      output_html = TRUE
    )

    estimate$point_null <- test_results$point_null
    estimate$interval_null <- test_results$interval_null

    if (!is.null(names(rope_upper))) {
      self$results$point_null$setVisible(TRUE)
      self$results$interval_null$setVisible(FALSE)
    }

  }

  return(estimate)

}