A 回答 (3件)
- 最新から表示
- 回答順に表示
No.3
- 回答日時:
※この回答は、“締め切られた質問への回答追加”として、2021/11/05 18:29 に回答者の方よりご依頼をいただき、教えて!gooによって代理投稿されたものです。
---
「正規分布で近似する」という話がどこから出てくるんだか、ちょっと不思議です。
「確率θの選択をN回試行してr回実現した、という場合のθの信頼区間の上限」とは、Nとθを母数とする二項分布を B[N,θ](k)とするとき、閾値pについて
Σ{k=0〜r}B[N,θ](k) ≦ p である最小のθ
ということであり、一方、自由度(n1,n2)のF分布をG[n1,n2](x)とすると、近似なしに
Σ{k=0~r}B[N,θ](k) = ∫{x=X〜∞} G[2(r+1),2(N-r)](x) dx
ただし、X = (r+1)θ/((N-r)(1-θ))
です。いわゆる「F分布表」は、この積分の値をlook upするとXが分かるように作ってある。だから、積分 ≦ p になるXを調べて
θ = (N-r)X / ((r+1) + X(N-r))
を計算するのが、当然のやり方。下限も同様です。
No.2
- 回答日時:
二項分布の信頼区間の推定方法はいくつもあり、Rで試すことができます。
やってみたところ、それぞれ微妙に違いますね。
> binom.confint(64, 752, conf.level=0.9, methods="all")
method x n mean lower upper
1 agresti-coull 64 752 0.08510638 0.06975483 0.1034326
2 asymptotic 64 752 0.08510638 0.06836911 0.1018437
3 bayes 64 752 0.08565737 0.06885151 0.1022490
4 cloglog 64 752 0.08510638 0.06935077 0.1028139
5 exact 64 752 0.08510638 0.06893687 0.1037474
6 logit 64 752 0.08510638 0.06979368 0.1034052
7 probit 64 752 0.08510638 0.06957288 0.1031040
8 profile 64 752 0.08510638 0.06936330 0.1028375
9 lrt 64 752 0.08510638 0.06936147 0.1028371
10 prop.test 64 752 0.08510638 0.06661902 0.1079505
11 wilson 64 752 0.08510638 0.06982041 0.1033671
これらの方法も含めて、勉強しつつプログラムをする予定です。
No.1
- 回答日時:
なかなか回答が付きませんが、私から。
難問ですよ。これは。
まず、これだけ項数が多いと(752)、コンビネーションや二項確率の計算ができません。そのため、95パーセンタイル値とか、99パーセンタイル値が計算できません。
信頼区間を議論するには、このようなすそ野の一致性を議論する必要があります。
普通は二項分布の分布形状が対称形に近ければ、平均値まわりの値は正規分布との近似性が高いので、npの値がまあまあ大きければ、という話になるのですが・・・。
まさか、平均値まわりの話をしているのではないと思います。
実は正規分布近似は、すそ野の近似性は極めて悪いんです。理由は「確率」は非負なのに正規分布の範囲はー∞から∞ですから。
そのため、裾野の議論をしたいときは、「ロジット変換」か「逆正弦変換」をしてすそ野の近似性を得ます。QCをやっている人にとっては常識の話です。QC検定に出題されますからね。
方針は、上記を踏まえて、
①まず、二項分布では計算ができません。二項分布の項数が大きい時はポアソン分布で近似せよ、といわれているので、ポアソンでnを変えて99パーセンタイルを調べる。pは64/752に固定する。
②常套手段である「ロジット変換」「逆正弦変換」「正規分布近似」のそれらと比較する。たぶん、正規分布近似がだんだん漸近してくるはず。期待する精度は指定されていませんが、±1%くらいですよね。その時のnを求める。
これはプログラムを書かないとできませんね。今度の土日にやってみようと思いますので、閉じずに待っていて下さい。
また、どの書籍か知りませんが、正解があるなら示して下さい。もしnpが・・・と書いてあったら、そんな本は即刻粗大ごみに出すべきです。
この回答へのお礼
お礼日時:2021/10/28 23:51
ご回答をしていただきありがとうございます。
大学の統計問題でまだ回答がわかりません。多分来週になるかと思います。
もし週末に解けたらご回答をよろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 統計学 こんな問題を使って教育するのは、文科省の方針ですか。 3 2022/06/17 09:14
- 統計学 統計学の問題です。教えてください(_ _) 数万人の有権者がいる選挙区で, 無作為に400人の標本を 2 2023/02/03 15:27
- YouTube Googleアカウント1個で複数のYouTubeチャンネルを作った場合、登録者数などは通算される? 4 2022/05/16 14:22
- 固定資産税・不動産取得税 登録免許税の計算について 1 2023/02/03 21:40
- 統計学 t検定について教えてください 2 2023/02/23 16:35
- 消費者問題・詐欺 お金を取り返すことは可能でしょうか? 4 2023/01/07 13:17
- いじめ・人間関係 精神障害者で生活保護の自宅療養者が、人として自分を見てもらうアイデアはありますか? 例えば、安い居酒 7 2022/08/23 07:32
- 統計学 信頼区間についての質問です。 6 2023/06/25 17:34
- 統計学 以下の問題が分からないので計算式を教えてください ある企業が製造している電球の寿命の母平均と,母標準 3 2023/01/14 00:43
- 統計学 統計学の質問【帰無仮説】 高校の新学習指導要領では、統計的仮説検定の基本的な考え方が必修単元となった 5 2023/05/23 21:00
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
EXCELでページに飾り枠を...
-
カンプって何ですか
-
紙コップの版下作り(Illustrat...
-
仕事が上手くいかない・上手く...
-
本を飾れる額縁ってありますで...
-
商品化してくれる企業はないか
-
セルフネイルをしています。イ...
-
かっこいいサイン(署名)をデ...
-
パワーポイントでデザインを削...
-
サザビーとシナンジュは同じ人...
-
商品デザインのマネ
-
会社のロゴデザイン
-
皆さんはぁぁぁ、 MacBook Pro...
-
デザインを盗むのは罪?
-
座り心地が悪い座椅子の改善法...
-
このスウェットは外で着るのは...
-
志望理由書
-
情報デザイン
-
28歳女、今からデザインの学校...
-
友達に誕生日プレゼントでお箸...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
EXCELでページに飾り枠を...
-
昔のグラニフに似てる服屋さん...
-
チラシやフライヤー、バナーを...
-
英文の改行位置について
-
紙コップの版下作り(Illustrat...
-
黄色い英語長袖tシャツ探してま...
-
女の子水着になる目的?
-
あなたはデザイン重視?機能性...
-
デファクトスタンダードとドミ...
-
OPPOってiPhoneに似てませんか...
-
カンプって何ですか
-
デザインにどれくらいの時間が...
-
WEBライターやるくらいならWEB...
-
セルフネイルをしています。イ...
-
発表のスライドの作り方のコツ...
-
自分がデザインした洋服を売り...
-
lone star state
-
本を飾れる額縁ってありますで...
-
FAXで透かし付文書を送ると、真...
-
おいしい牛乳 ロゴ 元になっ...
おすすめ情報