我正在尝试使用R估计以下具有 3 个内生变量的联立方程组:
y~a+B
B~c+D
D~f+g
实际上,我正在尝试重复来自Stata的结果。Stata中使用3SLS方法的代码是
reg3 (y=a B) (B=c D) (D=f g)
在 RI 中使用systemfit包。
library("systemfit")
systemfit(data, list(y~a+B, B~c+D,D~f+g), method = "3SLS")
但是这个函数是用于 3SLS 估计的工具变量。
我做错了什么?提前致谢。