Dear all,
I am working on household survey panel data (3 waves) to analyse the dynamic between 4 states in labour market. My data looks like:

Code:
* Example generated by -dataex-. To install: ssc install dataex
clear
input double(IDI year) float mobil_4 double age float(agesq female)
1 2012 4 50 2500 1
1 2013 4 51 2601 1
1 2015 1 53 2809 1
2 2012 1 20  400 1
2 2013 1 21  441 1
2 2015 1 23  529 1
3 2012 3 17  289 1
3 2013 1 18  324 1
3 2015 3 20  400 1
4 2012 4 10  100 1
4 2013 1 11  121 1
4 2015 1 13  169 1
My dependent variable is mobil_4 which contain 4 states in labour market. I tried to use -gllamm- to fit a dynamic multinomial logit model with random effects. The code used is adapted from the code in GLLAMM Manual, Rabe-Hesketh, Skondal, and Pickles (2004) and in Uhlendorff and Haan (2006):

Code:
global varx1 age agesq urban fondamental secondaire superieur form_prof 

mlogit     mobil_4 $varx1 , base(1)
matrix b=e(b)
matrix list b
scalar var=exp(0.5)
matrix start=b,var,var, var, 0.5, 0.5, 0.5
matrix colnames start= age agesq urban fondamental secondaire superieur form_prof _cons ///
age agesq urban fondamental secondaire superieur form_prof _cons ///
age agesq urban fondamental secondaire superieur form_prof _cons ///
infs other out _cons1 _cons2 _cons3
matrix coleq start= c2 c2 c2 c2 c2 c2 c2 c2 c2 ///
c3 c3 c3 c3 c3 c3 c3 c3 c3 ///
c4 c4 c4 c4 c4 c4 c4 c4 c4 ///
V2 V3 V4 COV12 COV13 COV23

sort IDI year age agesq urban fondamental secondaire superieur form_prof mobil_4
gen patt=_n
expand 4
sort patt
qui by patt:gen alt=_n
gen chosen=alt==mobil_4
sort patt alt
tab alt, gen(a)
gen dum=1
replace dum=0 if a1==1
eq dum:dum
eq a2:a2
eq a3:a3
eq a4:a4

gllamm alt $varx1 , expand(patt chosen m) i(IDI) link(mlogit) family(binom) nrf(3) /// 
eq(a2 a3 a4) nip(4) trace from(b) skip
The code didn't converge and took more than 12 hours in iteration 0, which means, I guess, that is something wrong in the code!
Code someone give me some insights on how to fix this issue?

I am using Stata 15.
Rferences:
Rabe-Hesketh, S., A. Skrondal, and A. Pickles. 2004.GLLAMM Manual. University of California–Berkeley, Division of Biostatistics, Working Paper Series. Paper No. 160.
http://www.bepress.com/ucbbiostat/paper160/.
Haan, P. and Uhlendorff, A. (2006), ‘Estimation of multinomial logit models with unobserverd heterogeneity using maximum simulated likelihood’, The Stata Journal 6(2), 229–245.