I’m running a multinomial regression model with random effects using gllamm and some doubts have come to me. I really would appreciate your help.

I have a panel data of 17785 observations of 2475 firms during a 14 years period with a polytomous dep_var with four categories. That would be a level 3 model, observation nested in firms nested in years, isn’t it?
When I try the next syntax, the level 1 and level 2 have the same number of observations!!

Moreover, I’m trying to include random intercepts for each alternative of the dependent variable at level 2 (firms) using eqs (a2 a3 a4) nrf (3 0) but this option is not allowed, since gllamm need to “split” the random effect in level 2 and 3. That is: