想请教一下，在双变量pribit（bivariate probit）模型引入一个工具变量的问题。

例如： biprobit (y1=x1 x2 x3) (y2=x1 x2 x3)
x1可能是内生性变量，x1与y1, x1与y2可能存在逆因果关系，现在想要引入工具变量z，请问应该如何写stata code呢？


是否有类似于ivprobit的命令，可以适用于双变量probit模型的情况呢？