同一の目的変数(従属変数)に対し、
複数の説明変数(独立変数)で単変量ロジスティック回帰分析を繰り返し、
odds ratio, p value, [95% CI]
をこの順で、Excelに貼り付けしやすい形で取り出す
ためのコード。
表示桁数も手ごろそうな感じに揃えてある。
多変量をやる前の準備的な感じ。
下記をdo fileに貼り付けて実行する。
→STATAの結果windowに出たものを、
選択・右クリック→表のコピー してExcelに貼り付け。
qui{
noi di "(Category.)説明変数", "OR", "Pvalue", "95%CI"
foreach v in 説明変数1 説明変数2 (説明変数3...){
/* カテゴリ変数もOK, 変数名の前にi.を付ける。*/
logistic 目的変数 `v'
matrix T=r(table)
local ColN = `=colsof(T)'
forvalues n = 1/`ColN'{
local ColName`n' = "`: word `n' of `: colnames T''"
noi di "`ColName`n''", %4.2f T["b",`n'], %4.3f T["pvalue", `n'], "[" %4.2f T["ll", `n'] "–" %4.2f T["ul", `n'] "]"
}
}
}