白黒の玉をランダムに並べたとき、白い玉がN個続く確率を、白の玉の数の割合pで書き表せますでしょうか?
そして、さらに、Nの平均値を求めたいのですがどうすればよいのでしょうか。
ぱっと思いついたのは、「N長鎖の前後が黒」、と考えると、確率は、
R(N,p)=(1-p)×(p^N)×(1-p) で、
平均値は、<N>=Σ_{N=0~∞} R(N,p)×N
という答えなのですが、実際に数値計算してみると、<N>がpに依存せず、どうもおかしいような気がしてきました。どこが間違っているのか、ご指摘いただければありがたいです。
A 回答 (7件)
- 最新から表示
- 回答順に表示
No.7
- 回答日時:
再びANo.5のコメントについてです。
ANo.6では「それは「分母」が違うからでしょうね。」と書いたけど、いや、むしろ違うのは分子の方じゃないかという気がしてきました。念のためにexcelでちょこっと実験してみましたら、ANo.5の結果と良く整合しています。ANo.5は間違っていないようだ。なのに、
> p=0.9801でm≡<N>≒49
> p=0.9025でm≡<N>≒9
> p=0.7569でm≡<N>≒3
ANo.5 だとこれらは50.3、10.3、4.1になる。ANo.5より見事に丁度1ずつ小さい値が並んでいます。
ということはですよ、ご質問とANo.5とでは「鎖の長さ」の定義が1だけ違っているんでしょうね。ANo.5では○の個数を数えているけれど、
●○●を長さ0
●○○●を長さ1
●○○○●を長さ2
という風に、○と○のスキマの数を長さとすれば話が合う。
しかしです。いやぁ、定義はともあれ「長さ0の鎖」を「鎖」と見るのはフツーじゃないでしょ。やっぱり変ですよねえ。
でも。でもでもしかしです。「長さ0だって鎖だ」と考えて分母に繰り入れるからこそ計算が合うんです。(もし「長さ0は鎖ではない」と言って勘定に入れなければ、答が再びANo.5と同じになってしまうのは自明ですから。)
No.6
- 回答日時:
ANo.5のコメントについてです。
> 結構ずれています。
それは「分母」が違うからでしょうね。つまり何に対する頻度を測っているかが違う。だとすると、まさにご質問の曖昧さがそのまま、数値の違いとして現れたわけで。
ご覧の文献(1次元イジング模型の臨界状態の話とか?)にあるR(N,p)はひょっとすると「単位長さ当たりにN長鎖が平均何個あるか」の期待値を「単位長さ当たりに鎖が平均何個あるか」の期待値で割ったものかも知れません。そうだとすると、Rは先の計算より小さく出ることになる。とは言っても、「単位長さ」がN以下ならN長鎖がそこに含まれる確率は0に決まっているから、これじゃ確率をきちんと定義したことになっていない。この方針で行くなら、なんらかの意味で極限(「単位長さ」Lを、L→∞としたときのN長鎖の出現頻度の期待値)を考える必要があるでしょう。
というわけで、まずは文献の数値が一体何なのか。すなわち、どう定義された(あるいは、どう測られる)数値なのか、という事をご確認なさって、再三指摘されたところの「曖昧さ」を解消してからお考えになってはいかがでしょうか。
この回答への補足
ランダムネスを確率pで入れていった時の連鎖長の話なので、元々の鎖の長さは∞のはずです。
文献の数値は、「簡単なestimationから求めると~」とあるだけでそれを実験結果と比べる話が主題なので、詳しい説明やrefは全くないのです。ですから、面倒な仮定をおいたり、複雑な計算をやっているとは思えないのです。
No.5
- 回答日時:
ランダムな無限列なら、Nが幾つだろうがN長鎖は必ず現れる。
つまりN長鎖の出現確率は1です。…という話では、どうもないらしい。だとすると、このご質問では、問うべき確率が定義できていないのです。(説明不足なのではなく、確率の概念についての理解不足でしょう。)
つまり、何に対するN長鎖の出現比率を指しているのかが不明。「全事象」に相当するものは何なのか。分母が決まっていないんですよ。
しかし、鎖の平均長を計算しようとしているのだから、「全事象」とはおそらく「全ての鎖」だろうと推察できます。つまり、
列を「n個目の要素は確率pで○、確率(1-p)で●」という規則に従ってランダムに作るものとする。この無限列に含まれる全ての鎖に占める、N長鎖の出現比率の期待値R(N,p)はいくらか。
これなら確率が考えられます。
まずは鎖だけに着目する。鎖とは「●の次に○が現れ、いくつか○が続いて、●で終わる」という部分列ですね。「●の次に○が現れる」までは、どんな鎖も同じ。
あとはこの後ろに○が幾つ続くかですが、確率(1-p)で●が来るから、その鎖の長さが1である確率は(1-p)。以下同様で、従って、
R(N,p)=(p^(N-1))(1-p)
(これは、ご質問にある式から(1-p)pを取り除いたものになってる。なぜ取り除いてあるかというと、鎖以外のものは最初から勘定に入れてないからです。)
で、Nの平均値mは
m=ΣNR(N,p) (ΣはN∈{1~∞})
=((1-p)/p)ΣN(p^N)
=1/(1-p)
この回答への補足
ご指摘ありがとうございます。おっしゃるとおり「相対的出現頻度」の意味で、確率という単語を連発してしまいました。すみません。
それで、結果なのですが、今読んでいる原稿に、
p=0.9801でm≡<N>≒49
p=0.9025でm≡<N>≒9
p=0.7569でm≡<N>≒3
と言う評価値が書かれていて、stomachamanさんの結果(50.3、10.3、4.1)とは微妙と言うか、結構ずれています。私が読んでいる原稿の評価値の方が、やはり間違いなのでしょうかな、、、。
ともかくありがとうございます。
No.4
- 回答日時:
なるほど。
それでしたら、細かな話はおいといて、黒 + N個の白 + 黒 と並ぶ確率は (1 - p) p^N (1 - p) となりますよね。しかし、連鎖長の平均値を
<N> = Σ[N=0,∞] N (1 - p) p^N (1 - p)
として計算したときに、これが p に依存しないというのはおかしいですね。
S(n) = Σ[N=0,n] N (1 - p)^2 p^N
とおいて
<N> = lim[n→∞] S(n)
S(n) = (1 - p)^2 Σ[N=0,n] N p^N = p - (n+1)p^(n+1) + n p^(n+2)
0 < p < 1
より
<N> = lim[n→∞] S(n) = p
という至って至ってシンプルな結論になりそうに思います。確認してみて頂けますか?ただし、収束していくに、n はかなりでかくしないと駄目みたいですが。
なお、興味をお持ちの話は、ベルヌーイ過程(離散マルコフ連鎖)の連鎖長の分布関数の話しと考えてよいと思います。もしまだ調べておられなければ、そのあたりを一度調べてみたらいかがでしょう。私はすっかり忘れてしまいましたが。
No.3
- 回答日時:
#1さんの「重要事項の欠落」という指摘に対して「説明不足の点はない」と言われますが、私は#1さんに1票。
数値計算がおできになるなら「モンテカルロ法」のプログラムを書いてみてはどうでしょうか。流れ図だけでけっこうです。重要事項の欠落に気がつくと思いますよ。
この回答への補足
var i,iw,isw:integer;
array a[1..1000] of integer;
const p=0.5;
(* definition: random takes value between 0 and 1.
definition: if random<p then white appears *)
iw:=0;
isw:=0;
for i:=1 to 1000 do
begin
if (isw=0)then
if (random<p)then
begin iw:=1; isw:=1 end
else
if (random<p)then
begin iw:=iw+1 end
else a[iw]:=a[iw]+1; iw:=0; isw:=0 end;
end;
for i:=1 to 1000 do write(i,' ',a[i]);
end.
すみませんが、どこが説明不足でしょうか?
pを固定値と宣言した段階で、玉の総数→∞の極限を
考えて欲しいなぁ、と思ったのですが、、、、。
ともかく、困ってるんでアイディアだけでもありませんか?
No.2
- 回答日時:
#1 に補足です。
玉が有限個とすると、白が N 個連続する確率を p^N とは書けないのだけど、例えば、白が出現する確率は常に p で、黒が出現する確率は 1 - p ということであれば、なんとなくそれらしくなりますね。
このように白、黒が出現する試行を M 回繰り返して、そのうち N 回連続で白、その他は全て黒となる確率ということならば、その確率は、M≧ N ≧ 1 で (M - N + 1) p^N (1 - p)^(M - N) ですね。
この回答への補足
いえ、そうではないのです。他は全部黒である必要はありません。
例えば、M=22個とすると、
○●●○○○●●●○○●○●●●○○●●●○
だと、白玉は、N=1が三回、N=2が二回、N=3が一回
となります。
これをM=∞に持っていったときに、Nの出現確率を求めたいのです。
Nの期待値の意味も通りますよね。
何とかなりますでしょうか。
(なお、もともと、物理のスピン系の問題に出てきた話です、、、。)
No.1
- 回答日時:
> 白い玉がN個続く確率を、白の玉の数の割合pで書き表せますでしょうか
白玉、黒玉の数が有限個ならば、pで単純に表すことはできません。
無限個ならば、白玉の比率が p というのはなんだかおかしな気がする。
> ぱっと思いついたのは、「N長鎖の前後が黒」、と考えると、確率は、
R(N,p)=(1-p)×(p^N)×(1-p)
残念ながら単なる勘違いでしょう。p^N として良い根拠が不明。また、前後を 1 - p で挟んだのは白の長鎖の前後に1個づつ黒が来ると言う意味ですか?合計N+2個の石を並べるとはどこにも書かれていません。合計N+2個であっても、(有限個の玉を並べるなら)間違い。直感と言えども、ずさん過ぎます。
この問題をまじめに解くにあたっては、問題の設定を細かくきちんとして頂く必要があります。白玉は何個あるのか、黒玉は何個あるのか分からなければこの問題は解けません。また、Nの定義が不明。Nが白玉の全数とするなら、「Nの期待値を求めたい」が意味不明。また、Nが白玉の全数でないなら、「N個連続する」の定義をきちんとしてもらわないと意味不明。例えば、白玉が10個あるとき、白玉が 2個、3個、3個、2個と4つの塊に分かれることもあるわけで、この場合、「N個連続する」をどう定義するのか。それとも、N個連続してあとはバラバラになる確率を求めたいということ?いずれにせよ、黒玉の数も重要です。有る一定数以上の黒玉が無ければ、白玉はバラバラになれない。
いずれにせよ、この問題を考えるにあたり重要な事項がすべて欠落しております。そこに気づくことからスタートしないと駄目かな・・・。
この回答への補足
考えていただきありがとうございます。
玉の総数をMとして、M>>Nの状況を考えています。M→∞と思っていただいてもOKです。
ですから、黒玉の個数はMpで、白玉の個数はM(1-p)です。
説明不足の点は無かったと思いますが、、、。
よろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 数学の問題です。 問1: ある(人数の非常に多い)集団から無作為に6名を選んで身長を測ったところ、そ 2 2022/12/09 12:03
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- 数学 確率 箱の中に赤玉が3個、白玉が3個、青玉が3個入っている。この箱の中から玉を一個ずつ取り出し全ての 4 2023/01/27 18:35
- 大学受験 学習院大 数A 確率の問題について質問です。 2 2023/06/02 15:53
- 数学 数学A 確率 白玉5個、赤玉n個の入っている袋がある。 この袋の中から、2個の玉をとりだすとき、白玉 4 2023/04/22 15:18
- 統計学 期待値を求める問題はとりあえず確率を全部出せばいいってことで答えは出せるのですが、なぜそれぞれの確率 2 2023/07/08 23:11
- 統計学 確率統計です。赤玉が 3 個、白玉が 7 個入った袋がある。 3 2022/05/07 08:36
- 数学 高校数学の確率についてです! 1 2022/05/26 17:29
- 数学 数学の質問です。 以下の2つの例題を①②と分けます。 ①2/4×2/4ではない理由。 ②0は分母が1 2 2023/07/05 15:46
- 数学 【 数A 条件つき確率 】 問題 2つの箱A、Bがあり、箱Aの中には赤玉が4個、白玉が3個入っている 6 2022/10/07 21:15
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
暗闇で2人が出会うには両方動...
-
確率0.02%って10000人に2人です...
-
一般常識を教えてください。1割...
-
30%の確率が5回連続で起きない...
-
「○○通りのパターンがある」の...
-
75%を3回連続で引かない確率
-
AとBが2回ジャンケンをします。...
-
数学の質問です。 一枚の硬貨を...
-
血液型がA型とO型の両親では...
-
五分を6回連続で外すのはなん...
-
「天文学的に低い確率」とは?
-
確率の問題について
-
確率の問題で困ってます。
-
子供が親より先に死ぬ確率は計...
-
中学校数学での確率問題への解...
-
確率
-
6人でジャンケンをした時、1人...
-
高校1年生です。 夏休みの宿題...
-
コインを3回投げた時の確率の...
-
生年月日が同じ人ってどのくら...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
「○○通りのパターンがある」の...
-
確率(高校数学A)
-
30%の確率が5回連続で起きない...
-
一般常識を教えてください。1割...
-
75%を3回連続で引かない確率
-
今の日本は自由社会なハズなの...
-
五分を6回連続で外すのはなん...
-
確率
-
数学の質問です。 一枚の硬貨を...
-
立体四目並べの引き分けになる確率
-
確率0.02%って10000人に2人です...
-
高校数学 確率の問題です。
-
反応速度や濃度は、大きいor小...
-
P(A|B)などの読み方
-
子供が親より先に死ぬ確率は計...
-
「天文学的に低い確率」とは?
-
当たりが4本入った10本のくじが...
-
AとBが2回ジャンケンをします。...
-
高校1年生です。 夏休みの宿題...
-
次の硬貨を同時に投げる時表の...
おすすめ情報