鉴于许多回归结果来自 xtable 不支持导出到 Latex 的类,我正在尝试创建自己的回归结果摘要。到目前为止,我已经能够创建以下矩阵:
> coef <- coefficients(summary(tob1))[,1]
> se <- sqrt(diag(vcov(tob1)))
> pvals <- 2*(1-pt(abs(coef)/se, df.residual(tob1)))
> ttob1 <- cbind(coef,se,pvals)
> ttob1
coef se pvals
X1 175.1252128 19.861699207 0.000000e+00
X2 -0.3377459 0.252531805 1.810945e-01
X3 -113.7189271 14.943601767 2.864375e-14
X4 50.8150461 13.863303952 2.475501e-04
X5 -244.2535872 12.739088876 0.000000e+00
X6 -244.4901871 13.447635109 0.000000e+00
X7 -3.5817039 1.675738225 3.257835e-02
现在我想创建一些函数来创建一个具有以下格式的表
X1 175.12
(19.86)
X2 -.3377
(.25)
X3 -113.72
(14.94)
我该怎么做?