create a random model for EM with three sided endogenous mobility with multinomial pr

m2.mixt.new(nk, nf, nb, fixb = F, stationary = F)