improved function summary all to match true and estimated parameters.

R/sim_all.R

@@ -169,7 +169,8 @@ sim_all <- function(nosim = 2,
   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::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"

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")

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"