選取云南白藥（000538）硫惕，基于深成指對云南白藥2015-2017年股票進(jìn)行估計
cntrade 000538
use "000538.dta",clear
global id="000538"
cntrade $id
global nn=stknme[1]
use "$id.dta",clear
tsset date
tsline clsprc,xlabel(,angle(60)) subtitle("$nn $id")
foreach i in "399001" "000001" "000300"{
cntrade `i',index
rename rmt retmkt_`i'
save index_`i'.dta,replace
erase `i'.dta
}
use "index_399001.dta",clear
merge 1:1 date using "$id.dta",nogen
gen year=year(date)
order date year index* stk* ret*
save "merge_$id.dta",replace
use "merge_$id.dta",clear
global yr_1=2014
global yr_2=2018
reg rit retmkt if year<$yr_2 & year>$yr_1? //對2015-2017年進(jìn)行回歸
aaplot rit retmkt if e(sample), xline(0,lp(dash) lc(red)) ///
? ? ? ? ? ? yline(0,lp(dash) lc(red)) msize(*0.6) ///?
? ? ? ? ? ? title("$nn ($id) 的 Beta 系數(shù), 2015-2017") ///
? ? ? ? ? ? xtitle("市場指數(shù)收益率(深成指)")
#graph export "估計貝塔.png" ?//做出各年貝塔值
bysort year: reg rit retmkt,noheader
statsby _b[retmkt], by(year) saving("beta_data.dta", replace): ///
? ? ? ? ? reg rit retmkt if year<$yr_2 & year>$yr_1? //對2015-2017年進(jìn)行逐年回歸
preserve
? ? use "beta_data.dta", clear
? ? keep if _stat!=.
? ? qui sum year
? ? local byr = r(min)
? ? local eyr = r(max)
? ? list
? ? #d ;
? ? twoway connect _stat year,?
? ? ? ? ? yline(1, lpattern(dash) lcolor(red*0.3))
? ? ? ? xlabel(`byr'(2)`eyr')
? ? ? ? ylabel(,angle(0) format(%2.1f))
? ? ? ? subtitle("$nn ($id) 的Beta系數(shù)")
? ? ? ? ytitle("beta") xtitle("");
? ? #d cr
? ? graph export "$out\beta $nn $id.png", replace
? restore
