Экспортируйте p-значения из reghdfe в Excel - PullRequest
Новая
       13

Экспортируйте p-значения из reghdfe в Excel

1 голос
/ 22 марта 2019

Я пытаюсь экспортировать p-значения в Excel из команды Stata , предоставленной сообществом reghdfe: 

// get data
use "http://www.stata-press.com/data/r9/nlswork.dta", clear
xtset idcode year

// run regression
reghdfe ttl_exp age not_smsa msp nev_mar, abs(idcode year) cluster(idcode)

// store numbers
local rmse = `e(rmse)'  

// export
putexcel A1 = (`rmse') using "export.xlsx" , modify // this works

Как я могу получить p-значение заданной переменной (скажем, msp), экспортированной в тот же документ Excel?

ereturn list не предлагает никакой полезной информации, e(V) содержит только ковариационную матрицу.

Крест опубликовал в Statalist .

1 Ответ

2 голосов
/ 22 марта 2019

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

use "http://www.stata-press.com/data/r9/nlswork.dta", clear
xtset idcode year

reghdfe ttl_exp age not_smsa msp nev_mar, abs(idcode year) cluster(idcode)

------------------------------------------------------------------------------
             |               Robust
     ttl_exp |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         age |   -.008044   .0618821    -0.13   0.897     -.129366    .1132781
    not_smsa |  -.1616256   .0812134    -1.99   0.047    -.3208472   -.0024039
         msp |   -.004893   .0634814    -0.08   0.939    -.1293505    .1195645
     nev_mar |    -.36351   .0955207    -3.81   0.000    -.5507816   -.1762384
------------------------------------------------------------------------------

local pval = (2 * ttail(e(df_r), abs(_b[msp] / _se[msp]) ) )

display `pval'
.93856498

putexcel A1 = (`pval') using "export.xlsx", modify
...