# NOT RUN {
# Simulation data
dt_strategies <- data.frame(strategy_id = c(1, 2, 3))
dt_patients <- data.frame(patient_id = seq(1, 3),
age = c(45, 50, 60),
female = c(0, 0, 1))
dt_states <- data.frame(state_id = seq(1, 3),
state_name = paste0("state", seq(1, 3)))
hesim_dat <- hesim_data(strategies = dt_strategies,
patients = dt_patients,
states = dt_states)
# Create StateVals object
fit_costs_medical <- stats::lm(costs ~ female + state_name, data = psm4_exdata$costs$medical)
dat <- expand(hesim_dat, by = c("strategies", "patients", "states"))
costs_medical <- create_StateVals(fit_costs_medical, input_data = dat, n = 5)
# Predict
head(costs_medical$sim(t = 5, type = "predict"))
# }
Run the code above in your browser using DataLab