I have the basic modechoice model estimated from the apollo_estimate()
model = apollo_estimate(apollo_beta, apollo_fixed, apollo_probabilities, apollo_inp