大致解决思路,还可以优化
下面:var1代表age var2代表incom
sort var1
gen s1=_n
sum s1
local Ns1=r(max)
forvalues t=1/`Ns1'{
gen m`t'=var1+5 if s1==`t'
gen n`t'=var1-5 if s1==`t'
quiet sum m`t'
gen k`t'=r(max)
quiet sum n`t'
gen l`t'=r(max)
quiet sum var2 if k`t'>=var1&var1>=l`t'
gen h`t'=r(mean)
}