By DevinChu, 2018/12/24
1 两种收益率的计算
例如 我们有下述股票的日度价格序列:
股票id | 日期 | 价格 |
---|---|---|
1 | 20180101 | 15 |
1 | 20180102 | 19 |
1 | 20180103 | 17 |
1 | 20180104 | 16 |
1 | 20180105 | 20 |
1 | 20180106 | 18 |
1 | 20180107 | 19 |
根据简单收益率和对数收益率计算公式:
简单算术收益率 EQ1:
连续对数收益率 EQ2:
我们可以计算出简单收益率和对数收益率如下:
股票id | 日期 | 简单收益率 | 对数收益率 | 1+简单收益率 |
---|---|---|---|---|
1 | 20180102 | 0.2667 | 0.2364 | 1.2667 |
1 | 20180103 | -0.1053 | -0.1112 | 0.8947 |
1 | 20180104 | -0.0588 | -0.0606 | 0.9412 |
1 | 20180105 | 0.2500 | 0.2231 | 1.2500 |
1 | 20180106 | -0.1000 | -0.1054 | 0.9000 |
1 | 20180107 | 0.0556 | 0.0541 | 1.0556 |
2 区间内收益率的计算
在实证研究中,很多时候我们希望得到周收益率、或者月度收益率,那么如何根据日度收益率得到这些指标呢?相应很多同学直觉上认为直接加总日收益率即可,那直接对简单收益率进行加总得到的周(月)收益率等于对数日收益率的加总吗?我们可以以上述数据来测试看一下(计算 2018 年 1 月 2 日至 2018 年 1 月 7 日的持有期收益率):
简单收益率的加总:
对数收益率的加总:
两者差值:
为什么两者会相差 0.0718 个单位呢?那么哪种计算方法才能准确反应投资者在 2018 年 1 月 2 日至 2018 年 1 月 7 日 期间内的持有收益呢?正确答案是对数日收益率的加总,即 EQ3:
因为区间内简单日收益率的加总意味着投资者每天都在进行交易,而非在区间内买入并持有至末期。当然,我们也可以通过简单算术日收益率计算区间内的持有期收益,但是需要采用下述方法:
持有期收益率 EQ4:
或者采用 EQ5:
3 Stata命令介绍:ascol
3.1 安装 ascol
在 Stata 命令框中输入下述命令并回车,即可实现安装:
ssc install ascol, replace
3.2 如何调用该命令进行收益率计算?
ascol
的完整命令格式如下:
ascol varname, returns([simple | log]) prices keep([all | vars]) frequency_options timevar(varname) panelvar(varname) generate(newvar)
其中:
- varname 主要是收益率序列或者价格序列。
- returns 中可选择 simple 或 log ,simple 代表varname是简单收益率,log 代表对数收益率。若 varname 设定的是价格序列,则这里的选择会指导具体收益率序列的计算方式。
- prices 指定了 varname 中序列的类型是否为价格序列,如果是收益率序列则无需填写。
- keep选项可保留计算后的重复数据,例如日收益率转换为月收益率后,在相同月会出现多个相同的收益率观测,如果选择了 keep 命令则保留了相同的观测,反之则只每个时间点上保留一个观测,剔除了重复观测
- frequency_options 中可设定 toweek tomonth toquarter toyear,即由日收益率转换为不同频率的收益率。
- timevar 和 panelvar 选项可设定数据的个体截面id和时间截面id。
- generate 可设定新生产的收益率序列的名称。
3.3 举个栗子
下述代码来源于
ascol
的官方帮助文档,读者可自行前往查看
/*生产收益率序列*/
clear
set obs 1000
gen date=date("1/1/2012" , "DMY")+_n
format %td date
tsset date
gen pr=10
replace pr=pr[_n-1]+uniform() if _n>1
gen simpleRi=(pr/l.pr)-1
gen logRi = ln(pr/l.pr)
save stocks,replace
/*Example 1: From Daily to weekly - simple returns*/
use stocks, clear
ascol simpleRi, toweek returns(simple)
/*Example 2: From Daily to weekly - log returns*/
use stocks, clear
ascol logRi, toweek returns(log)
/* Example 3: From Daily to monthly - prices*/
use stocks, clear
ascol pr, tomonth price
3.4 补充一个新命令 mtoq
上述的ascol
命令虽然可以对日收益率进行不同频率的转换,但是有些时候我们下载的收益率本身就是月度收益率,那如何转换为季度收益率呢?我们写了一个新的命令mtoq
可实现这一功能。
Step1 安装mtoq
:
/*1.没有安装github接口的请先通过下述命令进行安装*/
net install github, from("https://haghish.github.io/github/") replace
/*2.安装过github接口的同学可忽略上述命令,直接进行下述安装*/
github install zhbsis/mtoq
Step2 调用mtoq
:
mtoq
命令非常简洁,如下:
mtoq varlist, by(string) [Statistics(string)]
其中,
- varlist 中填写对数收益率序列,可同时填写多个序列。
- by 选项中填写个体和时间截面的id,例如股票类数据可填写 stkcd(股票代码) ym(年月时间标记)
- Statistics 选项中可填写 Stata官方命令
egen
中的常见 option ,例如 mean、sum等,默认为 mean 处理,在月度收益率转换为季度收益率的时候通常填写 sum 选项。
Tips:当然,
mtoq
命令默认了由月度频率向季度频率转换,因此不局限于收益率序列的转换,可转换其他变量。例如当你有月度的换手率(TurnOver)数据时,你想要通过加总月度换收益率得到季度的换手率,可输入:mtoq TurnOver, by(stkcd ym) s(sum)