- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Eviews编程
事件研究法EVIEWS编程
1. 编程前须知
编程过程用Reviews实现,编程重点在于循环语句和IF语句的使用,下面简单介绍循环语句的使用说明。
循环语句:
(1)循环语句块由关键字for开始,到关键字next结束,配对
(2)循环语句块以关键字for开头的语句称为循环控制语句,该语句由for开始,然后是控制变量的初始化。
例如:
For !i=1 to 41
~~~~~
~~~~~
Next
(3)for循环语句的执行过程为,判断循环控制变量是否穿越终止值。
(4)for循环允许嵌套
IF语句:
If语句以关键字if开始,end结束
例 if then
End
(2)If语句中If到Then之间的表达式称为If条件语句,只有当If条件语句为真时,If语句块才被执行。
@str():将控制变量的数值转化成字符串
@elem(a,b):取序列在给定时间的观测值,a一个序列,b指定的时间观测值,必须是字符串的形式。
Smpl %a=1 %b=100
@coefs(1)
@coefs(2)
Vector(41) car01
Matrix(41,3) nt01
@cnorm(J1/J2)
以估计窗口长度为90,事件窗(事后窗)为41,公司数220为例,给出编程过程及解释(一共收集了事件日前121天的数据,和事件日后61天的数据,每个样本含有183个观测值)。先对数据进行匹配处理,通过国泰安数据库,得到两个系列ser01和ser02,分别表示单个公司的收益率和市场综合的日收益率。Ser01、ser02共有220*183个数据。
2 估计窗 91 112 事件窗 152
事件日
编程过程如下:
第一步,提取系数
for !i=1 to 220 // 循环语句,以for开头,next结尾
%a=@str(183*(!i-1)+2) //将控制变量的数值转变成字符串变量
%b=@str(183*(!i-1)+91) //选取每个公司的估计窗口
smpl %a %b //设定样本区间
equation eq.ls ser01 c ser02 //市场模型回归,估计参数
a(!i)=@coefs(1) //提取系数的函数
b(!i)=@coefs(2) //提取市场模型的系数
next //循环结束
第二步,估计异常收益率及标准差
for !i=1 to 220
%a=@str(183*(!i-1)+112)
%b=@str(183*(!i-1)+152) //设定事件窗口
smpl %a %b
r01=ser01-(a(!i)+b(!i)*ser02) //r01存储异常收益率,格式41*220
sd(!i)=@stdev(r01) //存储每个公司的标准差
next //循环结束
第三步,对异常收益率进行横截面平均,仅计算时间窗口的,存储在向量ar01中
for !i=1 to 41 //
count=0 //计数器
for !j=1 to 220 //嵌套循环语句
count=count+@elem(r01,@str(41*(!j-1)+!i+1))
next //嵌套循环结束
ar01(!i)=count/220 //计算平均异常收益率
next //总循环结束
for !i=1 to 41
car01(!i)=ar01(!i)
if !i1 then
car01(!i)=ar01(!i)+car01(!i-1) //计算累计平均异常收益率
end if
next
第四步,计算J1统计量
for !i=1 to 220
for !j=1 to 41
r02(41*(!i-1)+!j)=@elem(r01,@str(41*(!i-1)+!j))/sd(!i) //r02存储标准化后异常收益率
next
next //标准化过程
for !i=1 to 41
count=0
for !j=1 to 220
count=count+r02(41*(!j-1)+!i) //标准化的平均异常收益路
next
ntest01(!1,1)=(@sqrt(90-4)/@sqrt(220*(90-2)))*count //J1统计量
next
for !i=1
文档评论(0)