model options maxthreads=all; algorithm tolerance=1e-008 emtolerance=0.01 emiterations=250 nriterations=50 ; startvalues seed=0 sets=16 tolerance=1e-005 iterations=50; bayes categorical=1 variances=1 latent=1 poisson=1; montecarlo seed=0 sets=0 replicates=500 tolerance=1e-008; quadrature nodes=10; missing excludeall; output parameters=first standarderrors estimatedvalues=model; variables caseid ID; dependent X, kT cumlogit; independent item nominal; latent Ability continuous, Speed continuous, Cluster nominal 2 dynamic; equations (1) Ability; Speed; Ability <-> Speed; Cluster <- 1; X <- 1 | Item Cluster + (+) Ability | Item Cluster; kT <- 1 | item + (-) Cluster + (aa) Speed; aa[1,1]=-1; end model