问题是这样，有面板模型Yit=ri+b1Xit+b2Zit+C1(Xit* Yit)+C2(Zit* Yit)+eit，经过差分消除个体效应后，差分模型同时存在随机扰动项相关和内生性问题，理论上应该使用工具变量估计。

      请问，如果使用stata估计，这类模型该使用什么命令？