Current parameter estimation techniques assume order invariance. Investigate discrete state-space parameter recovery of the form

θ_n = A θ _n-1 + B μ
m = C θ_n + D μ

y = expit(m)

This is an experimental issue, but adding for posterity