ʵ֤Ñо¿ÖÐÓ¦Óù㷺£¬Ê¹ÈËÃÇÄܸü¼Ó׼ȷµØ°ÑÎÕ·çÏÕ£¨²¨¶¯ÐÔ£©£¬ÓÈÆäÊÇÓ¦ÓÃÔÚ·çÏÕ¼ÛÖµ£¨VALUE AT RISK£©ÀíÂÛÖУ¬ÔÚ»ª¶û½ÖÊÇÈ˾¡½ÔÖªµÄ¹¤¾ß¡£
ËùνARCHÄ£ÐÍ£¬°´ÕÕÓ¢ÎÄÖ±ÒëÊÇ×ԻعéÌõ¼þÒì·½²îÄ£ÐÍ¡£´ÖÂÔµØËµ£¬¸ÃÄ£Ðͽ«µ±Ç°Ò»ÇпÉÀûÓÃÐÅÏ¢×÷ΪÌõ¼þ£¬²¢²ÉÓÃijÖÖ×ԻعéÐÎʽÀ´¿Ì»®·½²îµÄ±äÒ죬¶ÔÓÚÒ»¸öʱ¼äÐòÁжøÑÔ£¬ÔÚ²»Í¬Ê±¿Ì¿ÉÀûÓõÄÐÅÏ¢²»Í¬£¬¶øÏàÓ¦µÄÌõ¼þ·½²îÒ²²»Í¬£¬ÀûÓÃARCH Ä£ÐÍ£¬¿ÉÒÔ¿Ì»®³öËæÊ±¼ä¶ø±äÒìµÄÌõ¼þ·½²î¡£
ARCH(m)Ä£ÐÍ£º
yt?xt???t?t2??0??1?t2?1??2?t2?2????m?t2?mÆäÖУ¬?2ÊÇ²Ð²îÆ½·½ºÍ£¨²¨¶¯ÂÊ£© ?iÊÇARCHÄ£Ð͵ÄϵÊý GARCH(m,k)Ä£ÐÍ£º
£¨Ìõ¼þƽ¾ùÖµ£©
£¨Ìõ¼þ·½²î£©yt?xt???t???0???2t21t?1???22t?2?????2mt?m??1?2t?1??2?2t?2????k?2t?k
ÆäÖУ¬?iÊÇARCHÄ£Ð͵ÄϵÊý£»?iÊÇGARCHϵÊý 3.1 ARCHÄ£ÐÍÓ¦Óà Àý×Ó£º
. use http://www.stata-press.com/data/r11/wpi1,clear . regress D.ln_wpi
Source | SS df MS Number of obs = 123 -------------+------------------------------ F( 0, 122) = 0.00 Model | 0 0 . Prob > F = . Residual | .02521709 122 .000206697 R-squared = 0.0000 -------------+------------------------------ Adj R-squared = 0.0000 Total | .02521709 122 .000206697 Root MSE = .01438
------------------------------------------------------------------------------ D.ln_wpi | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- _cons | .0108215 .0012963 8.35 0.000 .0082553 .0133878 ------------------------------------------------------------------------------
. estat archlm,lags(1)
LM test for autoregressive conditional heteroskedasticity (ARCH)
--------------------------------------------------------------------------- lags(p) | chi2 df Prob > chi2 -------------+------------------------------------------------------------- 1 | 8.366 1 0.0038
--------------------------------------------------------------------------- H0: no ARCH effects vs. H1: ARCH(p) disturbance
ͨ¹ý¶ÔWPIµÄ¶ÔÊý²î·Ö½øÐг£Êý»Ø¹é£¬½Ó×ÅÓÃLM¼ìÑéÀ´ÅжÏARCH(1)ЧӦ,ÔÚ¸ÃÀý×ÓÖУ¬¼ìÑéµÄ½á¹ûPROB > CHI2£½0.0038<0.05£¬ËùÒԾܾøÃ»ÓÐARCH(1)ЧӦµÄÐéÎÞ¼ÙÉè¡£Òò´Ë£¬ÎÒÃÇ¿ÉÒÔͨ¹ýÖ¸¶¨ARCH(1)Ä£ÐÍÀ´¹À¼ÆARCH(1)µÄϵÊý¡£
. arch D.ln_wpi,arch(1) garch(1) ARCH family regression
Sample: 1960q2 - 1990q4 Number of obs = 123 Distribution: Gaussian Wald chi2(.) = . Log likelihood = 373.234 Prob > chi2 = .
------------------------------------------------------------------------------ | OPG
D.ln_wpi | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- ln_wpi |
_cons | .0061167 .0010616 5.76 0.000 .0040361 .0081974 -------------+---------------------------------------------------------------- ARCH | arch |
L1. | .4364123 .2437428 1.79 0.073 -.0413147 .9141394 | garch |
L1. | .4544606 .1866605 2.43 0.015 .0886126 .8203085 |
_cons | .0000269 .0000122 2.20 0.028 2.97e-06 .0000508 ------------------------------------------------------------------------------
ÕâÑù£¬ÎÒÃǾͿÉÒÔ¹À¼Æ³öÁËARCH(1)µÄϵÊýÊÇ0.436,GARCH(1)µÄϵÊýÊÇ0.454,ËùÒÔÎÒÃÇ¿ÉÒÔÄâºÏ³öGARCH(1,1)Ä£ÐÍ£º
yt?0.0061??t??0.436?2t2t?1?0.454?2t?1 ÆäÖУ¬yt?ln(wpit)?ln(wpit?1)
½ÓÏÂÀ´ÎÒÃÇ¿ÉÒÔ¶Ô±äÁ¿µÄ½øÐÐÔ¤²â£º
predict xb,xb /*¶Ô²î·Ö±äÁ¿µÄÔ¤²â*/ predict y,y /*¶Ôδ²î·Ö±äÁ¿µÄÔ¤²â*/
predict variance,var /*¶ÔÌõ¼þ·½²îµÄÔ¤²â */
predict res,residuals /*¶Ô²î·Ö±äÁ¿²Ð²îµÄÔ¤²â*/
predict yres,yresiduals /*¶Ôδ²î·Ö±äÁ¿²Ð²îµÄÔ¤²â*/
3.2 ARCHÄ£Ð͵ÄÈ·¶¨ÒÔ¼°¼ìÑé
Àý×Ó£º
use http://www.stata-press.com/data/r11/wpi1,clear
*- ¼ìÑé ARCH ЧӦÊÇ·ñ´æÔÚ£ºarchlm ÃüÁî regress D.ln_wpi archlm, lag(1/20)
regress D.ln_wpi L(1/3).D.ln_wpi archlm, lag(1/20)
* ͼÐη¨¡ª¡ª×ÔÏà¹Øº¯Êýͼ (ac) reg D.ln_wpi predict e, res gen e2 = e^2 ac e2, lag(40)
gen dlnwpi=D.ln_wpi gen dlnwpi2 = dlnwpi^2 ac dlnwpi2, lag(40) * ¾«¼òÄ£ÐÍ£ºARCH(1) * ±£ÊØÄ£ÐÍ£ºARCH(4)
*- Ô¤²âÖµ
arch D.ln_wpi, arch(1/4)
predict ht, variance /*Ìõ¼þ·½²î*/
* ht = c + a_1*e2_t-1 + a_2*e2_t-2 + ... + a_5*e2_t-5 line ht t
predict et, residual /*¾ùÖµ·½³ÌµÄ²Ð²î*/
*- Ä£ÐÍµÄÆÀ¹À * »ù±¾Ë¼Ï룺
* ÈôÄ£ÐÍÉ趨ÊǺÏÊʵģ¬ÄÇô±ê×¼»¯²Ð²î * z_t = e_t/sqrt(h_t)
* ӦΪһ¸ö i.i.d µÄËæ»úÐòÁУ¬¼´²»´æÔÚÐòÁÐÏà¹ØºÍARCHЧӦ£»
gen zt = et / sqrt(ht) /*±ê×¼»¯²Ð²î*/
gen zt2 = zt^2 /*±ê×¼»¯²Ð²îµÄƽ·½*/
* ÐòÁÐÏà¹Ø¼ìÑé pac zt
corrgram zt /*Ljung-Box ͳ¼ÆÁ¿*/
pac zt2
corrgram zt2
* Õý̬·Ö²¼¼ìÑé
histogram zt, normal wntestb zt wntestb zt2
* ÆÀÂÛ£º¾ùÖµ·½³ÌµÄÉ趨¿ÉÄÜÐèÒª¸Ä½ø£¬ÒòΪ zt ÈÔÈ»±íÏÖ³öÃ÷ÏÔµÄÐòÁÐÏà¹Ø¡£ * Ìõ¼þ·½²î·½³ÌµÄÉ趨»ù±¾Âú×ãÒªÇó£¬zt2 ²»´æÔÚÃ÷ÏÔµÄÐòÁÐÏà¹Ø¡£
3.3 ARIMA¹ý³ÌµÄARCHÄ£ÐÍ
ÎÒÃÇ¿ÉÒÔ¶ÔÌõ¼þ·½²îÄ£Ðͱ£³ÖARCH(1,1)Ä£ÐͶø¾ùֵģÐͲÉÓÃARMA¹ý³ÌµÄ×ԻعéÒ»½×ºÍÒÆ¶¯Æ½¾ùÒ»½×Å©ÒÔ¼°Òƶ¯Æ½¾ùËĽ×À´¿ØÖƼ¾½ÚÓ°Ï죺
. use http://www.stata-press.com/data/r11/wpi1,clear . arch D.ln_wpi,ar(1) ma(1 4) arch(1) garch(1)
ARCH family regression -- ARMA disturbances
Sample: 1960q2 - 1990q4 Number of obs = 123 Distribution: Gaussian Wald chi2(3) = 153.56 Log likelihood = 399.5144 Prob > chi2 = 0.0000
------------------------------------------------------------------------------ | OPG
D.ln_wpi | Coef. Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- ln_wpi |
_cons | .0069541 .0039517 1.76 0.078 -.000791 .0146992 -------------+---------------------------------------------------------------- ARMA | ar |
L1. | .7922673 .1072225 7.39 0.000 .582115 1.002419 | ma |
L1. | -.3417738 .1499944 -2.28 0.023 -.6357574 -.0477902 L4. | .2451725 .1251131 1.96 0.050 -.0000446 .4903896 -------------+---------------------------------------------------------------- ARCH | arch |
L1. | .2040451 .1244992 1.64 0.101 -.039969 .4480591 | garch |
L1. | .694968 .189218 3.67 0.000 .3241075 1.065829 |