set.seed(101010101)
utility <- data.frame(female = rbinom(size = 1,
n = 200,
p = 0.6))
utility[, 'age'] <- stats::rnorm(n = 200,
mean = 50 + utility$female*10,
sd = 15)
utility[1:50, 'eq5d'] <- stats::rnorm(n = 50,
mean = 0 - 0.1 *
utility[1:50, 'female'] +
0.0005 * utility[1:50, 'age'],
sd = 0.1)
utility[51:200, 'eq5d'] <- stats::rnorm(n = 150,
mean = 0.5 +
0.1 * utility[51:200, 'female'] +
0.0001*utility[51:200, 'age'],
sd = 0.2)
utility[utility$eq5d<(-0.594), 'eq5d'] <- -0.594
utility[utility$eq5d>0.883, 'eq5d'] <- 1
hist(utility$eq5d, breaks = 50)
Run the code above in your browser using DataLab