我正在面板数据上运行systemfit,看似无关的具有25个方程的回归。这将返回一个systemfit对象,该对象是每个方程式的25个systemfit.equation对象的列表。这些systemfit.equation对象包含系数,残差,方程标签等。我试图将它们提取到一个25行的“整洁”数据框中,其中每个方程和列用于我感兴趣的每个方程的结果。>
我已经能够使用
提取变量,估计值,误差和p值results = as.data.frame(summary(fgls)$coefficients)%>%
+ rownames_to_column(var = "variable")
人造结果表如下所示,
+-----------------------+
| Variable Estimate |
+-----------------------+
| 1_(intercept) -21.2 |
| X1_X 6.2 |
| X1_Y 1.2 |
| 2_(intercept) -86.9 |
| X2_X 3.5 |
| X2_Y 2.7 |
| .... .... |
+-----------------------+
我想为此添加一个变量,以标识每个eqnlabel。 方程标签可以通过
提取fgls$eq[[i]]$eqnLabel)
对于方程i
我可以通过使用dplyr
重复变量标签来“机械地”创建变量,但是如何从systemfit对象中提取变量?