I have unbalanced panel data on household cooking energy, with total observations equal to 1762 as follows:
Year
Count 2010 2013 2016 2019
1 62 12 17 422
2 57 55 20 8
3 307 308 302 7
4 46 46 46 46
Total 472 421 385 483
I would like to use the multinomial Probit to analyze the choice determinants in STATA 16 (I have 4 alternatives). I set my time variable, Year. First I need to know if I can still go ahead and use this unbalanced data, or if I should use only if the count==3 as it is the largest available. I thought of also analyzing using the count=4 data (46 households, 184 observations), just to be able to use the most recent data, but I am not sure if this will be enough for estimations.

Also, I haven't seemed to find a command for the multinomial Probit on panel data for STATA 16. Please assist with this as well.
Thank you