（可能几乎最快之）股票崩盘危机指标之指令 (I)

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

大家最常用之 statsby

1. set rmsg on
   
2. 
   
3. cap log close
   
4. log using "cr_statsby.log", replace
   
5. 
   
6. clear
   
7. set seed 123
   
8. 
   
9. // 2000 firms
   
10. set obs 2000
    
11. gen long id = _n
    
12. 
    
13. // 18 years
    
14. expand 18
    
15. bys id: gen year = 1999 + _n
    
16. 
    
17. // 250 daily / 50 weekly observations
    
18. expand 50
    
19. 
    
20. gen ri = runiform()
    
21. gen rm = runiform()
    
22. 
    
23. bys id (year): gen t = _n
    
24. xtset id t
    
25. 
    
26. gen L1rm = L1.rm
    
27. gen L2rm = L2.rm
    
28. gen F1rm = F1.rm
    
29. gen F2rm = F2.rm
    
30. 
    
31. // statsby
    
32. statsby _b r2 = e(r2), by(id year) saving("r2.dta", replace): reg ri L2rm L1rm rm F1rm F2rm
    
33. merge m:1 id year using "r2.dta"
    
34. gen e_statsby = ri - (_b_L2rm*L2rm + _b_L1rm*L1rm + _b_rm*rm + _b_F1rm*F1rm + _b_F2rm*F2rm + _b_cons)
    
35. 
    
36. // firm-specific returns
    
37. gen W = ln(1+e_statsby)
    
38. 
    
39. // de-mean
    
40. bys id year: egen W_mean = mean(W)
    
41. replace W = W - W_mean
    
42. 
    
43. * CRASH & CRASH_count
    
44. bys id year: egen W_sd = sd(W)
    
45. gen W_norm = (W-W_mean)/W_sd
    
46. 
    
47. gen crash_week = 1 if W_norm <= -3.2
    
48. replace crash_week = 0 if crash_week ==.
    
49. 
    
50. bys id year: egen crash_count = total(crash_week)
    
51. gen crash = 1 if crash_count != 0
    
52. replace crash = 0 if crash_count == 0
    
53. 
    
54. * NCSKEW
    
55. gen W2 = W^2
    
56. gen W3 = W^3
    
57. bys id year: egen TW2 = total(W2)
    
58. bys id year: egen TW3 = total(W3)
    
59. bys id year: gen n = _N
    
60. gen ncskew = -(TW3*n*(n-1)^1.5)/((n-1)*(n-2)*(TW2)^1.5)
    
61. 
    
62. * DUVOL
    
63. gen down = 0
    
64. replace down = 1 if W < 0
    
65. bys id year: egen down_sd = sd(W) if down == 1
    
66. bys id year: egen up_sd = sd(W) if down == 0
    
67. 
    
68. * ALL
    
69. collapse ncskew down_sd up_sd, by(id year)
    
70. gen duvol = ln(down_sd/up_sd)
    
71. drop down_sd up_sd
    
72. 
    
73. save "cr_statsby.dta", replace
    
74. 
    
75. sum
    
76. 
    
77. log close

复制代码

扫码加我 拉你入群

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

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

分享0 收藏30 回帖

关键词：observations observation observat Specific statsby

谢谢黄老师分享

谢谢你

你好 黄老师 我刚刚STATA入门 所以很多地方不懂 想请问一下 您的代码为什么没有数据也可以跑出结果呢？ 但是您最后跑出来的两个指标跟现在大家一般均值在-0.2到-0.4之间有所出入 想知道什么原因呢 我如何用这套代码跑自己的数据呢？

woodorange 发表于 2018-11-25 19:09
你好 黄老师 我刚刚STATA入门 所以很多地方不懂 想请问一下 您的代码为什么没有数据也可以跑出结果呢？ 但是 ...

因为我自己产生"假的"资料。

您好，黄老师，* NCSKEW gen W2 = W^2 gen W3 = W^3，这边的W已经变成了replace W = W - W_mean，但是按照公式不是应该是W = ln(1+e_statsby)这个最初的W吗？谢谢

lclzh 发表于 2018-12-9 11:51
您好，黄老师，* NCSKEW gen W2 = W^2 gen W3 = W^3，这边的W已经变成了replace W = W - W_mean，但是按照公 ...

你讲得没错，只是我画蛇添足，按照一般统计概念先 de-mean，不过结果应该不会差太太多。

黃河泉 发表于 2018-12-9 14:51
你讲得没错，只是我画蛇添足，按照一般统计概念先 de-mean，不过结果应该不会差太太多。

那您这边DUVOL的计算跟公式也是有一定细微差别吧

lclzh 发表于 2018-12-9 16:53
那您这边DUVOL的计算跟公式也是有一定细微差别吧

你觉得困扰，就将 demean 的部分去掉！

最后出来的结果怎么没有股票代码了呀， 完全看不懂的我只想做出结果来