統計解析：CV%の合成（？）について

Question

ある分析法AのCV%がx、別の分析法BのCV%がy、と判っている時、それらの分析法で得られたデータを合成した値のCV%は、どう計算すれば良いでしょうか？

例1）
分析法Aで得られたデータをa、分析法Bで得られたデータをbとする時、「a/b」のCV%は？

例2）
分析法Aで得られたデータをa、分析法Bで得られたデータをbとする時、「a*b」のCV%は？

...よろしければ、こういう事が判りやすく書かれている実践的な内容の書籍についても、アドバイスを頂ければ助かります。

grothendieck · Accepted Answer

申し訳ありませんん。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))
}

シミュレーションでやるのがどんな場合でも間違いなく確実にできるとも言えます。

grothendieck · Answer

CVa/b = √(CVa^2 + 1)/(CVb^2 + 1) - 1)

は

　CVa/b = √( (CVa^2 + 1)/(CVb^2 + 1) - 1)

に訂正させて頂きます。

grothendieck · Answer

「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が１に比べて小さく、
かつ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

grothendieck · Answer

無理なことなんかない。確率変数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
でした方が良いと思います（私は質問したことはありません）

endlessriver · Answer

Cv=(s/Xa)*100
sは標準偏差、Xaは平均ですから無理と思います。

統計解析：CV%の合成（？）について

申し訳ありませんん。

CVa/b = √(CVa^2 + 1)/(CVb^2 + 1) - 1)

「CVab = √(CVa^2・CVb^2 + Cva^2 + CVb^2)」

無理なことなんかない。

Cv=(s/Xa)*100

似たような質問が見つかりました

関連するカテゴリからQ&Aを探す

デイリーランキングこのカテゴリの人気デイリーQ&Aランキング

マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング

　CVa/b = √(CVa^2 + 1)/(CVb^2 + 1) - 1)