@junpenglao Thank you very much! Your code works for the entire data set which contains 100 households – 100 3-component alpha distributions, each alpha representing an individual household. I also added 2 fixed coefficients (using classical statistics terminology) – coefficient distribution with 2 components - a mixed logit model. Works well. Thanks again!
2 Likes