某ボラティリティ指数のスポット・先物(1か月後満期)の相関

ソースは御察しいただくとして、チャート観た感じ、すんげぇずれてるように見えたんだけど意外と相関してたわ。結構ヘッジできそう。

データ(df1, df2)は某所より入手してこんな感じで変形（リターン化）。データは直近250日分データで。

x <- df1[,1:2]
y <- df2[,1:2]
colnames(x) <- c("date", "value")
colnames(y) <- c("date", "value")
df <- cbind(x[(nrow(x)-250):nrow(x),], y[(nrow(y)-250):nrow(y),])
df <- df[,c(1,2,4)]
iv <- as.data.frame(apply(df[,2:3], 2, function(x)diff(x)/x[-length(x)]))
colnames(iv) <- c("v1", "v2")

相関、0.8位あるので結構イケてる。

> # 相関
> cor(iv)
          v1        v2
v1 1.0000000 0.7941028
v2 0.7941028 1.0000000

回帰してチャートにしてみる。

> # カイキー
> lm_iv <- lm(v2~v1+0, data=iv)
> plot(iv, xlim=c(-0.3, 0.3), ylim=c(-0.1, 0.1))
> abline(lm_iv,col=2)

ボラティリティの期間構造のパラレル成分だけをヘッジできるイメージなのかなぁ。