试试(红字部分请根据你的具体需要进行修改,希望可以满足你的要求):
// 示例数据
webuse klein, clear
// main
capture program drop bs_coef
program bs_coef, rclass
reg3 (consump wagepriv wagegovt) (wagepriv consump govt capital1)
nlcom [consump]_b[wagepriv]^2 + [wagepriv]_b[govt]^0.5 // 非线性组合
mat result = r(b) // nlcom's Coef, 目标指标
return scalar coef = result[1,1]
/* 计算95%置信区间上下限,但好像用不着吧
mat V_coef = r(V)
scalar se_coef = V_coef[1,1]^0.5 // Std. Err.
return scalar u_coef = result[1,1] + invnormal(0.975)* se_coef
return scalar l_coef = result[1,1] - invnormal(0.975)* se_coef
*/
end
bootstrap zhibiao=r(coef), rep(100) saving(jieguo, replace): bs_coef
use jieguo, clear
br