2.33σ
2σ
正規乱数の生成方法（一様乱数から作る方法）
様分布
（ⅰ） 一様乱数を作る（右図 ）。
Ｒａｎｄ（）
1
一様分布
×
Ｒａｎｄ（）
： ０以上で１より小さい乱数を発生させる。
（ⅱ） 一様乱数を標準正規乱数 に変換する（下図 ）
0 1
×
（ⅱ） 様乱数を標準正規乱数 に変換する（下図
×
） Ｎｏｒｍｓｉｎｖ（Ｒａｎｄ（））： 一様乱数の値を、標準正規分布の「分布関数の逆関数」に 代入すると 標準正規乱数に変換される
代入すると、標準正規乱数に変換される。
標準正規分布 1
分布関数
×
確率密度関数
（ⅲ） 標準正規乱数を（ⅱ）×σ＋μにより、正規乱数～Ｎ（μ、σ
2
） に変換する。0
（ⅳ） 正規乱数の生成方法には、様々なものがあり、どの方法が優れているか 研究の対象となっている。上記方法は一例に過ぎない
２項分布 結果が２通りある試行（実験）をＮ回繰り返したとき
２項分布： 結果が２通りある試行（実験）をＮ回繰り返したとき、２通りの結果のうち一方が起こる回数の確率分布
（例）サイコロを10回振って １の目が出る回数（Ｋ）
０回 ｆ（0）= 10 Ｃ 0 （１/６） 0 （５/６） 10 ・ ・ ・ ・ ・ ・ １回 ｆ（1）= 10 Ｃ 1 （１/６） 1 （５/６） 9
２回 ｆ（2）= 10 Ｃ 2 （１/６） 2 （５/６） 8 10回 ｆ（10）= 10 Ｃ 10 （１/６） 10 （５/６） 0
ｆ（Ｋ） 確率 Ｆ（Ｋ） 分布関数（累積確率）
・ ・ ・ ・ ・ ・
0.8 1 0.4
ｆ（Ｋ） 確率 Ｆ（Ｋ） 分布関数（累積確率）
N 10 １ /6 N 10 １ /6
0.4 0.6 0.2
N=10, ｐ = １ /6 N=10, ｐ = １ /6
0 0.2
0
0 2 4 6 8 10
Ｋ
0 2 4 6 8 10Ｋ
0 2 4 6 8 10 0 2 4 6 8 10
Ｋ
１の目が出る回数 １の目が出る回数
（例） ＶａＲを超過する損失が発生する回数（Ｋ）
（例） ＶａＲを超過する損失が発生する回数（Ｋ）
ＶａＲを超過する確率 ｐ ＝ １ ％
ＶａＲを超過しない確率 １－ｐ ＝ 99％（信頼水準）
ＶａＲを超過しない確率 １－ｐ ＝ 99％（信頼水準）
ＶａＲの計測個数 Ｎ＝250
発生確率 ｆ（Ｋ） Ｃ （0 01）
Ｋ
（0 99）250 Ｋ
0 4
ｆ（Ｋ） 確率 Ｆ（Ｋ） 分布関数（累積確率）
発生確率 ｆ（Ｋ） ＝
250
ＣＫ
（0.01）Ｋ
（0.99）250－Ｋ
0 6 0.8 0.4 1
N=250,ｐ=１％ N=250, ｐ = １％
0 2 0.4 0.6 0.2
0 0.2
0
0 2 4 6 8 10
Ｋ
0 2 4 6 8 10Ｋ
ＶａＲ超過損失の発生回数 ＶａＲ超過損失の発生回数
（４）確率変数の期待値
確率変数（Ｘ）は、平均的にみてどんな値をとるのか？（４）確率変数の期待値
（例）サイコロを振ったときに出る目の数 確率
P(X)
１／６
( )
１ ２ ３ ４ ５ ６ Ｘ
サイコロを振ったときに出る目の数の「期待値」
６
ＸＰ（Ｘ）
＝ １× （１／６） ＋ ２×（１／６） ＋ ３×（１／６）
Σ
X＝１１ （１／６） ＋ ２ （１／６） ＋ ３ （１／６）
＋ ４× （１／６） ＋ ５×（１／６） ＋ ６×（１／６）
＝ 3.5
例 変 率
（例） ＴＯＰＩＸの変化率（Ｘ）
確率密度関数 ｆ（Ｘ）
Ｘ Ｘ Ｘ
ＸX （現在値） X X
X
Ｘ
下落（－） 上昇（＋）
X
０
（現在値） X‐１
X‐２
X
‐３
ＴＯＰＩＸの変化率（Ｘ）の期待値
∫
＋∞Ｘ ｆ（Ｘ）ｄＸ（ ）∫
－∞
（５）確率変数の独立
（５）確率変数の独立
【定義】
確率変数 Ｘ、Ｙ が互いに影響されず、それぞれの確率分布にした がって値をとるとき、確率変数 Ｘ、Ｙは、互いに「独立」であるという。数式で表すと Ｐ（Ｘ＝ａ、Ｙ＝ｂ）＝Ｐ（Ｘ＝ａ）Ｐ（Ｙ＝ｂ）
【定理】
確率変数 Ｘ、Ｙ が互いに「独立」のとき、以下のことが 成り立つ。① 確立変数 ＸＹ の期待値は それぞれの確率変数の期待値の積になる
【定理】
① 確立変数 ＸＹ の期待値は、それぞれの確率変数の期待値の積になる。
Ｅ（ＸＹ）＝Ｅ（Ｘ）Ｅ（Ｙ）
② 確率変数 Ｘ＋Ｙ の分散は、それぞれの確率変数の分散の和に等しい。
② 確率変数 の分散は、それぞれの確率変数の分散の和に等し 。 Ｖ（Ｘ＋Ｙ）＝Ｖ（Ｘ）＋Ｖ（Ｙ）
③ 確率変数 Ｘ と Ｙ は無相関である。
ρ（Ｘ、Ｙ）＝０ （証明省略）
（例）サイコロを振ったときに出る目の数
１回目： Ｘ
１
＝ １、 ２回目： Ｘ２
＝ １ ３回目： Ｘ３
＝ ？回目
１
、 回目２
サイコロの目（Ｘ
３
） １ ２ ３ ４ ５ ６確 率 １／６ １／６ １／６ １／６ １／６ １／６
２回続けて１の目が出ても、３回目の結果には影響確 率 ／ ／ ／ ／ ／ ／
を及ぼさない。
ずれ が る確率も
３回目は、いずれの目が出る確率も１／６。（例）株価 金利 為替等リスクフ クタ の変化率
（例）株価、金利、為替等リスクファクターの変化率
過去の変化率（実績）が、将来の変化率（予想）に影響 を及ぼすことはないと考えて、互いに独立な確率変数と して捉えることが多い。リスクファクター（X）の推移と、その確率分布
X
０Ｘ Ｘ
Ｘｓ
０ Ｘ
Ｘ
Ｘｔ
？
ｔ０ Ｘｔ
現在 将来 47
過去
しかし リスクフ クタ の変化率が時点間で独立とは
しかし、リスクファクターの変化率が時点間で独立とは 限らず、相関関係が認められることも少なくないので 注意を要する。注意を要する。
－ 下図は、ＴＯＰＩＸ・日次対数変化率１期前の変化率との相関
をみたもの 独立の判定には 様々なタイムラグを置いて相関の
１期前
をみたもの。独立の判定には、様々なタイムラグを置いて相関の 有無をみる必要。
2 3 4
0 1
-4 -3 -2 -1 0 1 2 3 4
当期
-3 -2 -1
相関係数 0 037
相関係数 ρ ＝ 0.037 -4
４．推定と検定
（１）記述統計と推測統計
（２）推定
（３）検定
（３）検定
統 推 統
（１）記述統計と推測統計
記述統計 基本統計量の算定や図表 グラフを利用して
記述統計 ： 基本統計量の算定や図表、グラフを利用して 観測データが持つ特性を分析・記述する。均
（例）特定の集団（Ｎ人）の身長の平均と分散を計算する。
Ｘ ＝
Ｘ１＋Ｘ２＋・・・＋ＸＮ 平均
Ｎ
（Ｘ Ｘ）２＋（Ｘ Ｘ）２＋ ＋（Ｘ Ｘ）２ 分散
Ｖｐ ＝
Ｎ
（Ｘ１－Ｘ）２＋（Ｘ２－Ｘ）２＋・・・＋（ＸＮ－Ｘ）２ 分散
推測統計 ： 標本として集めた一部の観測データに基づき、母集団の特性について推測し、検証する。
母集団 特性 推測 、検証す 。
（例）任意に抽出したＮ人（標本）の身長を計測して、日本人 全体（母集団）の身長の平均と分散を推定する
Ｘ ＋Ｘ ＋ ＋Ｘ 平均
全体（母集団）の身長の平均と分散を推定する。
Ｘ ＝
Ｎ
Ｘ１＋Ｘ２＋・・・＋ＸＮ
（Ｘ１－Ｘ）２＋（Ｘ２－Ｘ）２＋・・・＋（ＸＮ－Ｘ）２ 分散（不偏標本分散）
Ｖａ ＝
Ｎ－１
１ ２ Ｎ
（注）上記定義（偏差平方和をＮ １で割る）による標本分散Ｖ については 理論上
（注）上記定義（偏差平方和をＮ－１で割る）による標本分散Ｖａについては、理論上、
「その期待値が母集団の分散となる」ことが知られている。Ｖａは母集団の分散を 偏りなく推定する統計量となるため、 「不偏標本分散」と言う。
（２） 推 定
母集団の確率分布、特性値は、誰にも分からない。標本 特性値から 集 特性値を統計的 推測する
（２） 推 定
母集団確率分布
標本の特性値から母集団の特性値を統計的に推測する。特性値 平均μ 標準偏差σ Ｖ Ｒ など ＶａＲ など.
× ×
× × ×
×
×
× 特性値
平均μ＊
×
×
×
×
×
× ×
× ×
× ×
×
× ×
× ×
標準偏差
σ
＊ ＶａＲ＊ など母集団 標本（実現値） 推定
（３） 検 定
一定の確率分布を前提にして推定した値について、値 確率 有意 が
（３） 検 定
その値をとる確率（有意水準α％）が十分に低いとき、
「偶然、珍しいことが起きた」と考えるのではなく、
「推定の際に置いた前提（帰無仮説） が誤 ていた
「推定の際に置いた前提（帰無仮説） が誤っていた」
と結論付ける。
× 真の確率分布
推定に利用した確率分布
② 推定の前提（確率分布）が
② 推定の前提（確率分布）が 誤っていたと結論付ける。
有意水準 α％
① 実現する確率が十分に低い と考えられることが起きた。
実現値
ＶａＲを超過する損失が発生する回数（Ｋ）とその確率 ＶａＲを超過する損失が発生する回数（Ｋ）とその確率
ＶａＲを超過する確率 ｐ ＝ １ ％
ＶａＲを超過しない確率 １－ｐ ＝ 99％（信頼水準）
ＶａＲを超過しない確率 １－ｐ ＝ 99％（信頼水準）
ＶａＲの計測個数 Ｎ＝250
発生確率 ｆ（Ｋ） Ｃ （0 01）
Ｋ
（0 99）250 Ｋ
発生確率 ｆ（Ｋ） ＝250
ＣＫ
（0.01）Ｋ
（0.99）250－Ｋ
0.4
2 項分布 Ｎ =250, ｐ = １％
0.2
0
Ｋ：ＶａＲ超過損失
0 2 4 6 8
10
Ｋ：ＶａＲ超過損失
の発生回数
バックテスト（２項検定）
バックテスト（２項検定）
観測データ数
250
Ｎ回 Ｎ回の観測で、Ｋ回、ＶａＲを超過する確率信頼水準
99% K N K
信頼水準
99%
１－信頼水準
1%
ｐ％ＶａＲ超過回数
(K回) 確率 確率
ＶａＲ超過回数 (K回以上）
２項分布
N C K p K (1-p) N-K
(K回) 確率 確率 (K回以上）
0
8.11% 100.00%
0回以上1 20.47% 91.89%
1回以上2 25.74% 71.42%
2回以上3 21.49% 45.68%
3回以上回以上4 13.41% 24.19%
4回以上5 6.66% 10.78%
5回以上6 2 75% 4 12%
6回以上6 2.75% 4.12%
6回以上7 0.97% 1.37%
7回以上8 0.30% 0.40%
8回以上9 0.08% 0.11%
9回以上10 0.02% 0.03%
10回以上11 0.00% 0.01%
11回以上12 0.00% 0.00%
12回以上13 0 00% 0 00%
13回以上13 0.00% 0.00%
13回以上14 0.00% 0.00%
14回以上15 0.00% 0.00%
15回以上バックテストは「検定」の考え方にしたがって行う
ＶａＲ計測モデルは正しい（帰無仮説）。
ＶａＲ超過損失の発生が、250回中、10回以上発生した。Ｖ Ｒ超過損失の発生が 250回中 10回以上発生する
ＶａＲ超過損失の発生が、250回中、10回以上発生する 確率は0.03％と極めて低い。
ＶａＲ計測モデルは誤っている（結論）２種類の過誤
「検定」では、次の2通りの「過誤」（エラー）が起きる可能性 がある したがって バックテストの結果も「過誤」（エラ ）２種類の過誤
がある。したがって、バックテストの結果も「過誤」（エラー）
を伴っている可能性がある点、注意を要する。
第１種の過誤（エラー）
本当は帰無仮説（ＶａＲ計測モデル）が正しいのに、
検定の結果、
帰無仮説（ＶａＲ計測モデル）が誤っていると結論付けてしまう。
第２種の過誤（エラー）
本当は帰無仮説（ＶａＲ計測モデル）が正しくないのに、
検定の結果、
帰無仮説（ＶａＲ計測モデル）が正しいと結論付けてしまう。
帰無仮説（ 計測 デ ） 結論付け まう。
真の確率分布 推定に利用した確率分布 ＝ 真の確率分布 推定に利用した確率分布 ＝
第１種の過誤
実現値 実現値
真の確率分布 推定に利用した確率分布 ＝
第２種の過誤 第２種の過誤
実現値
５．線形回帰分析
（１）線形回帰分析とは
（２）Ｅｘｃｅｌ分析ツールを利用した回帰分析
（３）チェック項目（決定係数、Ｐ値）
（１）線形回帰分析とは
Xｉ と Ｙｉ の間に 「直線的な比例関係」があることを前提に して、Xｉ と Ｙｉ の散布図の中の各点のなるべく近くに直線 を描くを描く。
Ｙ ｉ = ａＸ ｉ ＋ｂ＋ｅ ｉ Ｙ ｉ ａＸ ｉ ＋ｂ＋ｅ ｉ
変数 Ｙ を変数 Ｘ で説明する。
Ｙ ｉ ： 被説明変数（目的変数）
Ｘ ｉ ｉ ： 説明変数 ａ ： 回帰係数
ｂ ： 定数項（切片） （注）本例のように、説明変数が１つの場合、
単回帰分析という 説明変数が２つ以上
e i ： 残差 単回帰分析という。説明変数が２つ以上
の場合、重回帰分析という。
最小２乗法 最小２乗法
残差ｅ i i = Ｙ ｉ ｉ －ａＸ ｉ ｉ －ｂ
の２乗和を最小にするようにａ 、 ｂ
を推定する。それぞれの推定値をａ 、 ｂ
と表記する。Ｙ
＾ ＾
実測値
Ｙ Ｙ Ｙｉ
ｅ
ｉ＾
理論値
Ｙ
ａ
＾
＾ ＾ ｂ Ｙi＝ａ Ｘi＋ｂ
61
Ｘ Ｘ
（２）Ｅｘｃｅｌ分析ツ ルを利用した回帰分析
（２）Ｅｘｃｅｌ分析ツールを利用した回帰分析
【手順】
①「ツール」メニューから「分析ツール」を起動。
②ボックスの中の「回帰分析」を選択してＯＫをクリック。
③「入力Ｙ範囲」 「入力Ｘ範囲」に それぞれデータ範囲を入力
③「入力Ｙ範囲」、「入力Ｘ範囲」に、それぞれデータ範囲を入力。
チェックを入れると観測値、
残差のグラフ等をを表示
（注）ＰＣによっては、分析ツール のアドインが必要です。
（例）Ｅ ｃｅｌ分析ツ ル 回帰分析の出力結果
（例）Ｅｘｃｅｌ分析ツール・回帰分析の出力結果
概要
回帰統計
重相関 R 0 956320779
X 値 1 観測値グラフ 重相関 R 0.956320779 0 25
重決定 R2 0.914549432 補正 R2 0.90844582 標準誤差 0.022258115
観測数 16
分散分析表 0
0.05 0.1 0.15 0.2 0.25
Y Y
予測値 : Y 分散分析表
自由度 変動 分散 観測された分散比 有意 F
回帰 1 0.074233006 0.074233006 149.8374126 7.24E-09
残差 14 0.006935932 0.000495424
合計 15 0.081168938
係数 標準誤差 t P-値 下限 95% 上限 95% 下限 95.0% 上限 95.0%
-0.05 0
0 0.2 0.4 0.6 0.8
X 値 1
係数 標準誤差 t 値 下限 95% 上限 95% 下限 95 % 上限 95 %
切片 -0.047846512 0.013516678 -3.539813066 0.003266347 -0.07684 -0.018856096 -0.076836928 -0.018856096 X 値 1 0.37369024 0.03052823 12.24080931 7.24475E-09 0.308214 0.439166839 0.308213641 0.439166839
残差出力
観測値 予測値 : Y 残差 標準残差
1 -0.027293549 0.028293549 1.315772009 2 -0.023182956 0.024182956 1.124611728 3 0.009328095 -0.008328095 -0.387292319 4 0.051555092 -0.050555092 -2.351029759 5 0 104619106 -0 011619106 -0 540338532
X 値 1 残差グラフ
0.02 0.04 5 0.104619106 -0.011619106 -0.540338532
6 0.092287328 0.006712672 0.312168184 7 0.097145301 0.001854699 0.086251488 8 0.097145301 0.001854699 0.086251488 9 0.108729699 -0.009729699 -0.452472943 10 0.117698264 -0.018698264 -0.869549921 11 0.12629314 -0.02729314 -1.269248692
-0.06 -0.04 -0.02 0
0 0.2 0.4 0.6 0.8
X 値 1
残差
63
12 0.175993942 -0.003993942 -0.185735522 13 0.177862393 0.018137607 0.843476924 14 0.167399066 0.028600934 1.330066732 15 0.176367632 0.019632368 0.912989753