列名でStata行列要素を抽出する

Question

回帰を実行した後にvarlistにループするように、Stata行列要素を列名で読み取ってみたいと思います。

次のコードでは時々動作しますが、すべてではありません。そして私は理由を理解できません。

sysuse auto, clear 
scalar drop _all 
matrix drop _all 

regress price weight headroom trunk 
matrix b = e(b) 
local colnames : colnames b 
foreach v of local colnames { 
    if ("`v'" != "_cons") { 
     display "-------------------------" 
     display "`v'" 
     quietly summarize `v', detail 
     scalar beta = b[1,`v'] 
     scalar p75 = `r(p75)' 
     scalar p25 = `r(p25)' 
     scalar iqr = p75 - p25 
     scalar effect = beta * iqr 
     scalar dir 
    } 
} 

それはweightとtrunkではなくheadroomで失敗しました。 weightとtrunkはどちらも負ですが、私自身の別のコードでは正の係数では失敗します。

マニュアルカウンターを使用すると（下記）、なぜ列名で失敗するのですか？別の解決策がありますか？

regress price weight headroom trunk 
matrix b = e(b) 
local colnames : colnames b 
local i = 1 
foreach v of local colnames { 
    if ("`v'" != "_cons") { 
     display "-------------------------" 
     display "`v'" 
     quietly summarize `v', detail 
     scalar beta = b[1,`i'] 
     local ++i 
     scalar p75 = `r(p75)' 
     scalar p25 = `r(p25)' 
     scalar iqr = p75 - p25 
     scalar effect = beta * iqr 
     scalar dir 
    } 
} 

Answer 1

は、列名を使用して添字であなたの試みはうまくいかなかった理由の説明のためのhelp matrix subscriptingを参照してください。行を置換すると

 scalar beta = b[1,colnumb(b,"`v'")] 

が動作します。