I did a scenario-analysis with a simple ordered logit modle using apollo_prediction():
Code: Select all
pred <- apollo_prediction(Model.3,apollo_probabilities,apollo_inputs,prediction_settings = list(runs=100))
The issue is that I'm applying a bootstrap-like procedure like so:
Code: Select all
MS.save <- list()
for(j in 1:500){
wt.new <- manipulate_wt(wt.original,new.wt.distribution) # randomly assigns WaitingTime intervals to the observations according to new.wt.distribution using sample()
database$WaitingTime <- wt.new
apollo_inputs <- apollo_validateInputs()
pred <- apollo_prediction(Model.3,apollo_probabilities,apollo_inputs,prediction_settings = list(runs=100))
MS.save[[j]] <- colMeans(pred$at_estimates[,3:5]) # save market shares for each j
}
Best,
Lucas