estat icによって計算されたStataのBIC値にアクセスする方法は？

Question

This documentは、AICとBICの値がr(S)に格納されていると説明していますが、display r(S)を試してみると「型の不一致」が返され、sum r(S)を試してみると "r ambiguous abbreviation"が返されます。

私はこの誤解をおかけして申し訳ありません。r(S)ですが、計算されたBIC値にどのようにアクセスできるかご連絡いただければ幸いです。

Answer 1

参照先のドキュメントには、r(S)がマトリックスであることが記載されています。 displayコマンドは行列では機能しません。試してくださいmatrix list。 help matrixも参照してください。例えば

：

clear 
sysuse auto 
regress mpg weight foreign 
estat ic 
matrix list r(S) 
matrix S=r(S) 
scalar aic=S[1,5] 
di aic 

Answer 2

あなたが引用した同じ文書がr(S)が行列であると説明しています。これは、あなたのコマンドの失敗を説明します。は、変数を要約するためのものです。displayは、ヘルプが説明するように、文字列とスカラー式を表示するためのものです。行列もどちらでもない。アドバイスは同じいずれかの方法ですが、ドキュメントはあなたが

http://www.stata.com/manuals13/restatic.pdf

を引用

ノートには、最新バージョン

http://www.stata.com/manuals14/restatic.pdf

を書いていない時です。

r(S)を、次のr-classコマンドを実行したときに消えることのない行列にコピーし、直接一覧表示します。行列の基本的なヘルプについては、

help matrix 

ここは再現可能な例です。あなたの質問は、あなたがそのバージョンを使用している可能性があることをほのめかすので、私は、データセットのStataの13バージョンを使用します。

. use http://www.stata-press.com/data/r13/sysdsn1 
(Health insurance data) 

. mlogit insure age male nonwhite 

Iteration 0: log likelihood = -555.85446 
Iteration 1: log likelihood = -545.60089 
Iteration 2: log likelihood = -545.58328 
Iteration 3: log likelihood = -545.58328 

Multinomial logistic regression     Number of obs  =  615 
               LR chi2(6)  =  20.54 
               Prob > chi2  =  0.0022 
Log likelihood = -545.58328      Pseudo R2   =  0.0185 

------------------------------------------------------------------------------ 
     insure |  Coef. Std. Err.  z P>|z|  [95% Conf. Interval] 
-------------+---------------------------------------------------------------- 
Indemnity | (base outcome) 
-------------+---------------------------------------------------------------- 
Prepaid  | 
     age | -.0111915 .0060915 -1.84 0.066 -.0231305 .0007475 
     male | .5739825 .2005221  2.86 0.004  .1809665 .9669985 
    nonwhite | .7312659 .218978  3.34 0.001  .302077 1.160455 
     _cons | .1567003 .2828509  0.55 0.580 -.3976773 .7110778 
-------------+---------------------------------------------------------------- 
Uninsure  | 
     age | -.0058414 .0114114 -0.51 0.609 -.0282073 .0165245 
     male | .5102237 .3639793  1.40 0.161 -.2031626  1.22361 
    nonwhite | .4333141 .4106255  1.06 0.291  -.371497 1.238125 
     _cons | -1.811165 .5348606 -3.39 0.001 -2.859473 -.7628578 
------------------------------------------------------------------------------ 

. estat ic 

Akaike's information criterion and Bayesian information criterion 

----------------------------------------------------------------------------- 
     Model |  Obs ll(null) ll(model)  df   AIC  BIC 
-------------+--------------------------------------------------------------- 
      . |  615 -555.8545 -545.5833  8 1107.167 1142.54 
----------------------------------------------------------------------------- 
       Note: N=Obs used in calculating BIC; see [R] BIC note. 

. ret li 

matrices: 
        r(S) : 1 x 6 

. mat S = r(S) 

. mat li S 

S[1,6] 
      N   ll0   ll   df   AIC   BIC 
.   615 -555.85446 -545.58328   8 1107.1666 1142.5395 

BIC値がS[1,6]になりました。