全コードは右上のCodeを参照されたい。
各資産のリスクは(収益率の標準偏差)、以下のように求められる。
\[ Var[r_A] = E[r_A^2] - E[r_A]^2 \\ = 0.0101 - 0.01^2 \\ = 0.01 \\ \]
\[ SD[r_A] = \sqrt{Var[r_a]} \\ = \sqrt{0.01} \\ = 0.1 \]
\[ Var[r_B] = E[r_B^2] - E[r_B]^2 \\ = 0.0153 - 0.03^2 \\ = 0.0144 \\ \]
\[ SD[r_A] = \sqrt{Var[r_A]} \\ = \sqrt{0.0144} \\ = 0.12 \]
\[ Cov[r_A, r_B] = E[r_Ar_B] - E[r_A]E[r_B] \\ = 0.0087 - 0.01 \times 0.03 \\ = 0.0084 \]
\[ SD[r_A, r_B] = \sqrt{Cov[r_Ar_B]} \\ = \sqrt{0.0084} \\ = 0.0916515 \]
\[ Cor[r_A, r_B] = \frac{Cov[r_A, r_B]}{\sqrt{Var[r_A]Var[r_B]}} \\ = \frac{Cov[r_A, r_B]}{SD[r_A]SD[r_B]} \\ = \frac{0.0084}{0.1 \times 0.12} \\ = 0.01008 \]
p_mu <- 0.7*ra_mu + 0.3*rb_mu
p_var <- 2*0.7*0.3*rab_cov + 0.7^2*ra_var + 0.3^2*rb_var
p_sd <- sqrt(p_var)\[ E[r_p] = w_AE[r_A] + w_BE[r_B] \\ = 0.7\times0.01 + 0.3\times0.03 \\ = 0.016 \]
\[ Var[r_p] = w_Aw_ACov[r_A, r_A] + w_Aw_BCov[r_A, r_B] \\ + w_BW_ACov[r_B, r_A] + w_Bw_BCov[r_B, r_B] \\ = w_A^2Var[r_A] + 2w_Aw_BCov[r_A, r_B] + w_B^2Var[r_B] \\ = 0.7^2\times0.01 + 2\times 0.7\times0.3 \times 0.0084 + 0.3^2 \times 0.0144\\ = 0.009724 \]
\[ SD[r_p] = \sqrt{Var[r_p]} \\ = \sqrt{0.009724} \\ = 0.0986103 \]
ab_var <- a_var + b_var
ab_cov <- (2^2*a_var + 3^2*b_var - a2_b3_var) / 2
ab_cor <- ab_cov / sqrt(a_var * b_var)まず、以下の式を\(Cov[X_A, X_B]\)について解く。 \[ Var[2X_A + 3X_B] = Var[2X_A] + Var[3X_B] + 2Cov[2X_A, 3X_B] \\ Cov[2X_A, 3X_B] = \frac{Var[2X_A] + Var[3X_B] - Var[2X_A + 3X_B]}{2} \\ = \frac{2^2Var[X_A] + 3^2Var[X_B] - Var[2X_A + 3X_B]}{2} \\ = \frac{2^2\times600 + 3^2\times800 - 10^{4}}{2} \\ = -200 \qquad...(1) \]
(1)より、相関係数は、
\[ Cor[X_A, X_B] = \frac{Cov[X_A, X_B]}{\sqrt{Var[x_A]Var[X_B]}} \\ = \frac{-200}{\sqrt{600 \times 800}} \\ = -0.2886751 \\ \approx -0.29 \]
# 株式A
n_a <- 200
# 株式B
n_b <- 300
a_sd <- sqrt(a_var)
b_sd <- sqrt(b_var)
# 期待値
mu_x <- n_a * a_mu + n_b * b_mu
# 標準偏差
var_x <- 100^2 * a2_b3_var
sd_x <- sqrt(var_x)
# kの値
k_low <- abs(20000 - mu_x) / sd_x
k_high <- abs(80000 - mu_x) / sd_x
# チェビシェフの不等式の答
cheb1 <- 1/k_low^2
cheb2 <- 1/k_high^2\(\mu\) \[ E[X]=E[200X_A+300X_B]\\ =E[200X_A]+E[300X_B]\\ =200E[X_A]+300E[X_B]\\ =200\times70 + 300\times 120 \\ =5\times 10^{4} \]
\(\sigma\)
\[ Var[X]=Var[200X_A+300X_B]\\ =Var[100(2X_A+3X_B)]\\ =100^2Var[2X_A+3X_B]\\ =100^2 \times 10^{4} \\ =10^{8} \] \[ SD[x] = \sqrt{Var[X]} \\ = 10^{4} \]
\[ P(|X-\mu|\geq k\sigma) \leq \frac{1}{k^2} \\ P(|20000-5\times 10^{4}|\geq k\times10^{4}) \leq \frac{1}{k^2} \\ \] \(k \leq 3\)より、
\[ P(X \leq 20000) \leq \frac{1}{(3)^2} \\ P(X \leq 20000) \leq 0.1111111 \]
\[ P(|X-\mu|\geq k\sigma) \leq \frac{1}{k^2} \\ P(|80000-5\times 10^{4}|\geq k\times10^{4}) \leq \frac{1}{k^2} \\ \]
\(k \leq 3\)より、 \[ P(X \geq 80000) \leq \frac{1}{(3)^2} \\ P(X \geq 80000) \leq 0.1111111 \]
したがって、期待値\(\times\)株数を資産価値とした場合、このポートフォリオの資産価値が20000円を下回るか80000円を上回る確率は、0.2222222となる。
\[ P(X \leq 20000) + P(X \geq 80000) \leq 0.1111111 + 0.1111111 = 0.2222222 \]