inst/app/proto_widget_gen_psa_samp.R

In syzoekao/CEAutil: Utility function for Cost-effectiveness analysis at UMN

library(dampack)


dist_all <- c("normal", "log-normal", "truncated-normal",
              "beta", "gamma", "dirichlet", "bootstrap", "constant")

all_param_types <- c("mean, sd", "a, b", "shape, scale",
                     "value, mean_prop, sd", "value, n",
                     "value, alpha", "mean, sd, ll, ul", "val", "meanlog, sdlog")
all_param_types <- unique(strsplit(paste(all_param_types, collapse = ', '), ", ")[[1]])


df <- matrix(NA, ncol = length(all_param_types) + 2, nrow = 1)
df <- as.data.frame(df)
colnames(df) <- c("param", "dist", all_param_types)
df$param <- "par1"
df$dist <- "dirichlet"
df$alpha <- 0.2
df$value <- 4


gen_psa_samp_code(df, n_samp = 100)


library(dampack)
l_params_all <- gen_psa_samp(params = c("par2", "par3", "par4", "par4"),
                             dists = c("beta", "gamma", "log-normal", "truncated-normal"),
                             parameterization_types = c("a, b", "mean, sd", "meanlog, sdlog", "mean, sd, ll, ul"),
                             dists_params = list(c(3, 2), c(3, 4), c(1.3, 1), c(5, 2, 1, NA)),
                             nsamp = 100)