added aggegate_results_function

R/aggregateResults.R

+#' Aggregate Simulation Results
+#'
+#' Processes the simulation results to extract summaries, coefficients, and graphs.
+#'
+#' @param all_designs A list of simulation results from sim_choice.
+#' @param bcoeff A named list of true parameter values used in the simulation.
+#' @param designname A character vector of design names used in the simulation.
+#' @param reshape_type Method for reshaping data: "auto", "stats", or "tidyr".
+#'
+#' @return A list with aggregated results including summary, coefficients, graphs, and power.
+aggregateResults <- function(){
+
+
+
+powa <- purrr::map(all_designs, ~ .x$power)
+
+summaryall <- as.data.frame(purrr::map(all_designs, ~ .x$summary)) %>%
+  dplyr::select(!dplyr::ends_with("vars")) %>%
+  tibble::rownames_to_column("parname") %>%
+  dplyr::mutate(parname = stringr::str_remove(parname, "^est_")) %>%
+  dplyr::left_join(data.frame(truepar = unlist(bcoeff)) %>% tibble::rownames_to_column("parname") %>%
+                     dplyr::mutate(parname=stringr::str_replace_all(parname,"\\.","_")),   ##because parameters have been renamed for mixl, we have to make sure we substitute all . with _
+                   by="parname") %>%
+  dplyr::relocate(parname, dplyr::ends_with(c(
+    ".n", "truepar", "mean", "sd", "min" , "max", "range" , "se"
+  )))
+
+coefall <- purrr::map(all_designs, ~ .x$coefs)
+
+pat <- paste0("(", paste(designname, collapse = "|"), ").") # needed to identify pattern to be replaced
+
+
+preprocessed <- as.data.frame(coefall) %>%
+  dplyr::select(!dplyr::matches("pval|run")) %>%
+  dplyr::rename_with(~ sub("est_", "", .x), dplyr::everything()) %>%
+  dplyr::rename_with( ~ paste0(., "_", stringr::str_extract(., pat)),
+                      dplyr::everything()) %>%   # rename attributes for reshape part 1
+  dplyr::rename_with( ~ stringr::str_replace(., pattern = pat, replacement =
+                                               ""),
+                      dplyr::everything())
+
+
+s <- switch(
+  reshape_type,
+  "stats" = {
+    message("Using stats::reshape for reshaping...")
+    stats::reshape(
+      preprocessed,
+      varying = 1:ncol(preprocessed),
+      sep = "_",
+      direction = "long",
+      timevar = "design"
+    ) %>%
+      dplyr::select(-id)
+  },
+  "tidyr" = {
+    message("Using tidyr::pivot_longer for reshaping...")
+    tidyr::pivot_longer(
+      preprocessed,
+      cols = dplyr::everything(),
+      names_to = c(".value", "design"),
+      names_sep = "_",
+      values_drop_na = TRUE
+    )
+  },
+  "auto" = {
+    tryCatch(
+      {
+        message("Trying tidyr::pivot_longer for reshaping...")
+        tidyr::pivot_longer(
+          preprocessed,
+          cols = dplyr::everything(),
+          names_to = c(".value", "design"),
+          names_sep = "_",
+          values_drop_na = TRUE
+        )
+      },
+      error = function(e) {
+        message("tidyr::pivot_longer failed, falling back to stats::reshape...")
+        stats::reshape(
+          preprocessed,
+          varying = 1:ncol(preprocessed),
+          sep = "_",
+          direction = "long",
+          timevar = "design"
+        ) %>%
+          dplyr::select(-id)
+      }
+    )
+  },
+  stop("You need to specify either 'tidyr', 'stats', or 'auto' as the reshape_type")
+)
+
+
+p = list()
+
+for (att in names(dplyr::select(s, -c("design")))) {
+  p[[att]] <- plot_multi_histogram(s, att, "design")
+
+  print(p[[att]])
+
+}
+
+all_designs[["summaryall"]] = summaryall
+all_designs[["graphs"]] = p
+all_designs[["powa"]] = powa
+
+
+all_designs[["time"]] = time
+all_designs[["arguements"]] = list(
+  "Beta values" = bcoeff,
+  "Utility functions" = u ,
+  "Decision groups" = decisiongroups ,
+  "Manipulation of vars" = manipulations,
+  "Number Simulations" = nosim,
+  "Respondents" = resps,
+  "Designpath" = designpath
+)
+
+
+
+return(all_designs)
+}

man/aggregateResults.Rd

+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/aggregateResults.R
+\name{aggregateResults}
+\alias{aggregateResults}
+\title{Aggregate Simulation Results}
+\usage{
+aggregateResults()
+}
+\arguments{
+\item{all_designs}{A list of simulation results from sim_choice.}
+
+\item{bcoeff}{A named list of true parameter values used in the simulation.}
+
+\item{designname}{A character vector of design names used in the simulation.}
+
+\item{reshape_type}{Method for reshaping data: "auto", "stats", or "tidyr".}
+}
+\value{
+A list with aggregated results including summary, coefficients, graphs, and power.
+}
+\description{
+Processes the simulation results to extract summaries, coefficients, and graphs.
+}