本文介绍了如何在利润率后评估结果后将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计入表格的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持!

10-24 15:24