![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
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で質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・「I love you」 をかっこよく翻訳してみてください
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・昔のあなたへのアドバイス
- ・かっこよく答えてください!!
- ・あなたが好きな本屋さんを教えてください
- ・スタッフと宿泊客が全員斜め上を行くホテルのレビュー
- ・【大喜利】【投稿~8/27】 こんなガソリンスタンド二度と来るか!なぜそう思った?
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・【お題】動物のキャッチフレーズ
- ・【お題】甲子園での思い出の残し方
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・「それ、メッセージ花火でわざわざ伝えること?」
- ・自分用のお土産
- ・人生で一番お金がなかったとき
- ・一番好きなみそ汁の具材は?
- ・泣きながら食べたご飯の思い出
- ・ちょっと先の未来クイズ第1問
- ・ゴリラ向け動画サイト「ウホウホ動画」にありがちなこと
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
三重根号を簡単にする問題です...
-
「因数定理」は、いつ習います...
-
ベクトルの問題
-
(a+b)(a2-ab+b2)=a3+b3の途中式
-
abx²-(a²+b²)x+abを因数分解す...
-
異なる自然数a,bがあり, a+bを7...
-
(3)の三番目の式で、(1-cosx)2...
-
方程式と不等式の問題で計算が...
-
異なる2つの無理数の積について
-
計算が苦手です。 4a四乗−25a二...
-
統計解析:CV%の合成(?)につ...
-
数学の問題をといてください。
-
a^3+b^3=(a+b)(a^2-ab+b^2)とな...
-
数学
-
6の倍数であることを証明
-
a^n+b^nの因数分解の仕方
-
平方すると、-18i になる複素数...
-
不等式の証明が分からないです ...
-
lim{(a^x+b^x)/2}^1/x x→0 (a...
-
|a|-|b|≦|a-b| 等号成立
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
6の倍数であることを証明
-
(a+b)(a2-ab+b2)=a3+b3の途中式
-
a^3+b^3=(a+b)(a^2-ab+b^2)とな...
-
三重根号を簡単にする問題です...
-
「因数定理」は、いつ習います...
-
数学1についてです。 a=bはa^2+...
-
【数学】2√3の整数部分をa, 少...
-
異なる2つの無理数の積について
-
アーベル群
-
高二 数学 指数の拡張 写真の問...
-
lim{(a^x+b^x)/2}^1/x x→0 (a...
-
群の乗積表の作り方は?
-
2離れた奇数が互いに素なこと...
-
a^n+b^nの因数分解の仕方
-
この不定積分が解けません
-
他の式を利用した因数分解 x^3+...
-
ab>a+bは常に成り立つでしょうか?
-
数学
-
高一数学です。なるべくお早め...
-
行列式の因数分解
おすすめ情報