- ベストアンサー
統計学：サンプリング精度の応用についての疑問
- サンプリング精度の応用についての疑問を解決
- 推測したいグラウンド全体の黒の碁石の総数の95%信頼区間の計算方法について教えてください
- 統計学でサンプリング精度を適用する際の疑問について解説
- みんなの回答 （4）
- 専門家の回答
質問者が選んだベストアンサー
- ベストアンサー
その他の回答 (3)
> 「試行回数」は「総数」と同じと扱う理由はなんでしょうか？ 碁石をグラウンド内に一個ずつ入れていくことをイメージしてください。 碁石一個一個について、サンプリング率の確率でサンプリングエリア内に入れるか入れないかという試行を碁石の総数だけ行っているので、試行回数と総数は同じになります。
お礼
実際の試行回数ではないのですね。 いろいろありがとうございました。 長く悩んでいた問題が解決しました。
> やっぱり自力では不等式をnについて解くことはできませんでした。 |(x-na/A)/√{n(a/A)(1-a/A)}| <= z を二乗すると (x-na/A)^2/{n(a/A)(1-a/A)} <= z^2 (x-na/A)^2 <= z^2{n(a/A)(1-a/A)} あとはnについての二次不等式を解くだけです。 No.2でいきなりサンプリング率をpにしていましたね…… 修正しきれてませんでした。 > 「822278≦n≦1216135」は中央値が「1000000」ではないのですね。感覚として不思議です。 これは二項分布がp=0.5以外のとき非対称な分布だからでしょう。 > 数式上、p=100%では正確な区間がでないのは正規分布に近似したため なのでしょうか？ 確かに、近似のため正確な信頼区間ではありませんが、p=100%のときはx≦n≦xつまりn=xとなります。 計算間違いされてませんか？
お礼
お礼が遅れまして申し訳ございません。 色々ありがとうございました。 実は「計算間違いされてませんか？」と問われてExcelの式を確認した結果、修正して p=100%でも成立しました。 数日に渡ってありがとうございました。 最後に、ひとつ理解できなかった点があります 「「碁石総数 = 試行回数」なので、どちらもnです」。 「試行回数」は「総数」と同じと扱う理由はなんでしょうか？
面積A m^2のグラウンド中に1個の碁石があるとすると、面積a m^2のサンプリングエリア内に碁石が見つかる確率はa/Aとなるので、碁石がn個ならサンプリングエリア内に見つかる碁石の数xは試行回数n、成功確率a/Aの二項分布に従います。 xが正規分布で近似でき、標準正規分布の97.5%点をzとすると、95の確率で |(x-na/A)/√{n(a/A)(1-a/A)}| <= z となります。 これをnについて解けば、A m^2のグラウンド中に存在する碁石の数の95%信頼区間が得られます。
お礼
早速の回答ありがとうございます。スマートな解説を頂き感謝いたします。 「碁石総数がN個」「試行回数n」とすると |(x-Na/A)/√{n(a/A)(1-a/A)}| <= z でしょうか？数式の意味は良く理解できませんでしたが、z＝1.96として、Nで解いても解が得られませんでした。私の計算法が誤っているのかもしれません。 例えば、サンプリング率（a/A）0.01%（=1/10,000）で（x）100個カウントできた場合は、総数は1x10^6個と推定できる思うのですが、その際の95%信頼区間幅はどれくらいになるのでしょうか？ サンプリング率（a/A）とサンプリングエリア内のカウント数（x）のみから求める計算式となると思うのですが。
お礼
詳細な回答ありがとうございました。 やっぱり自力では不等式をnについて解くことはできませんでした。 (2x+z^2*(1-p)-√[{2x+z^2*(1-p)}^2-4*x^2])/(2p)≦n≦(2x+z^2*(1-p)+√[{2x+z^2*(1-p)}^2-4*x^2])/(2p) は何とか数値を代入してを確認しました。 「822278≦n≦1216135」は中央値が「1000000」ではないのですね。感覚として不思議です。また、サンプリング率（p）が0.001%でも10%でも、推定値に対する95%信頼区間の割合はあまり変動しないのですね。推定値に対する95%信頼区間の割合を±10%程度に抑えるには（p）に関係なく「400個」程度数えれば良いことが理解できました。数式上、p=100%では正確な区間がでないのは正規分布に近似したためなのでしょうか？ 対象数が多い場合に限定して考えたいと思います。 何度も質問にお付き合い頂き誠にありがとうございました。統計学を基礎から勉強したいと思います。