正規分布と対数正規分布について教えてください。統計学をかじり始めた初学者です。助けて下さい(´;ω;`)
例えば、元手1万円が正規分布に従う平均利率R・標準偏差sで増えていくとします。
そうすると、n=10年後のリターンVは自然対数eと各年でバラつくRnで書き表すと以下になるかと思います。
V= 1×e^R1 ×e^R2 ×e^R3 ×e^R4 ・・・×e^R10
ここで私の理解が間違っていなければ、 Rnは正規分布ですが、べき乗が入ったVは対数正規分布に従うのかなと理解しています。
実際、R=10%, s= 20%、Vを1万回シミュレーションした結果、
Vの標準偏差σ≒2.3, 平均値μ≒3.3 ・・・①
となり画像のように対数正規分布っぽい形になりました。
ちなみにExcelでは、e^RnをEXP(NORM.INV(RAND(),0.1,0.2)、標準偏差σはSTDEV(Vの1万回試行), 平均μはAVERAGE(Vの1万回試行)で出しました。
さて前置き長くなりましたが、ここからがわからないことです。
以下の参考URLによれば、対数正規分布に従う確率密度関数の平均値と標準偏差は以下の公式が成り立つようです。
http://valavg.com/risk-return-compound-1/
標準偏差σ= SQRT(LN((s/m)^2+1))
平均μ= LN(m)-LN((s/m)^2+1)/2
※m=1+R, LNはlogの底e、sは正規分布の標準偏差
で、この公式にm=1+R=1.1, s=0.2を代入すると、
σ=0.18, μ= 0.079 ・・・②
①≠②になってしまいます( ノД`)
これはエクセルによる数値計算の結果①と、公式計算の結果②が全然異なるということなので、どこかを根本的に勘違いしていると思いますが、何をやらかしているのかわかりません。
教えてください!!
よろしくお願いします。
A 回答 (7件)
- 最新から表示
- 回答順に表示
No.7
- 回答日時:
V= 1×e^R1 ×e^R2 ×e^R3 ×e^R4 ・・・×e^R10
は違うんじゃないですか?
毎年、平均利率R・標準偏差sの正規分布に従う利率で増えていくとすると、
V= 1 ×(1+R1) ×(1+R2) ×(1+R3) ×(1+R4) ×(1+R5) ×(1+R6) ×(1+R7) ×(1+R8) ×(1+R9) ×(1+R10)
ただし R1, R2, R3, R4, R5, R6, R7, R8, R9, R10 は N(R,s^2) に従う
...になるはずです。
指数関数が現れる理由がない。
No.6
- 回答日時:
#5です。
正規乱数の作り方がおかしいのではないでしょうか。と書きましたが、
NORM.INV(RAND(),0.1,0.2)という作り方は合っていました。でも、
V= 1×e^R1 ×e^R2 ×e^R3 ×e^R4 ・・・×e^R10
この式はおかしいですよね。
正しくは、
Ri=NORM.INV(RAND(),1.1,0.2)で(利率+1)を作って、
V= 1×R1 ×R2 ×R3 ×R4 ・・・×R10
でリターンVが出ます。これを1万回やればリターンVの分布が出ます。
これが非対称分布になるので、1万個のVの対数を取って対数正規分布に直します。その平均や標準偏差を調べれば良いと思います。
ただし、計算上、同じ利率が10回使用されてしまわないように、面倒でも、こまごまと式を書かなければなりません。私はforループで回しています。
No.5
- 回答日時:
#4です。
私のシミュレーション結果が合っているかどうかの確認として、シミュレーションで得られた平均値μ≒2.6の10乗根を取って1を引いたところ、
> mean(V)^(1 / 10) - 1
[1] 0.09918868
となり、1年分の利息0.1にほぼ等しかったです。
ご質問者様は3.3になったとのことですので、この検算を行うとシミュレーションは間違いだったと分かります。
分布形は合っていますので、正規乱数の作り方がおかしいのではないでしょうか。
No.4
- 回答日時:
#1です。
回答を間違えた責任を感じ、懺悔のつもりでシミュレーションをやってみました。そしたら、ご質問者さんの間違いに気づきました。
換算する年数を間違えています。
ご質問者が示された数値の、
・Vの標準偏差σ≒2.3, 平均値μ≒3.3 は10年分(しかし、このシミュレーション値は間違いでした)
・σ=0.18, μ= 0.079 は1年分です(これは合っていました)
実際のシミュレーション結果を以下に示します(R言語です)。
結果は、
10年分の利息は(ただし、ご質問者様に合わせる為、1は引いていません)
・Vの標準偏差σ≒1.6, 平均値μ≒2.6 (10年後)
この1万個の結果を全て自然対数に直すと、
・Vの標準偏差σ≒0.6, 平均値μ≒0.77 (10年分)
これを1年分に換算すると(標準偏差は分散の加法性から換算)
・Vの標準偏差σ≒0.19, 平均値μ≒0.077
となり、ご質問者が計算された理論値とほぼ一致しました。
> # 10年利息のシミュレーション
>
> r <- 0.1
> s <- 0.2
>
> V <- NULL
> i <- 0
> while(i < 10000){
+ V1 <- 1
+ for(j in 1:10){
+ V1 <- V1 * rnorm(1, (r + 1), s)
+ }
+ V <- append(V, V1)
+ i <- i + 1
+ }
>
> par(mfrow = c(2, 1))
> hist(V, breaks = 30, las = 1)
> abline(v = mean(V))
> hist(log(V), breaks = 20, las = 1)
> abline(v = mean(log(V)))
>
> mean(V)
[1] 2.574675
> sd(V)
[1] 1.606537
>
> mean(log(V))
[1] 0.7707654
> sd(log(V))
[1] 0.6011561
>
> mean(V)^(1 / 10) - 1
[1] 0.09918868
>
> mean(log(V)) / 10
[1] 0.07707654
> sqrt(var(log(V)) / 10)
[1] 0.1901023
No.1
- 回答日時:
企業で統計を推進する立場の者です。
参考URLの公式に、m=1.1、s=0.2 を代入されたそうですが、それは間違いです。mやsは対数正規分布の平均、標準偏差でなければなりません。
両者の分布間のパラメータの換算式を示した表(社内教育資料)を添付します。
もしm=1.1、s=0.2 を代入して対数正規分布のパラメータ理論値を求め、シミュレーション結果と比較するのであれば、一番左の式に代入して下さい。ただし、面倒なことに、一旦左下の式でSを求めてから、左上の式でEを求めることになります。
ちなみに参考URLの式は、一番右と同じです。
【表】対数化する前の平均μ、標準偏差σと、対数正規分布の平均E、標準偏差Sの関係
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 統計学 Excelによるサンプルの拡大について 6 2023/08/22 16:03
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- 数学 以下の数学の問題を教えてください。 確率変数Xは標準正規分布N(0、1)に確率変数Yは平均3のポアソ 3 2022/12/02 19:13
- 大学・短大 大学 統計学 1 2022/09/14 11:27
- 統計学 信頼区間についての質問です。 6 2023/06/25 17:34
- 統計学 以下の問題の解き方が分からないので式と使用した数字の求め方を教えてください 全国の中小企業の取締役か 8 2023/01/13 17:13
- 統計学 数学です 4 2022/07/20 23:06
- 数学 高校の数学Bの、確率分布と統計的な推測の、 正規分布の問題でわからない箇所がございます。問題文が、 2 2022/03/27 20:57
- 統計学 母集団分布を平均 μ, 分散 σ2 の正規分布と想定し, 母集団から無作為抽出した標本のデータ(標本 4 2023/01/30 20:25
- 統計学 統計学が分かりません!詳しい解説と回答を教えてくださる方お願いいします! 5 2022/08/23 03:10
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・【大喜利】【投稿~11/22】このサンタクロースは偽物だと気付いた理由とは?
- ・お風呂の温度、何℃にしてますか?
- ・とっておきの「まかない飯」を教えて下さい!
- ・2024年のうちにやっておきたいこと、ここで宣言しませんか?
- ・いけず言葉しりとり
- ・土曜の昼、学校帰りの昼メシの思い出
- ・忘れられない激○○料理
- ・あなたにとってのゴールデンタイムはいつですか?
- ・とっておきの「夜食」教えて下さい
- ・これまでで一番「情けなかったとき」はいつですか?
- ・プリン+醤油=ウニみたいな組み合わせメニューを教えて!
- ・タイムマシーンがあったら、過去と未来どちらに行く?
- ・遅刻の「言い訳」選手権
- ・好きな和訳タイトルを教えてください
- ・うちのカレーにはこれが入ってる!って食材ありますか?
- ・おすすめのモーニング・朝食メニューを教えて!
- ・「覚え間違い」を教えてください!
- ・とっておきの手土産を教えて
- ・「平成」を感じるもの
- ・秘密基地、どこに作った?
- ・【お題】NEW演歌
- ・カンパ〜イ!←最初の1杯目、なに頼む?
- ・一回も披露したことのない豆知識
- ・これ何て呼びますか
- ・初めて自分の家と他人の家が違う、と意識した時
- ・「これはヤバかったな」という遅刻エピソード
- ・これ何て呼びますか Part2
- ・許せない心理テスト
- ・この人頭いいなと思ったエピソード
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・好きなおでんの具材ドラフト会議しましょう
- ・餃子を食べるとき、何をつけますか?
- ・あなたの「必」の書き順を教えてください
- ・ギリギリ行けるお一人様のライン
- ・10代と話して驚いたこと
- ・大人になっても苦手な食べ物、ありますか?
- ・14歳の自分に衝撃の事実を告げてください
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・あなたの習慣について教えてください!!
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
統計学でいうRSD%とは何ですか。
-
正規分布に従わないと標準偏差...
-
相対標準偏差についてですが…
-
標準偏差について
-
幾何平均と幾何標準偏差について
-
正規分布について教えてくださ...
-
エクセルの標準偏差を用いた%計...
-
中央値の求め方ってどのように...
-
標準偏差
-
統計学における有効数字につい...
-
推定(統計)です。有効数字の...
-
偏差値60の人と偏差値50の人が...
-
今日数学のテストが返ってきて ...
-
(1)で分散を求める時、解答では...
-
工程能力のN数補正について
-
標準偏差バーをグラフに入れた...
-
統計でいう「n」は、何の略な...
-
パーセンテージのバラツキを表...
-
平均値と中庸値の違い
-
社会人の偏差値の調べ方
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
正規分布に従わないと標準偏差...
-
統計学でいうRSD%とは何ですか。
-
相対標準偏差についてですが…
-
中央値の求め方ってどのように...
-
幾何平均と幾何標準偏差について
-
算術平均に対する標準偏差は、...
-
標準偏差0.09ってどうですか?
-
変動係数は正規分布を前提?
-
ワイブル分布の標準偏差
-
貯蓄額の平均値と中央値の差額...
-
標準偏差の意味とは?
-
正規分布について教えてくださ...
-
対数正規分布
-
エラーバー(誤差棒)についてです。
-
困っています。標準誤差(エラ...
-
幾何正規分布を描くために、幾...
-
エクセルの標準偏差を用いた%計...
-
(1)「イギリス、日本、東京にお...
-
標準偏差から実際のばらつきを...
-
正規分布曲線のカーブについて
おすすめ情報