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

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

再来是用 (ssc install) rangestat:

1. set rmsg on
   
2. 
   
3. cap log close
   
4. log using "cr_rangestat.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. // rangestat
    
32. rangestat (reg) ri L2rm L1rm rm F1rm F2rm, interval(year 0 0) by(id)
    
33. gen e_rangestat = ri - (b_L2rm*L2rm + b_L1rm*L1rm + b_rm*rm + b_F1rm*F1rm + b_F2rm*F2rm + b_cons)
    
34. 
    
35. // firm-specific returns
    
36. gen W = ln(1+e_rangestat)
    
37. 
    
38. // de-mean
    
39. bys id year: egen W_mean = mean(W)
    
40. replace W = W - W_mean
    
41. 
    
42. * CRASH & CRASH_count
    
43. bys id year: egen W_sd = sd(W)
    
44. gen W_norm = (W-W_mean)/W_sd
    
45. 
    
46. gen crash_week = 1 if W_norm <= -3.2
    
47. replace crash_week = 0 if crash_week ==.
    
48. 
    
49. bys id year: egen crash_count = total(crash_week)
    
50. gen crash = 1 if crash_count != 0
    
51. replace crash = 0 if crash_count == 0
    
52. 
    
53. * NCSKEW
    
54. gen W2 = W^2
    
55. gen W3 = W^3
    
56. bys id year: egen TW2 = total(W2)
    
57. bys id year: egen TW3 = total(W3)
    
58. bys id year: gen n = _N
    
59. gen ncskew = -(TW3*n*(n-1)^1.5)/((n-1)*(n-2)*(TW2)^1.5)
    
60. 
    
61. * DUVOL
    
62. gen down = 0
    
63. replace down = 1 if W < 0
    
64. bys id year: egen down_sd = sd(W) if down == 1
    
65. bys id year: egen up_sd = sd(W) if down == 0
    
66. 
    
67. * ALL
    
68. collapse ncskew down_sd up_sd, by(id year)
    
69. gen duvol = ln(down_sd/up_sd)
    
70. drop down_sd up_sd
    
71. 
    
72. save "cr_rangestat.dta", replace
    
73. 
    
74. sum
    
75. 
    
76. log close

复制代码

扫码加我 拉你入群

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

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

关键词：observations observation observat Interval Specific

黄老师 想请教您一个问题
在计算股价崩盘风险收益上下波动比率DUVOL遇到一个问题 描述性统计计算出来其均值为正
但是看文献一般都是为-0.2至-0.4属于正常 这是为什么呢 看您之前回复的 说是可能是没有采用周资料 或者周顺序问题 但是这两个问题我都不存在 请问还会有其他什么原因会导致吗

小白不爱做白日梦 发表于 2019-4-15 10:03
黄老师 想请教您一个问题
在计算股价崩盘风险收益上下波动比率DUVOL遇到一个问题 描述性统计计算出来其均值 ...

我觉得是OK的！

黃河泉 发表于 2019-4-15 12:06
我觉得是OK的！

可是我看文献里面都是负值

黄老师 而且不知道为什么 相关性分析里面 会计稳健性与股价崩盘风险的两个指标均呈正相关并且显著

小白不爱做白日梦 发表于 2019-4-17 09:45
黄老师 而且不知道为什么 相关性分析里面 会计稳健性与股价崩盘风险的两个指标均呈正相关并且显著

无从评论！回归结果才重要！

黃河泉 发表于 2019-4-17 10:41
无从评论！回归结果才重要！

谢谢黄老师

小白不爱做白日梦 发表于 2019-4-15 10:03
黄老师 想请教您一个问题
在计算股价崩盘风险收益上下波动比率DUVOL遇到一个问题 描述性统计计算出来其均值 ...

你好 请问你解决了问题吗？我算出来的DUVOL也是均值为正值。