黒木玄 Gen Kuroki(@genkuroki)/2020年05月/Page 11

#統計頻度論とベイズが別だと思うのはおかしいので、そういうことをするなら、最尤法(やMAP法)側の用語も変えちゃった方がスッキリするかも。

ただし、Stanのmodelで記述されるモデルだけで、対応する最尤法・MAP法側が唯一に決まらないことに注意。予測分布が唯一に決まらないのと同じ。 twitter.com/katzkagaya/sta...

タグ：統計

posted at 00:44:52

#統計 MCMCに必要な情報を記述するためのモデル情報は、何から何を予測したいかに関する情報は持っていない。

だから、MCMCに必要な情報を記述するためのモデル情報だけに注目すると、予測分布の概念も自然に出て来なくなる。(同じ理由で対応する最尤法・MAP法も決まらなくなる。)

タグ：統計

posted at 00:44:52

#統計再度例で説明。Y_kは観測されたデータのとき、モデル

z ～ p₁(z)
w_k ～ p₂(w_k|z)
Y_k ～ p₃(y_k|w)
(k=1,…,n)

に対応する予測分布として、

① p(y|z)=∫p₃(y|w) p₂(w|z)dwの事後分布平均

② Π_{k=1}^n[p₃(y_k|w_k) p₂(w_k|z)]の事後分布平均

の少なくとも二通りが考えられる。

タグ：統計

posted at 00:50:37

#統計

① に対応するのは、Π_{k=1}^n p(Y_k|z) に関する最尤法・MAP法

② に対応するのは Π_{k=1}^n[p₃(Y_k|w_k) p₂(w_k|z)] に関する最尤法・MAP法

w_kで積分するかしないかの違い。

予測分布の定義を与えることと、対応する最尤法・MAP法の定義を与えることは同じ。

タグ：統計

posted at 00:56:33

#統計普通の状況で興味があるのは、w_kで積分した①の方です。

しかし、WAICやLOOCVを計算してくれる既存のパッケージを余り考えずに使うと、②を前提に計算してしまうことになる。自分で数値積分すれば同じパッケージで①も可能なのですが。

タグ：統計

posted at 01:01:32

Otepipi @Otepipipi

はてなブログに投稿しました #はてなブログ, #Julia言語, #微分方程式
JuliaでCobweb plotを作成する - システムとモデリング
otepipi.hatenablog.com/entry/2020/05/...

タグ： Julia言語はてなブログ微分方程式

posted at 07:15:14

io302 @io302

『米田きよし氏、自称「カナダの小児救命救急医」（2011年）』
『「医師法違反（医師名称使用）の疑い」で逮捕』

『森口尚史氏、自称「ハーバード大客員講師」（2012年）』
『騒動で懲戒解雇』

なつかしの医療デマの人達ｗ

president.jp/articles/-/352...

タグ：

posted at 07:53:38

シュー @shu_n148

普通でも危険を伴うランニング中にマスク着用をエチケット化するのは、体温の上昇などによる事故の増加の原因になりかねません。屋外で一人でできるランニングは距離を保てば安全です。マスクをしないように呼びかけた方がいいと思います。 twitter.com/kenjiimoriwaki...

タグ：

posted at 07:54:27

io302 @io302

『齋藤ウィリアム浩幸氏、自称「カリフォルニア大学ロサンゼルス校（UCLA）卒医師、サイバーセキュリティ専門家」（2017年）』
『経歴詐称疑惑が指摘』

『内藤理恵氏、自称「シンガポール医大卒の外科医、読者モデル」（2019年）』
『携行する医師免許が偽造』『婚約も解消』

皆いろいろ失うんだね

タグ：

posted at 07:54:47

データのロード以外はJuliaのLightGraphs.jlが最速？
www.timlrx.com/2020/05/10/ben... pic.twitter.com/pLcnbBsX2y

タグ：

posted at 08:30:20

みやびの@miyabino.py @miyabikno

プログラミング言語「Julia」の導入方法です。

Juliaの導入と実行方法〜計算と速度に優れたプログラミング言語を使おう〜 se.miyabikno-jobs.com/julia-introduc...

タグ：

posted at 08:34:44

LightGraphs.jlはマルチスレッドでの性能も良さそう（他のツールが載ってないのはサポートしてないから）。 pic.twitter.com/Nc95i91C57

タグ：

posted at 08:41:24

非公開

タグ：

posted at xx:xx:xx

OokuboTact　大久保中二病中年 @OokuboTact

タグ：統計

posted at 09:24:25

＃超算数

掛け算の順序が生まれた理由は、割算九九なのかもしれない

タグ：超算数

posted at 09:24:31

#統計さらに、データの項目で、観測された変数の値が上のモデル f に与えられます。例えばn=3でy_1=Y_1, y_2=Y_2が観測され、y_3はmissingならば、観測されたデータをfに代入して得られるgは

g(y_3,w_1,w_2,w_3,z)=f(Y_1,Y_2, y_3,w_1,w_2,w_3,z)

となる。続く

タグ：統計

posted at 09:29:02

#統計続き。このgをMCMCサンプラーに渡せば、gの定数倍で正規化して得られる確率密度函数を持つ分布のサンプル(i.i.d.の実現値)が(近似的に)得られます。

これは観測データY_1,Y_2から、y_3,w_1,w_2,w_3,zの分布を推定したことになります。ただし、解釈にはいろいろ注意が必要かも。

タグ：統計

posted at 09:35:04

Kazuho Oku @kazuho

（もう長らくいってない）オフィスが27インチだけど、34インチウルトラワイドのほうが断然いい。A4縦を３枚開けるから。真ん中にメインのドキュメント置いて、左右に補助広げられるのは、体にやさしい。
ウルトラワイドじゃない27インチだと、２枚になって、左右どちらかに首がむいちゃう twitter.com/kazuho/status/...

タグ：

posted at 10:08:00

#統計そうやって作った(y_3,w_1,w_2,w_3,z)のサンプルのy_3への射影は、観測データY_1,Y_2からy_3を予測するためのベイズ版予測分布のサンプルになっています。

もしかしたら、この例のように考える方が易しいと感じる人がいるかも。MCMCで直接予測分布のサンプルが得られる。

タグ：統計

posted at 10:08:33

これ。ウルトラワイドはずっと正面を向いてられる。

タグ：

posted at 10:11:36

#統計 #Julia言語のTuring.jlにおけるモデル作成の内部で何をやっているかについては

turing.ml/dev/docs/using...

を参照。lp = … の lp は上の方の例では

(a) log p₃(z)
(b_k) log p₂(w_k|z)
(c_k) log p₁(y_k|w_k)

に対応しています。

a ~ dist(b,c) ↦ log pdf(dist(b,c), a)

タグ： Julia言語統計

posted at 10:48:32

#統計 2つ上の続き。

(w_1, w_2, z)への射影は

尤度函数×事前分布
= p₃(Y_1|w_1)p₃(Y_2|w_2)p₂(w_1|z)p₂(w_2|z) p₃(z)

の場合の事後分布のサンプルになる。

twitter.com/genkuroki/stat...

タグ：統計

posted at 10:53:47

タグ：統計

posted at 10:55:13

@389jan @katzkagaya #統計ど素人過ぎて考えが足りなかった。

より普通で楽な考え方↓

twitter.com/genkuroki/stat...

タグ：統計

posted at 10:56:58

#統計 3つ上の続き

(a) z ～ p₁(z)
(b_k) w_k ～ p₂(w_k|z)
(c_k) y_k ～ p₃(y_k|w)
(k=1,…,n)

のようにモデルを与えることは

(a) log p₃(z)
(b_k) log p₂(w_k|z)
(c_k) log p₁(y_k|w_k)

を与えて足し上げること(その指数函数がすぐ上で訂正した f になる)と同じだという理解は基本的。

タグ：統計

posted at 11:01:08

#統計そういうスタイルでシンプルな良い考え方は、予測したい対象を最初からすべて上の意味でのモデルに含めておくこと。

仮にMCMCなどの方法で予測対象の分布のサンプルが得られたなら、そのサンプルが予測対象のモデルによるシミュレーションになっていると解釈すればよい。

タグ：統計

posted at 11:04:03

#統計で、残された問題は予測誤差の(差の)推定でモデル選択する方法。

上の流れでw_kに関する数値積分抜きにWAICやLOOCVを近似計算するうまい方法があればよいと思う。

単に予測しっぱなしで、予測の相対的妥当性の評価ができないと困る。

安易に交差検証したつもりになると危ないと思う。

タグ：統計

posted at 11:08:21

nagashima m. @_nagashimam

今期、諸事情により小学校免許用の授業を代理で一部持ってるんだけど、補足資料でつらつら書いてたらつい「かけ算の順序を強制するのは子どもの生きる力を削ぐ教育」とか「『はじき』なんて子どもに点を取らせて教師が自己満足に浸るだけの最悪の教育」とか書いちゃった(後悔はしていない

タグ：

posted at 11:48:27

こいなぎ @naginyan135

【級位者向け実戦型詰将棋part196】
この５手詰めなかなか実戦では見えなくない？ pic.twitter.com/EcXt86PDPO

タグ：

posted at 12:13:17

あおじるPPPP @kale_aojiru

ごまふあざらし(GomahuAzaras @MathSorcerer

それは割合や速度をそういう風に教えるからだよね twitter.com/kaburaya75/sta...

タグ：

posted at 12:35:35

非公開

タグ：

posted at xx:xx:xx

OokuboTact　大久保中二病中年 @OokuboTact

#Julia言語

Interact.jl が Jupyter Lab で使えたら誰か教えてください．（PythonのKernelではipywidgetsは動いている模様）

タグ： Julia言語

posted at 13:07:20

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

＃超算数

大正時代に小学校で総九九が採用された時に、割算九九の影響が・・・

＃掛け算の順序の固定の誕生の秘密を考えてみる pic.twitter.com/w8Df7FUG2H

タグ：掛け算の順序の固定の誕生の秘密を考えてみる超算数

posted at 15:33:54

稲葉振一郎 @shinichiroinaba

何かこう永井均の柄谷行人『探究』批判を思い出した。 / 他109件のコメント b.hatena.ne.jp/entry?url=http... “平田オリザ氏炎上にみる、コミュニケーション教育は「他者と相互理解する能力」を害しているという現実｜ヒラヤマ　タカシ｜note” htn.to/2C1Vvojcam

タグ：

posted at 16:08:01

Matt Henderson @matthen2

ごまふあざらし(GomahuAzaras @MathSorcerer

plotting Julia sets on the Riemann sphere pic.twitter.com/wGTmt0OWc5

タグ：

posted at 16:47:25

github.com/JuliaLang/juli...

タグ：

posted at 17:14:18

ReRTs ｴﾇﾇﾇﾇﾇﾇﾇﾝさんは仲間。 #統計

検査を受けた人の感染率と感染者を検査したときに陽性になる割合と非感染者を検査したときに陰性になる割合が分かっているときに、検査陽性者中の感染者の割合と検査陰性者中の非感染者の割合を計算する問題は

小学校レベルの割合の問題

です。続く

タグ：統計

posted at 17:33:41

#統計続き。実際に関連の大学入試問題を小学生に出したら解けました。小学校レベルの割合の概念を実用レベルで習得できているかどうかの確認に陽性的中率を求めさせる問題はちょうど良いレベルだと思う。

もちろんベイズの定理でも解けるけど、小学生レベルの問題であるという自覚は欲しい。

タグ：統計

posted at 17:33:44

#統計上記の{陽性,陰性}的中率を求める問題は、ベイズの定理を使って__も__でも正解できるという理由で、その問題をベイズ統計の話題だとみなすのだろうか？小学生レベルの割合の問題に過ぎないのだが。

小学生も当事者であり、小学生でも理解できる話に過ぎない。

タグ：統計

posted at 17:33:45

#統計ベイズ統計に限らず統計学は「未知の何か」を推測するために使われています。例えば、私がちょっと前にツイートで取り上げた論文

www.ncbi.nlm.nih.gov/pubmed/32046819
www.ncbi.nlm.nih.gov/pmc/articles/P...

ではベイズ統計の方法で(実際には他の方法も使える)、新型コロナの潜伏期間期間の分布を推測しています。続く

タグ：統計

posted at 17:40:32

#統計さすがにそこまで行けば小学校レベルの割合の問題とは言えないし、未知の何かの推測なので相当に非自明なことをやっています。

ベイズ統計の方法が未知の何かを推測するための適切な方法の1つになることの理解は、ベイズの定理のような自明でくだらない話題と違って、ずっと複雑で面倒です。

タグ：統計

posted at 17:43:49

#統計「どうしてベイズ統計の方法が現実の予測に役に立つ道具の1つ(他にも沢山ある！)になり得るのか？」を知らないと、みんなが「ベイズのすごいぞ」と言っているからベイズ統計を理解せずに使う迷惑な人になってしまいます。

ベイズの定理を知っていても必要な理解は全く得られません。

タグ：統計

posted at 17:46:28

#統計ぶっちゃけた話、現時点では、そういう方向の理解には、渡辺澄夫著『ベイズ統計の理論と方法』を読む以外に選択肢がないと思う。いきなりその本を読まずにググって渡辺澄夫さんのウェブサイトを見つけて、講義のスライド(結構かわいい😊)を見るとよいと思う。

タグ：統計

posted at 17:49:04

#統計一般向けのベイズ統計の解説のほとんどが、20世紀の時代遅れの考え方に基いているので要注意です。

タグ：統計

posted at 17:50:29

はっぴぃ @shiawaseahiru

大人猫との接触はNGとシッターさんにうまく伝わってなくて仲良しさせてしまった動画が送られてきました。次からは接触させないようにお願いしました。うぃるちょん人気者。早くみんなと仲良しさせてあげたいのだけど pic.twitter.com/ztBVvOe5ZK

タグ：

posted at 17:50:31

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

#統計 (y₃,w₁,w₂,w₃,z) のサンプルが得られ、その y₃ への射影が確率モデル p(y|z)=∫p₁(y|w)p₂(w|z)dwと事前分布p₃(z)と観測データY₁,Y₂からy₃を予測するためのベイズ予測分布になっていることの手書きのノート。

twitter.com/genkuroki/stat... pic.twitter.com/iTlgfX0Eyg

タグ：統計

posted at 18:24:56

こや @entropiajp

自宅で学校のチャイムを鳴り響かせたらメリハリが出すぎた話-Google Home活用-｜染谷優作 @usakufei #note note.com/usakufei/n/nfc...

家チャイムいいなぁ

タグ： note

posted at 18:25:20

#統計確率が絡むゲームをやりなれていれば、モンティ・ホール問題は自明な問題。ぶっちゃけ小学生レベル。

最初にランダムに3つのうちの1つを選ぶと外れる確率は2/3になる。実際に外れていて、もう1つの外れの位置が分かっていれば、残りの1つを選べば勝てる。この方法で2/3の確率で勝てる。

タグ：統計

posted at 19:03:25

砂糖水 @puzzlexical

万博さんより原理主義の非常に丁寧なレビューをいただきました。
ありがとうございます！ twitter.com/bampaku/status...

タグ：

posted at 19:07:21

#統計陽性的中率の計算もモンティ・ホール問題で3分の2の確率で勝つ方法も小学生レベルの問題です。ベイズの定理の公式を経由して解くしかないと思っている人の理解不足の度合いは結構要注意だと思います。

タグ：統計

posted at 19:09:13

WindowsのCIでだけ謎エラー…

タグ：

posted at 19:10:01

万博@盲学校マジック @bampaku

原理主義はいいぞ。

タグ：

posted at 19:10:55

クロイヨシツネ @YoshitsuneK

@happyduckshige 早く一緒に遊べるようになるといいですね、あっさり受け入れられていて驚きましたw 先輩方が大らかなのですね。

タグ：

posted at 19:15:13

ジョナサン @itai_daboku

結局、ちゃんと確率空間さえ与えられていればベイズの公式は確率論における単なる条件付期待値。

通常ベイズと言われているp(θ|x)∝p(x|θ)p(θ)ってのが論争の的となるのは、そもそも解析の対象となる確率空間が定まっていないってだけの話。

タグ：

posted at 19:46:08

ultraviolet @raurublock

「営業停止で生計を断たれ、もはや一家で野垂れ死ぬしかありません」とブチャラティに訴える店主の絵が頭に浮かんだ www.bbc.com/japanese/featu...

タグ：

posted at 19:51:09

RANGER @narcissina

【佐渡の野鳥紹介 25】オオコノハズク
樹洞や穀物倉庫で暮らすフクロウ。民家のそばでうずくまっている雛がよく観察されます。フォークリフトの運転席に雛が座っていたことさえありました。

昼間は寝ぼけているので、誤って保護されることが多い鳥です。雛を見つけたら、やさしく見守りましょう。 twitter.com/narcissina/sta... pic.twitter.com/0myIcNBe5R

タグ：

posted at 20:12:17

TaKu @takusansu

@sekibunnteisuu @sunchanuiguru #超算数算数教科書の著作関係者の数ですが、2015年度から2020年度では大幅に増えています。

東京書籍
41人→83人

啓林館
59人→92人

学校図書
49人→64人

教育出版
30人→34人

日本文教出版
24人→32人

大日本図書
23人→35人

タグ：超算数

posted at 20:13:31

Ninja DAO | CryptoNi @CryptoNlnjaNFT

各種の時事問題に関する話題と比較すると、統計学っぽい話題が相対的に圧倒的に「平和」なことが際立っていると思った。😊

タグ：

posted at 20:22:23

大体、自分はナショナリズムそのものがあんまり好きじゃないので、日韓両国間の議論はどちらにもあんまり好感を持っていない。ナショナリズムが過剰という意味では、どちらもあるベクトルの同じ端の方にいる。

タグ：

posted at 20:30:17

飯田泰之 @iida_yasuyuki

Ninja DAO | CryptoNi @CryptoNlnjaNFT

報道出たということは人はもう決まってるはずだが誰だろう？　緊縮派でありませんように

→コロナ諮問委員会に経済学者を追加 | 2020/5/11 - 共同通信 this.kiji.is/63252953648252...

タグ：

posted at 20:33:11

なので、「日韓のどちらが正しいか！」みたいな議論については、「そもそもその問いの建て方無意味ねーよ」と思っている。個々の問題により結論は違うし、議論の前提がナンセンスな場合も多い。

タグ：

posted at 20:35:42

Jiahao Chen 陈家豪 @acidflask

New @TheSIAMNews publication: Algorithms from The Book - all written in @JuliaLanguage #julialang my.siam.org/Store/Product/...

タグ： julialang

posted at 20:37:40

#統計 S市の小学6年生男子の体重の分布を推測するために、n人分の体重のサンプルX_1,…,X_nを得て(i.i.d.の実現値とみなせるとする)、確率モデルp(x|w)と事前分布φ(w)のベイズ統計の方法を使うときには、モデル内で次のようになっていると考えます。

続く

タグ：統計

posted at 20:44:20

#統計続き

(1) wが分布φ(w)に従ってランダムに決めて固定する。
(2) 分布p(x|w)の独立試行でサンプルが生成される。

しかし、現実に得たサンプルは、S市の小学6年生男子全員分の体重のリストからn人分をランダムに抽出しており、上の(1)のようなステップはない！

続く

タグ：統計

posted at 20:44:26

#統計続き

だから、以上のベイズ統計でのモデル内設定は、決して現実とは一致することがない非現実的な想定をしていることになります。

しかし、それにも関わらず、ベイズ統計の方法が現実の問題で役に立つことを示せます。

これは数学的に非自明な話。

タグ：統計

posted at 20:44:27

#統計そういう非自明な数学的な話への理解が不十分なままで、「ベイズ統計いいよ」と主張しても、説得される人はほとんど出て来るはずがないし、出て来るべきでもありません。

ベイズ統計の理解には「ベイズ確率」「主観確率」などの「主義」に関わる概念は有害なだけで無用であり、数学が必要。

タグ：統計

posted at 20:44:28

#統計まるで定説のように語られる事柄なのに、それに説得されると知的にお馬鹿になる方向に進んでしまう話題が、統計学入門では少なくない。

ベイズ統計を「異なる主義・異なる確率概念に基く別の統計学」だと解説しているのを見て、それが正しいと信じた人は知的に劣化してしまったとみなされます。

タグ：統計

posted at 20:48:00

Ninja DAO | CryptoNi @CryptoNlnjaNFT

#統計ふつうに考えれば、そういう怪しい話を信じちゃいけないものだと思う。

あと、もっと入門的な話題では「母集団分布は正規分布であると仮定する」というのを当然に疑う必要がない前提であるかのように誤解するのも知的劣化の1つのパターン。

そこは当然疑うべきで、疑った人だけがまとも。

タグ：統計

posted at 20:50:24

毎日新聞ニュース @mainichijpnews

言い換えると、「日韓どちらが正しい？！」に対して「そもそもその議論意味なくね／議論として成立してないよね」という時もあれば、誰も注目してなくても「これめっちゃ面白いやん」という時もある。それが研究者のお仕事です。

タグ：

posted at 20:52:17

中高生の妊娠相談過去最多　新型コロナによる休校影響か　熊本・慈恵病院 mainichi.jp/20200511/k00/0...

タグ：

posted at 21:14:29

積分定数 @sekibunnteisuu

先生フォロワさんに、駄目なことを教わってしまったようだorz twitter.com/okameclassical...

タグ：

posted at 21:21:32

#統計ところが人間は弱いもので(私もよく騙される)、「みんながそう言っていると、そういうものかなと思ってしまう」ということがあると思います。

しかもそれが単に人が言っているだけではなく、書物に書いてあったりするとさらに「そういうものなのかな」と納得してしまい易くなる。

タグ：統計

posted at 21:26:02

吉田弘幸 @y__hiroyuki

高校の先生，高３生の保護者の方は，文科省に来年の入試をどうするつもりなのか問い合わせてください。
僕も今週中に，予備校講師の立場で問い合わせてみます。
多くの方が声を上げることが意味を持ちます。

タグ：

posted at 21:31:23

Otepipi @Otepipipi

はてなブログに投稿しました #はてなブログ, #Julia言語, #微分方程式
Interact.jlでプロットをgifにする - システムとモデリング
otepipi.hatenablog.com/entry/2020/05/...

タグ： Julia言語はてなブログ微分方程式

posted at 21:32:22

Massimo @Rainmaker1973

Demonstration of momentum conservation [source of the gif: bit.ly/2oSftWc] pic.twitter.com/ph8Mz2x6Ma

タグ：

posted at 21:32:52

The Julia Language @JuliaLanguage

The Julia Language has been selected to participate in Googles Season of Docs! We are excited to hopefully host a few technical writers and welcome them to our community. Read more here: discourse.julialang.org/t/the-julia-la...

#JuliaLang #GoogleSeasonOfDocs #OpenSource

タグ： GoogleSeasonOfDocs JuliaLang OpenSource

posted at 21:42:20

非公開

タグ：

posted at xx:xx:xx

#統計 #超算数

陽性的中率(偽陽性の問題)への「ベイズの定理」の応用について知り、「高校で教えるべきだ」のように言う人はひどく誤解している。

陽性的中率を求める問題は小学生レベルの割合の問題に過ぎず、実際に小学生に問題を出したら解けました。

「算数を復習するべきだ」と言うべき。

タグ：統計超算数

posted at 22:43:31