用stata进行有序probit回归分析，该如何控制固定效应？

各位好，我在做10年内全国地级市市长环保绩效和晋升状况的相关度研究，并决定使用有序probit模型进行回归。但在查阅参考文献后，发现必须要控制年度和城市的固定效应。在查阅了stata的各种教程后，发现控制固定效应是单独的一个模型（通过xtreg指令实现），而stata里只有xtoprobit指令能实现控制随机效应的有序probit分析。那么现在我该如何再分析中控制以上固定效应呢？还望各位大神不吝赐教。

有同样的疑惑，请问楼主问题解决了吗？