Как включить AIC в таблицу после результатов оценки маржи - PullRequest
Новая
       18

Как включить AIC в таблицу после результатов оценки маржи

0 голосов
/ 26 июня 2018

У меня есть модель 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 в таблице результатов.

Есть идеи, как это сделать?

1 Ответ

0 голосов
/ 26 июня 2018

Вам необходимо вернуть скаляр aic из вашей программы getAIC и использовать его соответствующим образом.

У меня работает следующее: 

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
...