Including the time-averages as controls for the individual effects is very intuitive: it allows to estimate the effect of changing Zit, while holding the time-average fixed. The model can be consistently estimated by pooled ordered probit. For valid inference, the standard errors should be made robust to arbitrary time dependence.22 As in the linear model, time-invariant variables should be dropped.