[Stata]单因素资本资产定价模型CAPM在Stata中的实现

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

[Stata]单因素资本资产定价模型CAPM在Stata中的实现

// 单因素资本资产定价模型在stata中的实现 // 数据来自25个资产池账面市值、无风险利率以及SP&500指数，数据来自Ken French主页，详见附件 // http://mba.tuck.dartmouth.edu/pages/faculty/ken.french/data_library.html *version 13.1 clear all // 设置日志文件 capture log close qui log using capm,replace // 添加数据标签 use fama_french.dta, clear label variable rm "Market Return" label variable rf "Risk Free Rate" // 时间序列设定 tsset dateid01 // 描绘市场预期收益率和无风险利率 twoway (tsline rm) (tsline rf), name(factors0, replace) ///             xlabel(192607(000800)201312,angle(forty_five)) xtitle("") graph export "../factors0.pdf", as(pdf) replace //将图片保存为pdf格式 // plot return on market and rf on same graph using different y-axes twoway (tsline rm, yaxis(1)) (tsline rf, yaxis(2)), name(factors, replace) ///             xlabel(192607(000800)201312,angle(forty_five)) xtitle("")   graph export "../factors.pdf", as(pdf) replace                          drop if _n>1021 // 产生超额收益 local N     = 25 forvalues i = 1/`N' {   qui gen z`i' = r`i' - rf } // drop r1-r9 // 对第一个基础资产池回归 reg z1 rm_rf // 如果alpha=0,表示市面价格（book-market-value）准确反映其内在价值，未被高估也未被低估。 //如果beta=1,表示该组合收益风险与市场波动状况一致 test _cons    test rm_rf=1 test (_cons=0) (rm_rf=1) // 另外4种检验方法 local N=25 forvalues i = 1/`N' {   qui regress z`i' rm_rf } //使用mvreg OLS mvreg z* = rm_rf test _cons //使用似不相关回归 sureg z* = rm_rf test _cons // 使用statby命令,结果保存在simplecapm.dta中 reshape long z, i(dateid01) j(portfolio) statsby  _b _se, by(portfolio) saving(simplecapm, replace): reg z rm_rf qui log close 复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏40 回帖

关键词：资本资产定价模型 资本资产定价 Stata tata 资产定价 模型

好贴，谢谢分享。

多谢分享，收益匪浅

主贴内容参考自：

                      Financial Econometric Modelling (Using Stata)

——Professor Stan Hurn, QUT and NCER Director (Singapore Management University)

1. //讲义下载
   
2. copy "http://www.ncer.edu.au/events/documents/Singapore-lecture3.pdf" "Singapore-lecture3.pdf"
   
3. 
   
4. //do-file下载
   
5. copy "http://www.ncer.edu.au/events/documents/Singapore-lecture3.zip" "Singapore-lecture3.zip"
   
6. 
   
7. //sjlog 生成符合stata journal 排版格式的日志文件
   
8. net inst sjlatex, from(http://www.stata-journal.com/production)
   
9. 
   
10. //建立一个名为working文件夹与后续范例匹配
    
11. cap findfile working
    
12. if "`r(filename)'" == "" mkdir working
    
13. 
    
14. // Testing one factor pricing models using data from Ken French's website
    
15. // Looking at the data and the simple CAPM
    
16. 
    
17. 
    
18. sjlog using "./working/lec3_aa", replace
    
19. set more off
    
20. version 13
    
21. clear *
    
22. 
    
23. //clear * 与clear all 同义
    
24. 
    
25. // set up a log file
    
26. capture log close capm
    
27. log using capm, name(capm) replace
    
28. 
    
29. // set current working directory
    
30. * cd ~/Dropbox/Teaching/Singapore/do 可自行设定，文中仅示例
    
31. 
    
32. // load the data and add a few labels 注意数据放在当前文件夹下(数据在主贴)
    
33. use fama_french.dta, clear
    
34. label variable rm "Market Return"
    
35. label variable rf "Risk Free Rate"
    
36. 
    
37. sjlog close, replace
    
38. 
    
39. sjlog using "./working/lec3_bb", replace
    
40. // format date variable and set data set as time series
    
41. 
    
42. drop dateid01
    
43. gen dateid01 = tm(1926m7)+_n-1
    
44. format %tm dateid01
    
45. tsset dateid01
    
46. // plot return on market and rf on same graph using same y-axes
    
47. twoway (tsline rm) (tsline rf), name(factors0, replace) ///
    
48.             tlabel(,angle(forty_five) format(%tmCCYY)) xtitle("")
    
49. 
    
50. graph export "../factors0.pdf", as(pdf) replace
    
51. 
    
52. // plot return on market and rf on same graph using different y-axes
    
53. twoway (tsline rm, yaxis(1)) (tsline rf, yaxis(2)), name(factors, replace) ///
    
54.             tlabel(,angle(forty_five) format(%tmCCYY)) xtitle("")
    
55. 
    
56. graph export "../factors.pdf", as(pdf) replace
    
57. sjlog close, replace
    
58. 
    
59. 
    
60. // generate excess returns
    
61. local N = 25
    
62. forvalues i = 1/`N' {
    
63.   qui gen z`i' = r`i' - rf
    
64. }
    
65. 
    
66. // drop raw returns
    
67. drop r1-r9 // note use of hyphen r10 comes right after r1
    
68. 
    
69. sjlog using "./working/lec3_cc", replace
    
70. // estimate CAPM for first portfolio and test alpha = 0 and beta = 1
    
71. reg z1 rm_rf
    
72. 
    
73. // test the model
    
74. test _cons
    
75. test rm_rf=1
    
76. test (_cons=0) (rm_rf=1)
    
77. sjlog close, replace
    
78. 
    
79. reg z2 rm_rf
    
80. reg z3 rm_rf
    
81. 
    
82. // at least four ways to do this estimation
    
83. 
    
84. forvalues i = 1/`N'{
    
85.   qui regress z`i' rm_rf
    
86. }
    
87. 
    
88. mvreg z* = rm_rf
    
89. 
    
90. sureg z* = rm_rf
    
91. 
    
92. sjlog using "./working/lec3_dd", replace
    
93. reshape long z, i(dateid01) j(portfolio)
    
94. sjlog close, replace
    
95. 
    
96. statsby  _b _se, by(portfolio) saving(simplecapm, replace): reg z rm_rf
    
97. 
    
98. log close capm
    
99. 
    
100. // gen daten = tm(1926m7)+_n-1

复制代码

---

---

很好，很有用，谢谢楼主

niuniuyiwan 发表于 2016-1-4 20:52
主贴内容参考自：
                      Financial Econometric Modelling (Using Stata) ——Professor  ...

3楼内容应补：

1. set scheme s1color

复制代码

非常好~~有帮助

群主 你好 我找到了你给出网址的数据，但是我发现，你的代码针对的数据是你整理过后的数据，你介意把你整理过后的数据发给我吗？？想跟着你的代码以及数据过这个过一遍。
来自一个自学stata还没看懂你的代码的学渣~·

小小小漫 发表于 2017-4-25 16:35
群主 你好 我找到了你给出网址的数据，但是我发现，你的代码针对的数据是你整理过后的数据，你介意把你整理 ...

已经找到了  谢谢楼主

小小小漫 发表于 2017-4-27 21:46
已经找到了  谢谢楼主

你好，请问能分享下你拿到的数据么？楼主贴的数据似乎不对。