我有个实证研究，需要用到联立方程来估计：
1. Y=f(Z1, Z2, x1, x2, x3...)
2. Z1=f(Y, Z2, x1, x2, x3...)
3. Z2=f(Y , Z1, x1, x2, x3...)
请问：在stata或SAS如何实现上述simultaneous equation的估计。
如果只有两个方程的话，在stata里用reg3可出来，但三个方程不知怎么实现。非常感谢！