Я хотел бы запустить линейную регрессию с использованием категориальной переменной воздействия и вывести результаты на лист Excel с именами каждой ковариаты, включенными рядом с их результатами.
Код Stata ниже отлично работает для экспорта результатов:
sysuse auto.dta, clear
describe
summ
postfile temp str40 exp str40 outcome adjust N beta se lci uci pval using ///
"test.dta", replace
foreach out in price {
foreach exp in i.foreign {
foreach adjust in 1 2 {
if `adjust'==1 local adjusted ""
if `adjust'==2 local adjusted "mpg weight length displacement i.trunk"
reg `out' `exp' `adjusted'
local N = e(N)
matrix table=r(table)
forvalues i = 1 / `= colsof(r(table))-1' {
local beta = table[1,`i']
local se = table[2,`i']
local lci = table[5,`i']
local uci = table[6,`i']
local pval=table[4,`i']
post temp ("`exp'") ("`out'") (`adjust') (`N') (`beta') ///
(`se') (`lci') (`uci') (`pval')
}
}
}
}
postclose temp
use "test.dta", clear
Однако все строки помечены как i.foreign, поэтому трудно понять, какие результаты соответствуют другим ковариатам.
В идеале мне нужен столбец со строками, показывающими имена, соответствующие результатам, т.е. mpg, weight, length, displacement, i.trunk
