[求]Multivariate EGARCH的code

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大虾好啊！
最近在做关于volatility transmission的研究。看了不少文献（如KOUTMOS and Booth (1995)），想自己动手操作一下，无奈怎么样都找不到多变量EGARCH模型的code，自己编了半天也找不到方向。请问谁有或者写过，借小妹观摩一下:)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Multivariate multivariat Variate EGARCH multi Multivariate EGARCH code

楼主弄的怎么样了，我最近也在做这个～
想向楼主讨教一下～

ifallasleepjing 发表于 2010-1-15 10:13
楼主弄的怎么样了，我最近也在做这个～
想向楼主讨教一下～

Here is the univariate case using simulated data with SAS

%let var2 = 2.5;
%let nobs = 1000 ;
%let arch0 = 0.1 ;
%let arch1 = 0.2 ;
%let garch1 = 0.75 ;
%let intercept = 0.5 ;


data egarch ;
   lu = &var2 ;
   lh = &var2 ;
   theta = .65 ;
   lz = lu/sqrt(lh) ;
   lg = theta*lz + abs(lz)-sqrt(2/3.14159) ;
   do i = -500 to &nobs ;
          /* EGARCH */
      h = exp( &arch0 + &arch1*lg + &garch1*log(lh)) ;
      u = sqrt(h)*rannor(12346) ;
      z = u/sqrt(h) ;
      g = theta*z + abs(z) -sqrt(2/3.14159) ;
      y = &intercept + u ;
      lu = u ;
      lh = h ;
      lg = g ;
      if i > 0 then output ;
   end ;
run ;   


/* Estimate EGARCH Model with PROC AUTOREG */
proc autoreg data= egarch ;
   model y =  / garch=( q=1, p=1 , type = exp) ;
run;
quit;

/* Estimate EGARCH Model with PROC MODEL */     
proc model data = egarch ;
   parms  earch0 .1 earch1 .2 egarch1 .75 theta .65 ;
   /* mean model */
   y = intercept ;
   /* variance model */
   if (_obs_ = 1 ) then
      h.y = exp( earch0 + egarch1 * log(mse.y)  );
   else h.y = exp(earch0 + earch1*zlag(g) + egarch1*log(zlag(h.y))) ;
   g = theta*(-nresid.y) + abs(-nresid.y) - sqrt(2/constant('pi')) ;
   /* fit the model */
   fit y / fiml method = marquardt ;
run;
quit;

谢谢呢，还想问问多元的怎么办呢？