如何在利润率后评估结果后将AIC计入表格

问题描述

我有一个GLM模型，我在Stata中进行了估计。感兴趣的系数是边际效应，这是我用margins命令得到的。

然而，后估计表不包括像AIC这样的汇总统计数据，我希望在那里有这样的统计数据。

我已经通过编写辅助程序进行了尝试getAIC：

program getAIC
    estat ic
    matrix list r(S)
    matrix S = r(S)
    scalar aic = S[1,5]
end

然后按如下方式进行评估：

 qui glm y x, fa(bin) link(probit)
 getAIC
 qui margins, dydx(x) post
 estadd loc AIC aic

输出命令为：

 esttab using output.tex, s(aic, fmt(0))

但是，我的结果表中没有AIC。

有什么办法吗？

推荐答案

您需要从程序getAIC返回标量aic并相应地使用它。

以下内容适用于我：

program getAIC, rclass
    estat ic
    matrix list r(S)
    matrix S = r(S)
    scalar aic = S[1,5]
    return scalar aic = aic
end

sysuse auto, clear

glm foreign price, fa(bin) link(probit)

getAIC
local AIC = round(`r(aic)', .01)

margins, dydx(price) post

estadd local AIC `AIC'

esttab using output, s(AIC) replace

type output.txt
----------------------------
                      (1)   

----------------------------
price          0.00000766   
                   (0.43)   
----------------------------
AIC                 93.89   
----------------------------
t statistics in parentheses
* p<0.05, ** p<0.01, *** p<0.001

这篇关于如何在利润率后评估结果后将AIC计入表格的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！