(完整word版)如何用stata做事件研究(word文档良心出品)

1、用 stata 做事件研究时间研究通常被用来检验市场对附带相应利益事件的反应。事件研究通常包括以下几 步：（1） 净化数据（提出无关和无法研究以及其他原因）和计算事件窗口（2） 估计正常表现（3） 计算异常表现和累积超额回报（4） 显著性检验（5） 全部事件交叉检验（相当于稳健性检验）这个文章旨在帮助你用 stata 来开展事件研究。假设你已经拥有一个时间变量（date）和 公司标识（company_id），如果你需要准备你的数据或者想要利用我们的数据尝试一下相应的命令，去data preparation页面。我们同时也假设你对 stata 有一个基本的了解。如果你需要和 stata 命令相关

2、的帮助，你 可以从 here 找到更多。如果你从一个可执行文件中进入命令，你的任务就会简单的多，因 为可执行文件中包含了一系列的 stata 命令。净化数据并计算事件窗口和估计窗口你很可能获取了超出你需要的每个公司的观察值，也有可能有一些公司的观察值不充 分。在你开展下一步前，你必须确保你的分析是建立在正确的观察值之上。为了实现这一点， 你需要设立一个变量 dif ，它将会计算从观察（相当于估计期）到事件期的天数，这即有 可能是节假日，也有可能是交易日。就交易日天数：sort company_id dateby company_id: gen datenum=_nby company_id:

3、gen target=datenum if date=event_dateegen td=min(target), by(company_id)drop targetgen dif=datenum-td就节假日：gen dif=date-event_date由上可以看出，计算交易日天数比计算节假日天数稍微复杂一点。对于交易日，我们首 先需要设立一个变量，这个变量是用来计算每一个公司标识范围内的天数，然后我们再确定 哪些观测值发生在事件期内。我们设立一个和事件天数有关的变量。这个事件天数建立在 company_id 范围内的所有 观察值的基础上。最后，我们选取二者之间的差异来设立一个变量 dif

4、，以计算每一个个体 观测值和事件期之间的天数。下一步，我们需要确定我们事件期前后的最小观察期天数以及 事件窗口之前的估计窗口的最小观察期天数。比如我们想要一个事件期前后两天的窗口（总 共 5 天的事件期）以及一个 30 天的估计窗口（你可以改变这些数字以适应你的分析）。by company_id: gen event_window=1 if dif=-2 & dif=2egen count_event_obs=count(event_window), by(company_id)by company_id: gen estimation_window=1 if dif=-60egen coun

5、t_est_obs=count(estimation_window), by(company_id)replace event_window=0 if event_window=.replace estimation_window=0 if estimation_window=.确定事件窗口和估计窗口的程序是相同的。首先，我们设立一个变量，当观测值在特定 的日期内等于 1；其次，我们设立另一个变量来计算，在每一个公司标识范围内，有多少观 测值的变量值等于 1；最后，我们用 0 代替所有的缺省值，设立一个哑变量。此时，你就可 以确定哪个公司没有充分的观测值。tab company_id if c

6、ount_event_obs5tab company_id if count_est_obs30tab 命令将会产生一个公司标识符的列表，这个列表中包含的就是没有足够的事件窗口 及估计窗口观察值窗口的公司，同时还会报告这些公司的总的观测值的数量，为了排除这些 公司使用如下命令：drop if count_event_obs 5drop if count_est_obs 30你必须确保在删除任何观察值之前，样本已经用另一个名字保存。此时，你也可以剔除一些你不需要的变量。估计正常表现现在我们到了真正分析的地方。首先，我们需要一种途径和方法来估计正常表现。为了实现这一点，我们将会利用估计窗口中的数据

7、来对每个公司进行单变量回归，并保存a（常数项）和 b （解释变量的系数），随后我们会使用这个保存的回归系数来预测事件窗口的正 常表现。就回报这个回归中的被解释变量来讲，它仅仅是针对既定股票回报的 crsr 变量，而 我们用来预测 ret 的解释变量，使一个对于任何股票的指数加权回报。对你的样本集使用等 价变量。set more off /* this command just keeps stata from pausing after each screen of output */gen predicted_return=.egen id=group(company_id)/* for m

8、ultiple event dates, use: egen id = group(group_id) */forvalues i=1(1)n /*note: replace n with the highest value of id */l id company_id if id=i & dif=0reg ret market_return if id=i & estimation_window=1predict p if id=ireplace predicted_return = p if id=i & event_window=1drop p这儿，我们设立了一个变量 id，以计算公司

9、数量（从 1 到实际数）。n 代表的是包含完 整数据的公司事件联合体数目（即事件有关的有完整数据的公司的数目）。这个过程通 过在公司间重复，并对每一个公司在估计窗口开展回归，然后用回归结果来预测事件窗口的 正常回报。异常表现和累积超额回报现在我们可以用我们的数据计算异常表现和累积超额回报。日异常回报等于事件窗口中 的每天的实际回报减预期回报。整个事件窗口的超额回报的总额就是累积超额回报。sort id dategen abnormal_return=ret-predicted_return if event_window=1by id: egen cumulative_abnormal_ret

10、urn = sum(abnormal_return)此处我们仅仅计算了事件窗口中每个样本观测值的超额回报，随后我们将 car 等于所 有公司的超额回报之和。显著性检验我们将计算一个检验统计量（t 值）来确定每支股票平均的超额回报是否显著不等于 0。 ar 即超额回报，ar_sd 是超额回报的标准差。如果 t 统计量的绝对值大于 1.96，那么平均 的超额回报在 5%的水平显著不等于 0.1.96 这个值来自于均值为 0，标准差为 1 的标准正态 分布，其中正态分布 95%的部分分布在1.96 之间。sort id dateby id: egen ar_sd = sd(abnormal_retu

11、rn)gen test =(1/sqrt(number of days in event window) * ( cumulative_abnormal_return /ar_sd)list company_id cumulative_abnormal_return test if dif=0note: this test uses the sample standard deviation. a less conservative alternative is to use the population standard deviation. to derive this from the sample standard deviation produced by stata, multiply ar_sd by the square root of n-1/n; in our example, by the square root of 4/5.这会将你的事件研究的结果生成一项 excel 文件。outsheet company_id event_date cumulative_abnormal_return test using stats.csv if dif=0, comma name全部事件交叉检验除了观察每个公司的平均超额回报，你可能还想计算将所有公司视为一个整