ようこそ ゲスト さん、新規登録(無料)して気になる疑問を解決しませんか?

質問

QNo.4150487 超幾何分布
質問者:nattin_111 超幾何分布と二項定理について教えてください
(1+x)^(m+n)
=(1+x)^m*(1+x)^n
=Σ{i=0〜m}mCi*x^i*Σ{j=0〜n}nCj*x^j
=Σ{i=0〜m}*Σ{j=0〜n}nCj*mCi*nCj*x^(i+j)
 
(1+x)^(m+n)=Σ{m=0〜m+n}(m+n)Cm*x^m
に持って行きたいのですが、問題の3つ目の=から先が理解できません。
シグマ二つをどう処理すればよいのでしょうか?
アドバイスいただけませんか。
よろしくお願いいたします
困り度:
  • すぐに回答を!
質問投稿日時:
08/07/04 12:57
この質問に対する回答は締め切られました。

回答良回答20pt

ANo.3 どういう順でお答えするのが良いでしょう・・・。

まず最初に、
> Σ(k=0〜p) (m C k) (n C (p-k)) = (m+n) C p
> というのは公式なんですか?当たり前のように使用していいものなのでしょうか?
私は当たり前として良く使います。説明するとなると、#2 でお答えしたように説明されることが多いのではないかと思うのですが、逐一説明する必要は無いでしょう。

> Am = Σ(i=0〜m) (m C i) (n C m-i) のΣ(i=0〜m) が残るのはどうしてでしょうか。
> この問題においては、i=kのときをさしていると思うので、
> Am=(mCk)(nCm-k)かと思ってしまったのですが・・
まず、(1 + x)^(m+n) において x^m の係数を Am とすると、
Am = (n + m) C m ・・・(1)
は良いですね。
これと、(1 + x)^m (1 + x)^n の x^m の係数とは一致するはずですよね。
(1 + x)^m (1 + x)^n
= (mC0 + mC1 x + ... + mCm-1 x^(m-1) + mCm x^m) (nC0 + nC1 x + ... + nCn-1 x^(n-1) + nCn x^n)
この計算で得られる x^m の係数は、左の括弧の x^k の項と 右の括弧の x^(m-k) の項の積を k=0 〜 m で加算したものになることは分かりますか?
つまり、左の括弧の x^0 の項と右の括弧の x^m の項の積と、左の括弧の x の項と右の括弧の x^(m-1) の項の積と、・・・左の括弧の x^k の項と右の括弧の x^(m-k) の項の積と、・・・左の括弧の x^m の項と右の括弧の x^0 の項の積とを全て加算したのが、(1 + x)^m (1 + x)^n の x^m の項になるのです。そうして得られるx^m の係数が、Σ(k=0〜m) { (m C k) (n C (m-k)) } です。
これが、(1) で求めた x^m の係数 Am と等しいのだから、
Σ (i=0〜m) {(m C i) (n C m-i)} = (m+n) C m

これがどうしても分からなければ、具体的に数字を入れて自分で手を動かしてみる。
m=2, n=3 とすると、(1+x)^5 の x^2 の項の係数 5C2 = 10 。
一方で、
(1 + x)^5 = (1+x)^2 (1+x)^3
= (2C0 + 2C1 x + 2C2 x^2) (3C0 + 3C1 x + 3C2 x^2 + 3C3 x^3)
これを展開して得られる x^2 の項は、
左括弧の x^0 の項(2C0) と 右括弧の x^2 の項( 3C2 x^2) の積 = 2C0 3C2 x^2
左括弧の x^1 の項(2C1 x) と 右括弧の x^1 の項( 3C1 x) の積 = 2C1 3C1 x^2
左括弧の x^2 の項(2C2 x^2) と 右括弧の x^0 の項( 3C0) の積 = 2C2 3C0 x^2
のこれら全ての和だから、x^2 の係数は Σ(k=0〜2) { 2Ck * 3C(2-k) }
です。これが (1 + x)^5 の x^2 の係数 5C2 なんだから、
Σ(k=0〜2) { 2Ck * 3C(2-k) } = 5C2
ってことです。

> Σ(k=0〜m)(mCk)(nCm-k)=m+nCmとなり、両辺をm+nCmで割るのだと思いますが、
> 割るときにm+nCmはシグマのなかに入れてしまっていいのでしょうか?
Σ(k=0〜m) P(X = k)
= Σ(k=0〜m) [ {mCk*nC(m-k)} / (m+n)Cm ]
= (1 / (m+n)Cm) Σ(k=0〜m) {mCk * nC(m-k)}  (∵ (m+n)Cm は k とは無関係に定数 )
= (1 / (m+n)Cm) ((m+n) C m)   (∵ Σ(k=0〜m) {mCk * nC(m-k)} = (m + n) C m )
= 1

ラストアンサーでいい?
回答者:kumipapa
種類:アドバイス
どんな人:一般人
自信:参考意見
回答日時:
08/07/06 11:46
この回答へのお礼理解できました!!
事細かに丁寧にお教えいただいて本当に助かりました。
もう一度自分で考えて、モノにできるようにしてみます。
本当にありがとうございました。

回答

ANo.2 問題は分かりました。

> (1+x)^(m+n)
> =(1+x)^m*(1+x)^n
> =Σ{i=0〜m}mCi*x^i*Σ{j=0〜n}nCj*x^j
> =Σ{i=0〜m}*Σ{j=0〜n}nCj*mCi*nCj*x^(i+j)
ここでなぜ n C j が 2 度出てくる(n C j を2乗している)のですか?

> (1+x)^(m+n)=Σ{m=0〜m+n}(m+n)Cm*x^m
> で分母の部分を作って、
細かいようですが、相変わらず、右辺で m が混乱しています。
(1+x)^(m+n)=Σ{k=0〜m+n}(m+n)Ck*x^k
ですね。

さて本題。
> P(X=k)={mCk*nC(m-k)}/(m+n)Cm であるとき、
確かに超幾何分布です。
m個の不良品、n個の良品、合わせてm+n個の製品の中からm個を取り出したとき、そのm個の中にk個の不良品が含まれている確率です。
ここで、
Σ(k=0〜m)P(X=k) = Σ(k=0〜m) { (m C k) (n C (m-k) } / ((m+n) C m) = 1
を示せというわけですね。
まず
Σ(k=0〜p) (m C k) (n C (p-k)) = (m+n) C p
というのは、今まで使ったことはないですか?
m+n 個の物から p 個の物を取り出す場合の数は (m+n) C p であるが、これは、m+n 個の物を m 個のグループ A と n 個のグループ B の2つにわけ、
グループ A から 0 個と グループ B から p 個を取り出して組み合わせる場合 m C 0 × n C p 通り
グループ A から 1 個と グループ B から p-1 個取り出して組み合わせる場合 m C 1 × n C p-1 通り
グループ A から 2 個と グループ B から p-2 個取り出して組み合わせる場合 m C 2 × n C p-2 通り
・・・
グループ A から p 個と グループ B から 0 個取り出して組み合わせる場合 m C p × n C 0 通り
の場合の数を全て足した数と等しい。
∴ Σ(k=0〜p) (m C k) (n C (p-k)) = (m+n) C p
この問題の場合は、p=m だから
Σ(k=0〜m) (m C k) (n C (m-k)) = (m+n) C m
∴ Σ(k=0〜m) P(X=k) = Σ(k=0〜m) { (m C k) (n C (m-k) } / ((m+n) C m) = 1

数式で示すならば、質問者さん流を踏襲してみますが、(1+x)^(m+n) の一般項をごちゃごちゃいじっても仕方がない。問題なのは、x^m の項の係数だけなんだから計算のターゲットはそこに絞りましょう。

(1 + x)^(m + n) = Σ(k=0〜m+n) ((m+n) C k) x^k
より、(1 + x)^(m+n) の x^m の係数は (m+n) C m ・・・(1)
一方で、
(1 + x)^(m+n)
= (1 + x)^m (1 + x)^n
= ( Σ(i=0〜m) m C i x^i ) ( Σ( j =0〜n) n C j x^j)
ここで x^m の係数を Am としてそれのみを計算すれば
(上式の左側括弧の x^i の項の係数と 右側括弧の x^(m-i) の項の係数との積を i=0〜mで加算すればよいから)
Am = Σ(i=0〜m) (m C i) (n C m-i) ・・・ (2)
(1), (2) より
Σ(i=0〜m) (m C i) (n C m-i) = (m+n) C m
∴ Σ(k=0〜m) P(X=k) = Σ(k=0〜m) { (m C k) (n C (m-k) } / ((m+n) C m) = 1

なお、
if p≧q, p C q = p!/(q! (p-q)!)
otherwise, p C q = 0
としています。念のため。
回答者:kumipapa
種類:アドバイス
どんな人:一般人
自信:参考意見
回答日時:
08/07/05 01:46
この回答への補足質問者です。
先ほどお礼で、質問をさせていただいたのですが、
すみませんが、もう1点あるので補足に書かせていただきました。

(1)(2)より、
Σ(k=0〜m)(mCk)(nCm-k)=m+nCmとなり、両辺をm+nCmで割るのだと思いますが、
割るときにm+nCmはシグマのなかに入れてしまっていいのでしょうか?
というかこのシグマはどこに対するシグマかがわかりません・・

申し訳ありません・・お教えください。
この回答へのお礼> (1+x)^(m+n)
> =(1+x)^m*(1+x)^n
> =Σ{i=0〜m}mCi*x^i*Σ{j=0〜n}nCj*x^j
> =Σ{i=0〜m}*Σ{j=0〜n}nCj*mCi*nCj*x^(i+j)
ここでなぜ n C j が 2 度出てくる(n C j を2乗している)のですか?
打ち間違えです・・すみません。

ご丁寧にありがとうございます!
ひとつ質問なのですが、
Am = Σ(i=0〜m) (m C i) (n C m-i)
のΣ(i=0〜m) が残るのはどうしてでしょうか。
この問題においては、i=kのときをさしていると思うので、
Am=(mCk)(nCm-k)かと思ってしまったのですが・・
基本的なことかもしれませんがご教授ください。

あと、
Σ(k=0〜p) (m C k) (n C (p-k)) = (m+n) C p
というのは公式なんですか?当たり前のように使用していいものなのでしょうか?

たびたびで申し訳ありませんが、よろしくお願いいたします。

回答

ANo.1 > (1+x)^(m+n)=Σ{m=0〜m+n}(m+n)Cm*x^m に持って行きたいのですが・・・
式をよーく見てください。右辺の m が混乱しています。{ m=0 〜 m+n } ってことは、m は変数? 定数? 何者ですか?

もともとが、普通に二項定理で
(1 + x)^(m+n) = Σ[i=0,m+n] ((m+n) C i) x^i
ですよね。何を主張したくて、(1 + x)^(m+n) を考えたのか、さらにそれをなぜ (1 + x)^m (1 + x)^n と分けて展開しはじめる必要があったのかが良く分かりません。最終目標が何なのか、補足していただけませんか?超幾何分布と関係があるなら、超幾何分布の何を示したいのでしょうか?
回答者:kumipapa
種類:アドバイス
どんな人:一般人
自信:参考意見
回答日時:
08/07/04 13:55
この回答への補足P(X=k)={mCk*nC(m-k)}/(m+n)Cm であるとき、
Σ{k=0〜m}P(X=k)=1を示せ。

なのですが、
(1+x)^(m+n)
=(1+x)^m*(1+x)^n
=Σ{i=0〜m}mCi*x^i*Σ{j=0〜n}nCj*x^j
=Σ{i=0〜m}*Σ{j=0〜n}nCj*mCi*nCj*x^(i+j)
で分子の部分を、

(1+x)^(m+n)=Σ{m=0〜m+n}(m+n)Cm*x^m
で分母の部分を作って、最後に割って証明しようとしたのです。

方針が立たないときに、この問題は超幾何分布だと聞いたので、そのまま質問してしまいました。
ネット等で調べて考えてみたのですが、
・・・
=Σ{i=0〜m}*Σ{j=0〜n}nCj*mCi*nCj*x^(i+j)
=ΣmCi*nCj*x^(i+j)
=Σ{l=0〜(n+m)}x^lΣ{i+j=l}mCi*nCjとなっており、
シグマ2つの処理の仕方が理解できなかったため
質問させていただいたわけです。

わかりにくくてすみません。
再度アドバイスいただけると助かります。
この回答へのお礼この回答にお礼をつける(質問者のみ)