《金融数据分析导论：基于R语言》习题答案（第一章）



《金融数据分析导论：基于R语言》是芝加哥大学的教授Ruey S.Tsay所著，李洪成、尚秀芬、郝瑞丽翻译，机械工业出版社出版，是一本学习R语言和金融数据分析的很好的参考书籍。

** 注：这些答案都是本人自己做出的结果，可能有错，仅供参考，发现有错的地方欢迎大家指出。 **

** **

** 第一章 **

1.

首先，将数据包放在当前工作目录下

library(fBasics)

da = read.table("d-axp3dx-0111.txt",header=T) % 读出数据

head(da) % 显示数据的前 6 行，可以观察数据格式

mmm = da[,2:5] % 取出简单收益率的数据，即把日期去掉

basicStats(mmm)

对 mmm 中的数据做基础分析，就可得到各列简单收益率序列的样本均值 (Mean) 、标准差 (Stdev) 、偏度 (Skewness) 、超额峰度 (Kurtosis) 、最大值 (Maximum) 、最小值 (Minimum)

rs = log(mmm+1) % 根据简单收益率求出对数收益率，公式为 r=In(R+1)

basicStats(rs) % 基础分析求出一系列统计值

ss = rs[,1]% 取出 AXP 股票的对数收益率

t.test(ss)

p 值大于 0.05 ，接受原假设，即 AXP 对数收益率的均值为 0

2.

首先，将数据包放在当前工作目录下

library(fBasics)

da = read.table("m-ge3dx-4011.txt",header=T) % 读出数据

head(da) % 显示数据的前 6 行，可以观察数据格式

mmm = da[,2:5] % 取出简单收益率的数据，即把日期去掉

basicStats(mmm)

提示错误如下：

经检查发现在数据的第 58 行（不包括标题） 19441031 的 sp 数据为 ”.” （非数值数据），因此系统报错。通过上网查阅资料，只查到如何删除含有 ”NA” 的行数据： na.omit(mmm) 或 mmm[complete.cases(mmm), ] [1] ，所以手工删除第 58 行：

mmm = mmm[-58,]%x = x[-m,] 表示删除 x 矩阵中的第 m 行

然而删除了还是报同样的错误，于是我手工去 txt 文件里删除了第 58 行，就 ok 了，不知道程序怎么改才可以 = =

然后继续

对 mmm 中的数据做基础分析，就可得到各列简单收益率序列的样本均值 (Mean) 、标准差 (Stdev) 、偏度 (Skewness) 、超额峰度 (Kurtosis) 、最大值 (Maximum) 、最小值 (Minimum)

rs = log(mmm+1) % 根据简单收益率求出对数收益率，公式为 r=In(R+1)

basicStats(rs) % 基础分析求出一系列统计值

ss = rs[,1]% 取出 AXP 股票的对数收益率

t.test(ss)

p 值小于 0.05 ，拒绝原假设，即 AXP 对数收益率的均值不为 0

3.

首先，将数据包放在当前工作目录下

library(fBasics)

da = read.table("m-ge3dx-4011.txt",header=T) % 读出数据

head(da) % 显示数据的前 6 行，可以观察数据格式

mmm = da[,5] % 取出 S&P 月收益率的数据

t.test(mmm) % 做 t 检验

因为 p 值为 2.436e-05<0.05 ，拒绝原假设，即 S&P 综合指数的月股票收益率均值不为 0.

s3 = skewness(mmm)

T = length(mmm)

t3= s3/sqrt(6/T)

t3

t3 为 -7.057059 ，绝对值大于1.96 ，故拒绝原假设，即 S&P 综合指数的月股票收益率偏度不为 0.

s4 = kurtosis(mmm) % 计算出超额峰度

t4 = s4/sqrt(24/T)

t4

t4>1.96 ，故拒绝原假设，即 S&P 综合指数的月股票收益率的峰度不为 3.

4.

首先，将数据包放在当前工作目录下

library(fBasics)

da = read.table("d-axp3dx-0111.txt",header=T) % 读出数据

head(da) % 显示数据的前 6 行，可以观察数据格式

mmm = da[,2] % 取出第 2 列 axp 的日简单收益率的数据

rs = log(mmm+1) % 求出对数收益率

s3 = skewness(rs)

T = length(rs)

t3= s3/sqrt(6/T)

t3

t3 的绝对值小于 1.96 ，故接受原假设，即 axp 日对数收益率的偏度度量等于 0.

s4 = kurtosis(rs) % 计算出超额峰度

t4 = s4/sqrt(24/T)

t4

t4>1.96 ，故拒绝原假设，即 axp 日对数收益率的月股票收益率的超额峰度不为 0.

5.

首先，将数据包放在当前工作目录下

library(fBasics)

da = read.table("d-fx-ukus-0711.txt",header=T) % 读出英镑对美元汇率数据

head(da) % 显示数据的前 6 行，可以观察数据格式

发现最后一列是英镑对美元的汇率，取出最后一列数

mmm = da[,4]

将其转换为美元兑英镑的汇率

mmm = 1/mmm

根据每天的汇率值求出每天的日对数收益率，rt = In(pt/p(t-1)) ，rt 是第 t 期的日对数收益率，pt 是第 t 期的汇率值，设第 1 期r1 = In(pt/p(t-1)) = In1 = 0 ，循环计算从第 2 期到第 length(mmm) 期的日对数收益率。

r1 = 1:length(mmm) % 定义一个数组表示日对数收益率

r1[1] = 0 % 设定第 1 天的日对数收益率为 0

for (i in 2:length(r1)) {

r1[i]=log(mmm[i]/mmm[i-1])

}

r1 即为美元对英镑的汇率从 20070102-20111130 的日对数收益率。

同样的办法求出美元对日元的汇率从 20070102-20111130 的日对数收益率。

da = read.table("d-fx-usjp-0711.txt",header=T)

head(da)

r2= 1:length(mmm)

r2[1] = 0

for (i in 2:length(r2)) {

r2[i]=log(mmm[i]/mmm[i-1])

}

r2 即为美元对日元的汇率从 20070102-20111130 的日对数收益率。

对 r1 和 r2 分别作基础分析，即可得到它们的样本均值 (Mean) 、标准差 (Stdev) 、偏度 (Skewness) 、超额峰度 (Kurtosis) 、最大值 (Maximum) 和最小值 (Minimum) 。

对于数据集 r2 ，求其密度分布的公式是 d = density(r2) ，画图函数为 plot(d$x, d$y, xlab=’log- return’, ylab=’density’, type=’l’)( 这两个公式见书中第 35 页的代码 )

d = density(r2)

plot(d$x,d$y,xlab='log-return',ylab='density',type='l')

得到密度图如下：

假设检验（ t 检验）：

t.test(r2)

p 值大于 0.05 ，接受原假设，即美元 / 日元汇率的日对数收益率的均值为 0.

[1] 人大经济论坛 http://bbs.pinggu.org/thread-2331183-1-1.html

吐槽： R 的矩阵定义好复杂啊