No.5ベストアンサー
- 回答日時:
申し訳ありませんん。
No3,4の回答は誤りのようです。b = (1+ε)mb
より
1/b ≒ (1-ε)/mb、 1/b^2 ≒ (1-2ε+ε^2)/mb^2
としましたが、1/b^2をε^2のオーダーまで計算するには1/bもε^2のオーダーまでとる必要がありました。そこで
1/b ≒ (1-ε+ε^2)/mb、1/b^2 ≒ (1-2ε+3ε^2)/mb^2
に修正します。
σ^2[a/b] = E[(a/b - E[a/b])^2]
= E[a^2]・E[1/b^2] - (ma・E[1/b])^2
に
E[a^2] = (1 + CVa^2)ma^2
E[1/b^2] = (1 + 3 CVb^2)/mb^2
E[1/b] = (1 + CVb^2)/mb
を代入すると
σ^2[a/b] = (CVa^2 +CVb^2 + 3CVa^2・CVb^2)ma^2/mb^2
よって
CVa/b =√(CVa^2 +CVb^2 + 3CVa^2・CVb^2) /(1+CVb^2)
すると下のようにシミュレーションとの一致が良くなりました。(CVb<CVaという条件は必要ありません)
ma sda mb sdb 乱数 上の式
500 9 1000 3 0.0183 0.0182
500 15 1000 5 0.0305 0.0304
500 20 1000 5 0.0403 0.0403
1000 3 200 9 0.0453 0.0450
間違えやすい点を多々含んでおり、(E[1/b]=1/E[b]とかついやってしまいそう)これを考えるのは良い演習問題になると思います。上の計算に使ったRの関数を次に示します。
cvawb <- function(ma,sda,mb,sdb){
a <- rnorm(500000,mean=ma,sd=sda)
b <- rnorm(500000,mean=mb,sd=sdb)
ab <- a/b
cab <- sd(ab)/mean(ab)
ca <- sda/ma
cb <- sdb/mb
cv <- sqrt(ca^2 +cb^2 +3*ca^2*cb^2)/(cb^2 +1)
return(list(cab,cv))
}
シミュレーションでやるのがどんな場合でも間違いなく確実にできるとも言えます。
再度の詳細な御説明、誠に有り難う御座います。大変助かっております。
素人目には、「3CVa^2」の「3」が、そういう数字が出るのだなと面白く感じられました。
また「シミュレーションでやるのが...」は、なるほどなと思わされました。今は、どこにでも「パーソナル」コンピュータがあり、こういう計算やシミュレーションも個人で行えてしまう、大変良い時代なのですね。
締め切ってしまうと、重ねての御質問がしにくくなり、惜しい気がするのですが... ともあれ、誠に有り難う御座いました。また、どうぞ宜しく御願い致します。
No.4
- 回答日時:
CVa/b = √(CVa^2 + 1)/(CVb^2 + 1) - 1)
は
CVa/b = √( (CVa^2 + 1)/(CVb^2 + 1) - 1)
に訂正させて頂きます。
再び詳細に御回答頂き、大変助かります。有り難う御座いました。
素人ながら薄々思っていたのですが、「a/bは一般には複雑になると思います」との事、やはりそうなのか...と感じました。簡単な計算では求められないのですね。
この中の「かつCVb<CVaのときの近似式」という点についてですが、この近似式はCVb>CVaだとどこを変えなくてはいけないのでしょうか。そもそも、お示し頂いている計算の過程が良く判っていないのでお恥ずかしいのですが、この式(最終の公式や、計算の過程)の中でCVb<CVaの場合とCVb>CVaの場合では何が違うのか、知っておかなくてはと思いました。もしも御時間が頂けましたら、恐縮ですがお教え下さい。
No.3
- 回答日時:
「CVab = √(CVa^2・CVb^2 + Cva^2 + CVb^2)」
は変動係数を
Cv=(s/Xa) sは標準偏差、Xaは平均
としたときのabの変動係数です。CV%だと
「CVab%=√(CVa%^2・CVb%^2/10000 + Cva%^2 + CVb%^2)」
となります。
a/bは一般には複雑になると思います。簡単には下の回答に書いた様に乱数を使ってシミュレートして評価することができると思います。CVbが1に比べて小さく、
かつCVb<CVaのときの近似式を次ぎに導きます。E[b]=mbとして確率変数bがとる値を
b = (1+ε)mb
とすると
1/b ≒ (1-ε)/mb
この近似の下でE[1/b]=1/mbであり、
E[a/b] = ma/mb
σ^2[a/b] = E[(a/b - ma/mb)^2]
= E[a^2]/E[b^2] - (ma/mb)^2
(正確には E[1/b]=1/mb でない)などより
CVa/b = √(CVa^2 + 1)/(CVb^2 + 1) - 1)
前と同様に乱数でシミュレートした結果とこの式の比較例は
ma sda mb sdb 乱数 上の式
500 9 1000 3 0.0182 0.0177
500 15 1000 5 0.0305 0.0296
500 20 1000 5 0.0403 0.0397
No.2
- 回答日時:
無理なことなんかない。
確率変数aとbが独立でE[a]=ma, E[b]=mb とします。このときE[ab] = ma・mb
σ^2[ab] = E[(a・b - ma・mb)^2]
= E[a^2]E[b^2] - (ma・mb)^2
これから、aの変動係数をCVa、bの変動係数をCVb、abの変動係数をCVabとすると
CVab = √(CVa^2・CVb^2 + Cva^2 + CVb^2)
aとbが正規分布に従うとして、平均と標準偏差を与えた時、乱数で計算したabの変動係数と上の式で計算したCVabを比較する関数を統計ソフトR
http://cse.naro.affrc.go.jp/takezawa/r-tips/r.html
で作りました。
cvab <- function(ma,sda,mb,sdb){
a <- rnorm(500000,mean=ma,sd=sda)
b <- rnorm(500000,mean=mb,sd=sdb)
ab <- a*b
cab <- sd(ab)/mean(ab)
ca <- sda/ma
cb <- sdb/mb
cv <- sqrt(ca^2*cb^2 + ca^2 + cb^2)
return(list(cab,cv))
}
aの平均ma、標準偏差sda、bの平均mb、標準偏差sdbを与えたときの計算例を示します。
ma sda mb sdb 乱数 上の式
10 1 20 2 0.1414726 0.1417745
10 1 20 4 0.2249844 0.2244994
10 3 20 2 0.3179832 0.3176476
10 3 20 4 0.3651448 0.3655133
なお統計に関する質問はここでするより
http://aoki2.si.gunma-u.ac.jp/lecture/mb.html
でした方が良いと思います(私は質問したことはありません)
御回答、大変有り難う御座いました。
済みません、素人なものでまだ理解が不足しているのですが、ここで
「CVab = √(CVa^2・CVb^2 + Cva^2 + CVb^2)」
としてお示し頂きましたのは「a*b」のCV%の計算式であるという理解でよろしいでしょうか。...まだ、これを応用出来るまでの力が無いもので、「a/b」のCV%ならどういう計算式になるか、併せて御教示頂けましたら大変助かります。
統計ソフトRや、「統計学関連のメッセージボード」の件、気になってはいながらそれぞれに敷居が高かったのですが、今後活用して行きたいと思います。
こういう分野は、困った時に助けを求められる相手が周囲にいませんので、grothendieckさんにこうして教えて頂ける様な事があると、ウェブというのは有り難いものだなと思います。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
(a+b)(a2-ab+b2)=a3+b3の途中式
-
分数に未知数を含む等式の変形...
-
「因数定理」は、いつ習います...
-
数学
-
高校数学です。 (-a2乗bx3乗)...
-
ab-b-ac+c この式を因数分解し...
-
ルートの公式と複素数について
-
必要条件・十分条件・必要十分条件
-
【数学】2√3の整数部分をa, 少...
-
a³+b³の解き方
-
三重根号を簡単にする問題です...
-
6の倍数であることを証明
-
因数分解 松山大・秋田大 入...
-
数学が得意な方
-
x^2-16x-512の因数分解をすぐに...
-
数学の展開
-
高校1年レベルの問題 因数分解
-
9b-9-3ab+a² の因数分解の仕...
-
因数分解の基本
-
a^4+3a^2-2ab+4-b^2の因数分解...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
異なる2つの無理数の積について
-
(a+b)(a2-ab+b2)=a3+b3の途中式
-
lim{(a^x+b^x)/2}^1/x x→0 (a...
-
三重根号を簡単にする問題です...
-
6の倍数であることを証明
-
「因数定理」は、いつ習います...
-
a^n+b^nの因数分解の仕方
-
高校生です! 数学の質問です。...
-
【数学】2√3の整数部分をa, 少...
-
ベクトルの問題
-
他の式を利用した因数分解 x^3+...
-
群の乗積表の作り方は?
-
異なる自然数a,bがあり, a+bを7...
-
立方完成,N乗完成は存在するの?
-
ab>a+bは常に成り立つでしょうか?
-
不等式の問題
-
楕円の変数変換
-
ax + by (a,bは自然数で互いに...
-
この問題の解き方教えてくださ...
-
2離れた奇数が互いに素なこと...
おすすめ情報