X
トップページ⇒数学
1002コメント451KB
ベイズの統計学を学び始めたんだけど
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
垢版 |
2017/12/03(日) 00:52:27.23ID:v3VGsge3
信用に値するのか疑問です。
人工知能とかではなく日々の動機付けに利用する予定です
0675１３２人目の素数さん
垢版 |
2018/07/04(水) 20:33:40.52ID:VKaeCrJi
時系列データでデータが下降傾向にあるってことを確認する統計手法ってあるの?
二項分布でp＝1/2として検定するとか?

例えば5回連続以上で下降になれば微妙なさげでも下降傾向があると判定できる?
0676１３２人目の素数さん
垢版 |
2018/07/04(水) 20:35:41.15ID:VKaeCrJi
5回連続で下降になると(1/2)^5＜5%になるからって意味。
0678１３２人目の素数さん
垢版 |
2018/07/05(木) 10:06:33.58ID:VNxH3Ygp
同窓会に各人、景品を持ち込む。
全体を集めてクジで持ち帰る景品が決まる。
自分の景品を持ち帰ることになる人数の期待値はいくらか？
0685１３２人目の素数さん
垢版 |
2018/07/05(木) 19:39:42.78ID:q7mJJjLt
>>678
kは0 からｎまでとして

k*nCk*(1/n)^k*(1-1/n)^(n-k)

がｎによらず１になることが示せれば終わり。
0687１３２人目の素数さん
垢版 |
2018/07/05(木) 23:55:11.14ID:0BMJ/kwg
そもそも一人あたり何個の景品を持って来るの？
一人一個持ってきて、分割せず一個を渡すなら、期待値は一人じゃないの？
0688１３２人目の素数さん
垢版 |
2018/07/06(金) 14:37:02.37ID:qCH50YBb
インフルエンザの迅速キットは特異度は高いが感度は検査時期によって左右される。
ある診断キットが開発されたとする。
このキットは特異度は９９％と良好であったが、
感度については確かな情報がない。
事前確率分布として一様分布を仮定する。
５０人をこの診断キットで診断したところ４０人が陽性であった。
この診断キットの感度とその95%CI、及び母集団の有病率とその95%CIは？
0689１３２人目の素数さん
垢版 |
2018/07/06(金) 20:27:38.96ID:qCH50YBb
１回のじゃんけんで決まる勝者の数が最大になるのは何人でじゃんけんをしたときか？

計算していたら４人になったのは意外。
0691１３２人目の素数さん
垢版 |
2018/07/06(金) 22:07:41.15ID:qCH50YBb
>>689
あいこのときは勝者が出るまでやり直すなら
ｎ人のジャンケンなら期待値は n/2人
0692１３２人目の素数さん
垢版 |
2018/07/07(土) 00:29:03.30ID:+asIctIM
>>689
２〜１０人で全員で一度だけじゃんけんをしたときの勝者の数の期待値は

6 / 9
27 / 27
84 / 81
225 / 243
558 / 729
1323 / 2187
3048 / 6561
6885 / 19683
15330 / 59049

[1] 0.6666667 1.0000000 1.0370370 0.9259259 0.7654321 0.6049383 0.4645633 0.3497942
[9] 0.2596149

４人のときが84/81で最大値。
0693１３２人目の素数さん
垢版 |
2018/07/07(土) 00:39:16.72ID:+asIctIM
>>691
あいこはやり直しでの
ジャンケンのシミュレーション、１は3に勝ち、3は2に勝ち、2は1に勝つ
ｘは１，２，３　の並び
1回のジャンケンでの勝者の数を返す
Win <- function(x){# 1 beats 3, 3 beats 2, 2 beats 1
if(length(unique(x))!=2 ) return(0) # no winner
u=sort(unique(x))
if(all(u==c(1,2))) return(sum(x==2))
if(all(u==c(2,3))) return(sum(x==3))
if(all(u==c(1,3))) return(sum(x==1))
}


Jnk.sim <- function(n){
x=sample(1:3,n,replace = TRUE)
while(Win(x)==0){ # repeats when no winner
x=sample(1:3,n,replace = TRUE)
}
Win(x)
}
mean(replicate(1e3,Jnk.sim(10)))

10人でのシミュレーションで期待値は約５
> mean(replicate(1e3,Jnk.sim(10)))
[1] 4.99
0695１３２人目の素数さん
垢版 |
2018/07/29(日) 02:17:03.15ID:YeVWV6wk
　　　　　　　　,,＿＿,,
　　　　　　 ／ 　　　 ｀､
　　 　　　/　　　　　　　ヽ
　　　　　/　●　　　 ●　|
　　　 ／l　 '''''　し　 ''''''　|
　　 /　　l　　　＿＿.　　 |
　　 l　　/ヽ＿ ` --' ＿ノ
　　 ＼　　　　　 ￣　　ヽ∩
　　　　⌒l 　　　　　　　l三 |
　　　　　 |　　　　　　　 ヽ.__|
0696１３２人目の素数さん
垢版 |
2018/07/31(火) 15:04:18.63ID:1kdENkvf
エーザイのアルツハイマー新薬のベイジアン解析の資料読んで
俺に分かりやすく説明してください
0698１３２人目の素数さん
垢版 |
2018/08/04(土) 18:11:25.52ID:1Abx8ltf
高校で一般人も習う条件付き確率になぜ変な名前をつけて喜んでんの?
0699１３２人目の素数さん
垢版 |
2018/08/05(日) 10:32:08.32ID:2T+Vs/Wj
>>698
なんでも確率変数にできるから。
p値の信頼区間(ベイズでは信用区間と呼ぶ)もだせるよ。
0700１３２人目の素数さん
垢版 |
2018/08/05(日) 12:24:22.57ID:txWpcULJ
嘘には三つの種類がある
ただの嘘と
真っ赤な嘘と
統計だ
ってアカポス取れない数学科が統計に流れて薬学だの保険屋になって詐欺をするんだな
0701１３２人目の素数さん
垢版 |
2018/08/05(日) 13:19:41.73ID:hRb7cZuC
統計は嘘をつくための道具だよ
いかに上手く丸め込むか
尤もらしいことを言えた者が勝ち
最尤推定なんて
「もっとももっともらしい」
ってまさにそれだし
0702１３２人目の素数さん
垢版 |
2018/08/05(日) 13:20:30.13ID:hRb7cZuC
それでも推定はまだ可愛げがある
検定お前はダメだ
0706１３２人目の素数さん
垢版 |
2018/08/10(金) 08:07:29.43ID:aM6s4khI
1964年、茨城県南の国道バイパスで突如自家用車が
消失するという事件が起こった
毎日新聞でも取り上げられたこの不可解な失踪事件は
今もなお未解決である
0707１３２人目の素数さん
垢版 |
2018/08/11(土) 02:40:19.11ID:oC/pP3cc
検定はだめだめ言われるけど
なんで未だにどの教科書にも基礎事項としてでてくるんだ?
0710１３２人目の素数さん
垢版 |
2018/08/11(土) 12:56:56.92ID:/7veEAAF
>>707
嘘が上手かったから
0711１３２人目の素数さん
垢版 |
2018/08/12(日) 07:03:58.48ID:lXwEoafU
>>710
メンデルの法則の実験データは法則に合致し過ぎて捏造の疑いを指摘したのがフィシャーだったかな。
0712１３２人目の素数さん
垢版 |
2018/08/12(日) 07:11:07.05ID:lXwEoafU
統計による「検定」
サイコロを２回振ったら順に1,2であった。
その確率は(1/6)*(1/6)=1/36=0.027 < 0.05だから
このサイコロはイビツである。
100人に一人が当たるくじを1本太郎君が引いたら当たった。
この確率は0.01 < 0.05だからこのくじはイカサマである。
0713１３２人目の素数さん
垢版 |
2018/08/12(日) 07:45:35.34ID:ns/pkk0J
>>712
正にこれ
なんで5％しか起こらないことが起こったからといって
仮定を否定する根拠になるわけがない
しかも
改めた結論を元に別のことをまたまた検定してとか
検定に検定を重ねまくることが良くある
1枚でも眉唾な曇りガラスを
何枚も重ねて得られる結論って意味ないだろ
0714１３２人目の素数さん
垢版 |
2018/08/12(日) 11:32:43.28ID:b2zN8ddo
>>712
まさにコレか?

検定に問題あるといえこれは検定を理解してない
ダメな統計の典型でしょコレ
検定力全く考慮してないし
というか科学を理解してないような
0715１３２人目の素数さん
垢版 |
2018/08/12(日) 11:41:59.92ID:b2zN8ddo
＞検定に検定を重ねまくることが良くある
これはやってる本人が問題性に気付かないでやったり悪用するなら問題だけど
問題性に気づいてて探索的にやるのは別に問題ない
検討したい仮説がみつかれば
それを確かめるための実験をくむなり新しいデータをとって追試することが大事なわけで
検定でやる必要はないけど
再現性こそ科学をささえるものなのではないかと思うけど
0716１３２人目の素数さん
垢版 |
2018/08/12(日) 11:44:36.37ID:b2zN8ddo
検定の問題ってのは結局
使う側にも結果を見る側にも誤解をあたえやすいってのが問題なのではないかと
0717１３２人目の素数さん
垢版 |
2018/08/12(日) 12:12:05.87ID:tcAPP7kH
コインを５回投げたら全部、表であった。

0.5^5　<　0.03125なのでこのコインはイカサマ。

コインを５回投げたら表裏表裏表であった。
0.5^5 < 0.03125なので　このコインもイカサマ。
0718１３２人目の素数さん
垢版 |
2018/08/12(日) 12:16:07.59ID:tcAPP7kH
(脱字修正)


コインを５回投げたら全部、表であった。

0.5^5　=　0.03125 < 0.05 なのでこのコインはイカサマ。


コインを５回投げたら表裏表裏表であった。

0.5^5 = 0.03125 < 0.05 なので　このコインもイカサマ。
0719１３２人目の素数さん
垢版 |
2018/08/12(日) 12:36:00.12ID:tcAPP7kH
p値とは、その事象以下の確率でしか起こらない現象の確率の総和である。
　　　　　　　　　　　　↑↑

サイコロをふって１の目が１０回でた。その確率は(1/6)^10

6 2 4 5 2 5 1 3 6 2　の順にでる確率も(1/6)^10である。

1 4 2 4 2 2 4 4 1 1　の順にでる確率も(1/6)^10である。

....

この順列は6^10通り存在するｔ。

これを全部加算すると ｐ値は(1/6)^10 * 6^10 = 1

ゆえに、いかなるサイコロも統計的にイカサマとはいえない。
0720１３２人目の素数さん
垢版 |
2018/08/12(日) 12:40:57.07ID:b2zN8ddo
>>718
仮説が立ってない
データで検証しうる仮説は何をたてたの?
検定はまず仮説からたてないとダメ
君がやってるのは統計的仮説検定ではないよ
0721１３２人目の素数さん
垢版 |
2018/08/12(日) 12:44:40.42ID:b2zN8ddo
検定に置いて大事なのは
サンプルサイズ、検定力、有意水準、効果量
仮説もきちんと事前にたてるのプロセスとして大事

このあたりが欠けてる検定は何も意味ないよ
0722１３２人目の素数さん
垢版 |
2018/08/12(日) 12:46:12.27ID:b2zN8ddo
あれ?このやりとり前にもやったような。。。
同じやつが書いてるのならループするだけだな
落ちるわ
0723１３２人目の素数さん
垢版 |
2018/08/12(日) 12:47:29.69ID:62kxiEQs
帰無仮説は表がでる確率=1/2
サイコロの各々の目がでる確率=1/6
じゃね？
0724１３２人目の素数さん
垢版 |
2018/08/12(日) 14:11:00.62ID:ns/pkk0J
そこが問題じゃない
1/6であることが起こったから
帰無仮説を棄却するというのが根拠ないってこと
それから
前に上げたこと以外に
両側検定なら対立仮説が棄却されるのに
片側検定なら採択されるとか
恣意的に両側片側を選んで上手く欺すことが可能なことも多い
だいたい
片側検定では正の範囲あるいは負の範囲だけ考えるんだから
つまりは条件付き確率を考えるようなもの
であるからして有意水準は半分にしないとおかしい
両側検定の有意水準÷２＝片側検定の有意水準
であるべきなのに同じ有意水準でいいとしてるのが
この欺瞞の元凶
0725１３２人目の素数さん
垢版 |
2018/08/12(日) 14:14:37.21ID:tcAPP7kH
>>719

サイコロの各々の目がでる確率=1/6を事前分布として

サイコロをふって１の目が１０回でたら、１の目のでる確率のモード値は0.7

6 2 4 5 2 5 1 3 6 2　の順にでたら、１の目のでる確率のモード値は0.07

98％CIとともに図示すると、http://i.imgur.com/CTvq0lL.jpg

ｐ値なんぞ判断には不要。
0726１３２人目の素数さん
垢版 |
2018/08/12(日) 15:14:16.30ID:YrquRjB0
ベイス？？
0727１３２人目の素数さん
垢版 |
2018/08/12(日) 16:01:37.80ID:tcAPP7kH
サイコロのある目のでる確率の９５％CIが０．１０〜０．２０に含まれたときにその目に関しては歪ではない、と定義する。

サイコロを１０００回投げて１〜６の目の回数が
　
　309　 251 　196　 151 　 49 　 44

であったときの各々の目のでる確率分布は以下のようになる。

95%HDI　と 0.10-0.20を比べて１、２、５、６は歪、３はどちらとも言えない、４は歪でないと判断できる。

https://i.imgur.com/w4YD0t9.jpg
0728１３２人目の素数さん
垢版 |
2018/08/12(日) 17:54:58.31ID:b2zN8ddo
>>723
1のでる確率＝1/6じゃなくて
それぞれの確率＝1/6?
全部の確率について仮説たててるの?
P1=P2=P3=P4=P5=P6＝1/6が仮説?
それどう考えても自由度足りないよね?

その仮説考えたら次にすべきは
どういうデータを取れば論証できるだろうかになるわけだが
2回ふってでた目をみてそれの確率求めることで一体なにが論証できると思う?
それって何が出ても1/36になるけどそれが検証可能なデータだと思う?
そこに何の意味があるの?

そもそもなぜそこで、n＝2の母比率検定じゃなくて、
N=1の一回目の確率x二回目の確率という事象の確率を検討するの?

検定いぜんに実験の企画段階がめちゃくちゃ可笑しいとおもうんだけどな
0729１３２人目の素数さん
垢版 |
2018/08/12(日) 18:04:24.47ID:b2zN8ddo
>>724
そうかね自分は科学的なしゅほうだとおもうけどな

科学的な手法ってのは
問をたて仮説をたて
それを検証できるデータをとり
そこから仮説がどうなのか検証し
どうやらこの仮説ははいろいろ説明できるようだぞというなら
とりあえず残しておく
いやいろいろ説明できない点があるぞというときは
とりあえず研究対象から外す

ほんでほかのやつが残った仮説にもどうやら問題あるようだ
こっちのほうがうまく説明できるようだと別の仮説を立てる

しかしかし、以前、は研究対象から外したものがいろいろ調べてみると正しいようだというデータをあつめるやつがいる


これって科学的仮説がたどっていく過程そのもの何だがな
仮説検定はその一旦を担ってるいる似すぎないだけであって
それで全てというものではない。

科学なんてあくまで仮説の積み上げで
仮説が多くの合意を得ることもあるが、その仮説が正しいとするのは結構曖昧にすぎず
不安定なものに過ぎない

根拠がないと言うが統計によってえられた結論ってそもそもどんなものも大した根拠なくて主観がどこかにはいってるよね。
0730１３２人目の素数さん
垢版 |
2018/08/12(日) 19:39:29.57ID:b2zN8ddo
片側両側については一理ないこともないとおもう

しかし、片側にするってことは片側は起きないだろうという論理的説明を一応出来る状態にあるということ
となると対立仮説も片側に発生する可能性がたかくなるわけでしょ。
だったらかりに推定する母数が平均値なら、きむ仮説でたてた平均値にちかい真の平均値も増えるだろう
だったら採択いきは狭くして棄却域増やさざるを得なくなる
という仮説

こういう宗教的論争を多数生み出すてんで結局仮説検定は問題あるんだけどね。
0731１３２人目の素数さん
垢版 |
2018/08/12(日) 21:19:55.98ID:tcAPP7kH
帰無仮説をサイコロの各々の目がでる確率=1/6とするとき

１の目が続けて何回でたら、帰無仮説を棄却する？

棄却されたときに１の目がでる確率はどれくらいのなのか？

1/6.0000001 なのか 1/2なのか？　p値は何かを語るだろうか？

むしろ、１の目のでる確率の95%CIが1/12から1/3なら許容範囲、1/12以下や1/3以上なら１の目に関して歪なサイコロと判断する方が実用的だと思う。


１の目が２〜９回連続したときの１の目のでる確率をグラフにすると次のようになる。

http://i.imgur.com/nUd5UMO.jpg

ROPE：Range Of　Practically Equalの略　（1/6の半分および2倍の1/12〜1/3とした）

この結果から７回以上連続すれば、歪なサイコロと呼べると思う。
0733１３２人目の素数さん
垢版 |
2018/08/12(日) 22:14:39.76ID:b2zN8ddo
>>732
まじで?
よく知らんけど
そんだけのおもい仮説検証で
N=1で自由度足りるの?
0735１３２人目の素数さん
垢版 |
2018/08/12(日) 22:17:27.59ID:b2zN8ddo
N=1で6個の仮説を検証できるて凄くないか?
天才っているんだなw
どう考えても無理を可能にする天才
すげぇわ
0736１３２人目の素数さん
垢版 |
2018/08/12(日) 22:18:50.58ID:b2zN8ddo
自由度どう考えても足りないんだけどどうやって補完してるの?
0737１３２人目の素数さん
垢版 |
2018/08/12(日) 22:21:19.95ID:tcAPP7kH
５回1の目が出たときの事後分布はこんな感じでディリクレ分布がMCMCされるよ。

もちろんpi.1+pi2+...+pi.6=1の縛りはある。

pi.1 pi.2 pi.3 pi.4 pi.5 pi.6
[1,] 0.3283620 0.051967675 0.03618622 0.16401712 0.367065379 0.052401570
[2,] 0.6824722 0.009042027 0.02910249 0.03210041 0.018967898 0.228314983
[3,] 0.5938329 0.015613992 0.04579314 0.22754563 0.027786023 0.089428325
[4,] 0.5801030 0.033073118 0.01888069 0.19814067 0.001983984 0.167818557
[5,] 0.5501332 0.056247586 0.04896485 0.02131748 0.090452926 0.232883935
[6,] 0.3149563 0.051999120 0.04049730 0.32030618 0.039353349 0.232887755
[7,] 0.5496275 0.154893875 0.23899680 0.02146069 0.028756420 0.006264765
....
0738１３２人目の素数さん
垢版 |
2018/08/12(日) 22:21:43.52ID:b2zN8ddo
N=1で統計的に正しいですって結果出されても説得力まったくないなぁ
統計は詐欺だと立証されてるようなもんだけどw
0739１３２人目の素数さん
垢版 |
2018/08/12(日) 22:23:22.50ID:b2zN8ddo
ごめん詳しいこと知らないので分かるようにせつめして
N=1で六個の仮説についてはんだんできるってすごくね?
どういう理屈?
0740１３２人目の素数さん
垢版 |
2018/08/12(日) 22:24:15.64ID:tcAPP7kH
>>735

それがベイズ統計の醍醐味でもあり、胡散臭さでもある。

こういうのね

事前分布を決めてしまえば
まだ一発も撃ったことのない0発0中のゴルゴ16の命中期待値
のような、データ数が少ないどころか0個の場合でも算出・結論できる
0741１３２人目の素数さん
垢版 |
2018/08/12(日) 22:27:04.16ID:b2zN8ddo
>>740
よくわからんけどそうなると統計的仮説検定とあんまかわらんような気がするけど
0742１３２人目の素数さん
垢版 |
2018/08/12(日) 22:29:07.04ID:b2zN8ddo
N=1ってさ人文科学の事例研究なんだよね
そこから一般論を導き出すようなかんじ
人文科学ってこれだからばかにされるんだけど
コレをうわまわる手法が統計にはあるのかすげぇなw
0743１３２人目の素数さん
垢版 |
2018/08/12(日) 22:37:10.89ID:tcAPP7kH
サイコロの目の出る確率が目の数に比例する、


2は1の目の出る確率の2倍

3は１の目の出る確率の3倍
....

とかいう事前分布のとき

1の目が続けて5回でたときの1の目のでる確率の事後分布のグラフも書けるよ。
0744１３２人目の素数さん
垢版 |
2018/08/12(日) 22:40:11.41ID:b2zN8ddo
>>743
グラフ書けるかどうかじゃなくて
バカにも分かるように説明して下さい
0745１３２人目の素数さん
垢版 |
2018/08/12(日) 22:41:21.70ID:b2zN8ddo
N=1で6個の仮説の真実説位を検証できる理屈についておしえて
0746１３２人目の素数さん
垢版 |
2018/08/12(日) 22:45:25.16ID:tcAPP7kH
何を持って歪なサイコロとするか、各々の目のでる確率の分布をどう仮定するか、というだけの話。
0747１３２人目の素数さん
垢版 |
2018/08/12(日) 22:48:05.81ID:tcAPP7kH
これを多項分布に応用しただけの話。

ゴルゴ13は100発100中
ゴルゴ14は10発10中
ゴルゴ15は1発1中
とする。

各々10000発撃ったとき各ゴルゴの命中数の期待値はいくらか？
0748１３２人目の素数さん
垢版 |
2018/08/12(日) 22:48:28.76ID:b2zN8ddo
>>746
よくわからんが
有意水準の設定と変わらんね

しょうじきN=1で六個の仮説を検定できるなんて
統計的検定より胡散臭いけど結論がよく分からん
0749１３２人目の素数さん
垢版 |
2018/08/12(日) 22:51:42.49ID:b2zN8ddo
>>747
それとN=1で六個の仮説を論証出来るのと同関係あるの?

N=1と多項分布がどうかんけいあるのよくわからんえけど
0750１３２人目の素数さん
垢版 |
2018/08/12(日) 22:55:25.81ID:tcAPP7kH
>>747
これは命中確率の事前分布を一様分布として

ゴルゴ13は100発100中だったときに

命中確率の分布がどう変わるか、を計算することになる。

ベイズ統計って　事前と事後で　確率分布がどう変わるか（relocation of credibility）を探る手法。
0751１３２人目の素数さん
垢版 |
2018/08/12(日) 23:01:14.39ID:tcAPP7kH
>>749
すべての目の出る確率が等しいと事前確率分布を設定して
1回サイコロをふったら１の目がでた。
各々の目の出る事後確率分布はどうなるか？

これだけの話。

>727で数値を

1 0 0 0 0 0

に置き換えるだけ。
0752１３２人目の素数さん
垢版 |
2018/08/12(日) 23:01:25.62ID:b2zN8ddo
>>750
自分に返答するのは意味不明だけど
けっきょくN=1で六個の仮説を論証できる原理って何なの?
バカにも
IME okashiku natta
sayonara
0753１３２人目の素数さん
垢版 |
2018/08/12(日) 23:38:31.35ID:b2zN8ddo
直ったわ
とりあえずN=1で6個の仮説証明できる原理を教えて
0755１３２人目の素数さん
垢版 |
2018/08/13(月) 08:20:30.92ID:vxHeSM3M
１から６の目の出る確率がパラメータα1=α2=...=α6=1のディリクレ分布に従うを事前確率分布とする（これは信仰）。

1回サイコロをふったら１の目がでた。
各々の目の出る事後確率分布はどうなるか？

をMCMCして出すだけ。

１の目の出る確率は平均0.287[95%CI 0.0161-0.593]
２〜５の目の出る確率は各々平均0.143[95%CI 0-0.392]

と計算されたが、この分布で歪がどうかは、

何を歪と判断基準にするのか、というだけのお話。
0756sage
垢版 |
2018/08/17(金) 03:36:43.23ID:MEYBPptz
確率はすべて条件つき確率であり、この点、事後確率をかんがえるベイズの理論は正しい。
0757１３２人目の素数さん
垢版 |
2018/08/17(金) 09:12:30.98ID:5cP2NSmE
>>756
おまえなあ、sageを名前のほうに書いてるのは故意にやってる？
笑えないんだけど？

＞　ベイズの理論は正しい。

その場合の「正しい」の定義が不明なので意味はない。
意味はないことを書くのが趣味なのかね？
数学で示すべき。
数学板なんだから。
0758１３２人目の素数さん
垢版 |
2018/08/17(金) 09:51:51.87ID:Xs+I9BdE
誤差の分布が正規分布という事前確率で頻度主義統計も議論していると思う。
0759１３２人目の素数さん
垢版 |
2018/08/18(土) 10:01:05.05ID:w7FxWks2
どの面も出るのが同様に確からしい
6面ダイスを独立に2回振った時に
少なくとも一回は1の目が出る確率は
いくらですか？
0760１３２人目の素数さん
垢版 |
2018/08/18(土) 12:43:32.76ID:mXMMyLFy
>>759
＞同様に確からしい

どの程度、同様に確からしいのを事前確率分布して計算するのがベイズ統計。

ディリクレ分布でパラメータを(1,1,1,1,1,1)とするのか、(10,10,10,10,10,10)とするのか、(100,100,100,100,100,100)とするのかで

少なくとも一回は1の目が出る確率分布は変わる。

図示すると、以下の通り、http://i.imgur.com/J1XUpAw.jpg
0761１３２人目の素数さん
垢版 |
2018/08/18(土) 15:36:03.95ID:zXb74bfP
>>754
つまり宗教なの?
0762１３２人目の素数さん
垢版 |
2018/08/18(土) 16:32:44.71ID:+ZAzv04a
>>759
＞同様に確からしい

どの程度、同様に確からしいのかを事前確率分布して計算するのがベイズ統計。

ディリクレ分布でパラメータを(1,1,1,1,1,1)とするのか、(10,10,10,10,10,10)とするのか、(100,100,100,100,100,100)とするのかで

少なくとも一回は1の目が出る確率分布は変わる。

図示すると、以下の通り、http://i.imgur.com/J1XUpAw.jpg
0763１３２人目の素数さん
垢版 |
2018/08/18(土) 16:34:57.61ID:+ZAzv04a
>>761
日本人女性の身長の平均値は1〜２ｍの間にある、
というのもまあ、信仰と言えなくもない。
0766１３２人目の素数さん
垢版 |
2018/08/20(月) 15:01:08.82ID:2dQLFNqz
ジョーカーを除いたトランプ５２枚の内、ダイヤが３９枚
ハートが１３枚あるとする
この中から１枚のカードを抜き出し、表を見ないで箱の中にしまった
そして、残りのカードをよく切ってから３８枚抜き出したところ、
３８枚すべてがダイヤであった
このとき、箱の中のカードがダイヤである確率はいくらか
0767１３２人目の素数さん
垢版 |
2018/08/20(月) 19:26:00.92ID:jCZwkgjR
>>766
> 39/52 * choose(38,38)/choose(51,38) / ( 39/52 * choose(38,38)/choose(51,38) + 13/52 * choose(39,38)/choose(51,38) )
[1] 0.07142857
0769１３２人目の素数さん
垢版 |
2018/08/20(月) 19:54:20.26ID:jCZwkgjR
D=39
H=13
T=D+H
# choose(n,r) : nCr

39/52 * choose(38,38)/choose(51,38) / ( 39/52 * choose(38,38)/choose(51,38) + 13/52 * choose(39,38)/choose(51,38) )

(D/T * 1/choose(T-1,D-1)) /((D/T * 1/choose(T-1,D-1)) + H/T * choose(D,D-1)/choose(T-1,D-1))

(D/T)/(D/T+H/T*D)

D/(D+H*D)

1/(1+H)
0770１３２人目の素数さん
垢版 |
2018/08/20(月) 19:57:51.64ID:jCZwkgjR
>>766
これ面白いな。

ダイヤがD枚
ハートがH枚あるとする
この中から１枚のカードを抜き出し、表を見ないで箱の中にしまった
そして、残りのカードをよく切ってからD-1枚抜き出したところ、
すべてがダイヤであった
このとき、箱の中のカードがダイヤである確率は　1/(H+1)
0771１３２人目の素数さん
垢版 |
2018/08/20(月) 20:51:49.87ID:2dQLFNqz
ハート１３枚、ダイヤ３９枚の合計５２枚のトランプカードから
一枚のカードを表を見ないで箱に入れる

Ωの部分集合を事象と言う
Ω自身は全事象と言う

Ω＝｛ハート，ダイヤ｝となる

各 i （１≦i≦４） が根元事象である

ハートが出るという事象A＝｛ハート｝で確率P(A)は

P(A)＝１／４　となる

最初に箱に入れた時を i
山札をシャッフルしてダイヤがｎ枚出た後を ｊ として

箱の中のカードがハートであるという事象Aを考える.

A＝｛（i，j）|　i または j がハート｝

Ω＝｛（i，j）|１≦i≦４，１≦j≦５２−ｎ｝となり

この２０８−４ｎ通りの各要素が根元事象

ダイヤが出る枚数はｎ＝３８

Ω＝｛（i，j）|１≦i≦４，１≦ｊ≦５２−ｎ｝から

#A＝４ｘ（５２−ｎ）−３ｘ（５１−ｎ）

　　＝２０８−４ｎ−１５３＋３ｎ

　　＝５５−ｎ

#Aは事象Aに含まれる要素の個数

ハートである確率は

P(A)＝（５５−ｎ）／（２０８−４ｎ）＝１７／５６

ダイヤ以外である確率P(Ｘ)＝P(A)

ダイヤである確率は

∵ｑ＝１−｛（５５−ｎ）／（２０８−４ｎ）｝＝３９／５６

最初に箱の中にしまったカードが
ダイヤである確率は

P（B）＝３／４＝４２／５６
■ このスレッドは過去ログ倉庫に格納されています
ニューススポーツなんでも実況