Prob = ordered_logit(
continuous_value=V,
list_of_discrete_values=[1, 2, 3],
tau_parameter=tau1_V,
)
Prob_chosen = Elem(Prob, Relax)
condlike = (Prob_chosen * P_Car1 * P_Car2 * P_Car3)
#Panel data
condprobIndiv = PanelLikelihoodTrajectory(condlike)
loglike = log(MonteCarlo(condprobIndiv))