Hello,
I am performing a sensitivity analysis using a Latin Hypercube sampling.
However, I have difficulty to draw a Hypercube sample for one variable. I�ve
generated this variable from a Poisson distribution as follows:
set.seed(5)
mortality_probability <- round(ppois(seq(0, 7, by = 1), lambda = 0.9), 2)
barplot(mortality_probability, names.arg = seq(0, 7, by = 1), xlab = "Age
class", ylab = "Probability")
How can I draw a Hypercube sample for the variable �mortality_probability� so
that this variable exhibits the same pattern as the observed distribution?
Here is a reproducible code to draw Hypercube samples (my sensitivity analysis
includes several parameters and the variables �var1� and �var2� follow a
uniform distribution):
library(lhs)
set.seed(1)
parm <- c("var1", "var2", "mortality_probability")
X <- randomLHS(100, length(parm))
Any suggestions would be much welcome.
Thanks for your time
Marine
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
