From b432a950170173d1e1a9ac216485901f4cfa6595 Mon Sep 17 00:00:00 2001 From: Julian Sagebiel <julian.sagebiel@idiv.de> Date: Sun, 29 Dec 2024 14:10:07 +0100 Subject: [PATCH] improved function summary all to match true and estimated parameters. --- R/sim_all.R | 19 ++++++++++--------- tests/manual-tests/SE_Drive.R | 2 +- tests/manual-tests/feedadditives.R | 2 +- 3 files changed, 12 insertions(+), 11 deletions(-) diff --git a/R/sim_all.R b/R/sim_all.R index 9893ba2..c2ff269 100644 --- a/R/sim_all.R +++ b/R/sim_all.R @@ -165,15 +165,16 @@ sim_all <- function(nosim = 2, - 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"), - by="parname") %>% - dplyr::relocate(parname, dplyr::ends_with(c( - ".n", "truepar", "mean", "sd", "min" , "max", "range" , "se" - ))) +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) diff --git a/tests/manual-tests/SE_Drive.R b/tests/manual-tests/SE_Drive.R index 8943e7b..9b844e3 100644 --- a/tests/manual-tests/SE_Drive.R +++ b/tests/manual-tests/SE_Drive.R @@ -2,7 +2,7 @@ rm(list=ls()) devtools::load_all() - +library(rlang) designpath<- system.file("extdata","SE_DRIVE" ,package = "simulateDCE") diff --git a/tests/manual-tests/feedadditives.R b/tests/manual-tests/feedadditives.R index 23f91a7..7d54a51 100644 --- a/tests/manual-tests/feedadditives.R +++ b/tests/manual-tests/feedadditives.R @@ -7,7 +7,7 @@ designpath<- system.file("extdata","feedadditives" ,package = "simulateDCE") notes <- "Three heuristics" -resps =300 # number of respondents +resps =396 # number of respondents nosim=2 # number of simulations to run (about 500 is minimum) destype = "ngene" -- GitLab