[R] LHS random number generator

[Shah Alam]

Hello everyone,

I am using the LHS package to generate a combination of a
set of parameter values. However, I am not sure how to fix some parameter
values to 2 decimal points.
For example:

I want to make combinations in such a way that values for parameter
c("lmp", 0.40, 0.43) are taken as 0.40, 0.41, 0.42,0.43.

My codes are:

prior_lhs <- list(c("r_mu", 0.00299, 0.0032),
                  c("r_sd", 0.001, 0.002),
                  c("lmp", 0.40, 0.43),
                  c("gr_mu", 0.14, 0.16),
                  c("gr_sd", 0.01, 0.020),
                  c("alpha1", 0.0001, 0.0018),
                  c("alpha2", 0.0017, 0.0028),
                  c("alpha3", 0.005, 0.009),
                  c("beta", 0.69, 0.75)
)


### Simulation with priors from LHS ####
nb_simul <- 10000
nparam <- length(prior_lhs)
random_tab = randomLHS(nb_simul, nparam)
lhs_index = 1
param <- matrix(rep(0, nparam), nrow = 1, ncol = nparam)

for (i in 1:nb_simul) {
  temp_par <- matrix(rep(0, nparam), nrow = 1, ncol = nparam)
  for (j in 1:nparam) {
    temp_par[j] = as.numeric(prior_lhs[[j]][2]) +
      (as.numeric(prior_lhs[[j]][3]) - as.numeric(prior_lhs[[j]][2])) *
random_tab[lhs_index, j]
  }
  param <- rbind(param, temp_par)
  lhs_index <- lhs_index+1
}
param <- param[-1,]

Best regards,
Shah

	[[alternative HTML version deleted]]