黒木玄 Gen Kuroki(@genkuroki)/2019年09月/Page 24

安いエアを乗り継いで行くので荷物重量制限が厳しい。と聞いたので嫁が急遽作った内ポケットカーディガン。パンパンにポケットに詰めて3キロ太って飛行機に乗り込んだらしい。万引き犯か。 pic.twitter.com/3agfinsXI8

タグ：

posted at 23:40:25

#数楽一般に「権威的に響く」もしくは「魅力的に感じられる」可能性のある用語に引っ張られる「理解」を目指すとろくなことがないと思う。

「Kullback-Leibler」「情報量」「大偏差原理」などの用語自体に理解のための重要性はない。

どれも簡単で素朴で一貫した普遍的な考え方で理解できる。

タグ：数楽

posted at 23:40:11

James Tanton @jamestanton

How many essentially distinct ways are there to draw five circles on a page? pic.twitter.com/b6jmH5cIpW

タグ：

posted at 23:30:55

#統計別の議論で正規分布モデルの使用が妥当であることを示すことができない場合には

【100回研究を繰り返した場合、95回以上は真の値～を含んでいる区間が95%信頼区間といえます】

という主張は誤りになります。

おそらく、論文査読者の多くはこの事実を無視しているのではないか？大丈夫？ pic.twitter.com/9ADBSzwKgc

タグ：統計

posted at 23:27:40

#統計正規分布モデルを使って求めた平均の95%信頼区間の定義は「正規分布モデルM(μ)が有意水準5%で棄却されないμの範囲」なので、仮説検定の場合と同じ注意が必要です。

タグ：統計

posted at 23:27:35

#統計続き。もしも本当に棄却したい仮説が

(2) 母集団分布の平均はμである。

であるならば、別の議論で

(1) 母集団分布に正規分布モデルを適用することが妥当。

が正しいことを示さなければいけない。(母集団分布がほぼ正規分布になっていることまで示す必要はない。それは普通は無理)

タグ：統計

posted at 23:21:05

#統計「母集団分布は正規分布でかつその平均はμである」という帰無仮説が棄却される場合には

(1) 母集団分布に正規分布モデルを適用することが妥当。
(2) 母集団分布の平均はμである。

の少なくともどちらかが棄却されることになる。続く

タグ：統計

posted at 23:16:56

#統計個人的な意見では、仮説検定での帰無仮説の中には用いた統計モデルの仮定も入れてしまった方が誤解が少ないと思う。例えば「母集団分布の平均はμである」ではなく、「母集団分布は正規分布でかつその平均はμである」を帰無仮説にしてしまう。続く

タグ：統計

posted at 23:13:43

#統計 #信頼区間続き

信頼区間についても、95%信頼区間を求めるために使った統計モデルの現実の母集団への適用が不適切な場合には、95%という数値が信頼できないものになってしまうことも強調した方が良かったと思う。

タグ：信頼区間統計

posted at 23:03:59

#統計 #信頼区間続き

p値の解説では、現実の母集団に(近似的にも)適用できない統計モデルを使った場合には(具体的には例えば正規母集団の仮定が不適切な場合には)、p値は信頼できない数値になることを強調してある点が非常に良いです。

タグ：信頼区間統計

posted at 23:03:58

#統計 #信頼区間

 www.dr-kid.net/pand95ci
疫学からみたP値と95%信頼区間について
2019-08-31

の非常によろしくない極めて不快な点は内容ではなく、マウスでコピーして引用できないようにしてあること。

批判するための引用に開かれていることはこの手のサイトでは絶対に必要なことだと思う。

タグ：信頼区間統計

posted at 23:03:58

#統計 #信頼区間最近の事例

www.dr-kid.net/pand95ci
疫学からみたP値と95%信頼区間について
2019-08-31

* p値の解説は非常に良い。

* 95%信頼区間の説明も良いが、「使用した統計モデルが正しい」という仮定との関係で、現実には95%という数値の信頼性の問題があることに触れていない。続く

タグ：信頼区間統計

posted at 23:03:53

非公開

タグ：

posted at xx:xx:xx

こなみひでお @konamih

狸穴猫／松村りか @mamiananeko

記事を書いたのですがあまりにも先方にとって刺激的なものになってしまっているので，サーバーを止めろという仮処分を出されると理科探みたいなことになりかねません。公開の仕方を考慮中です。著作権侵害やデマとは無縁なのですけど，彼らは見境がないですからね。どうするかなあ。 twitter.com/konamih/status...

タグ：

posted at 22:45:41

@AS_Insects @genkuroki リンク先のスライドの61~65も興味深いですよ~。
mutsumino.info/ppt/Dr.N_world...

タグ：

posted at 22:22:38

Norio Nakatsuji @norionakatsuji

新帯秀樹 Hideki Shintai @hs_heddy

グーグルが、医学的根拠がない「幹細胞治療」の宣伝を禁止するらしい。日本の新聞や雑誌やTVやネットメディアも見習って欲しい。「認可なき/根拠なき治療」は、病気で苦しみ藁をもすがる気持ちの患者を騙して金儲けをするだけでなく、患者の命を危険にさらす、社会的だけでなく人間として最低の罪悪。 twitter.com/washingtonpost...

タグ：

posted at 21:16:31

How my Julia coding style changed
#Julia言語

 www.juliabloggers.com/how-my-julia-c... @juliabloggersさんから

タグ： Julia言語

posted at 20:22:30

S (ﾂｲｰﾄはｽﾚｯﾄﾞ全体をご確認く @esumii

#数楽 Laplace近似の弱バージョンもlogsumexp関係の話。

log Σ_i exp(n x_i) = n max{x_i} + o(n) as n→∞.

Σ は ∫ でもよい。

大偏差原理もLaplace近似の弱バージョンなので、大偏差原理もlogsumexp関係の話。

logsumexpパターンは普遍的。

タグ：数楽

posted at 19:15:52

中身がないのにプレゼンテーション能力だけ高い、は最悪のパターン。

タグ：

posted at 18:53:03

おばけ @triwave33

コメントいただいてPython (Numba) の方が3倍速いよって言われてマジか！ってなってる。なるべく早く検証する

qiita.com/triwave33/item...

タグ：

posted at 18:24:14

#統計以上に書いたようなことは、

p(x|a,b)=(1-a)N(x)+aN(x-b), N(x)=exp(-x²/2)/√(2π)

の尤度函数をサイズ2000程度のサンプルで作って、プロットしてみれば実感できると思います。うまく計算できれば、正則ではない確率モデルの尤度函数の様子を観察できる。

複数の意味で非常に良い演習問題。

タグ：統計

posted at 17:56:26

#統計 #Julia言語ただし、logsumexpのくそ真面目な実装では最後のlogを取るときのunderflowの問題も気にする必要がある。詳しくは以下のリンク先のスレッドを参照。

Juliaではlog(1+a)の代わりにlog1p(a)を使えばよい。

log(1+a)のaの函数としての直接的実装を使う。

twitter.com/genkuroki/stat...

タグ： Julia言語統計

posted at 17:50:24

#統計欲しい計算が

expを取る → 和もしくは積分 → logを取る

のパターンになっていることはよくある。logsumexpのパターン。和ではなく積分版も結構出て来る。

タグ：統計

posted at 17:44:12

#統計そういう場合には、積分領域K内での対数和の最大値をMとして

g(w) = exp(log f_1(w)+…+log f_n(w) - M)

とおけば、g(w)は真っ当に使える函数になりやすい。そしてf_i(w)の積の積分を

Z = ∫_K g(w) dw + M

で数値計算。これの対数log Zが重要な量になっているというパターンがよくある。

タグ：統計

posted at 17:41:24

#統計正の実数値函数f_i(w)の積

f_1(w)…f_n(w)

ではなく、対数和

log f_1(w) + … + log f_n(w)

なら浮動小数点数の範囲で安全に計算できる。

しかし、f_1(w)…f_n(w)が被積分函数の場合には対数和の exp を取って元の積に戻さなければいけない。そしてexpを取った瞬間に0や∞になったりする。

タグ：統計

posted at 17:36:43

#統計 #Julia言語 Float64では 10^{-383}～10^{384} 程度の正の浮動小数点数を扱えます。log_{10} を取ると -384～384 と非常に狭い！

沢山の正の実数をかけるとこの範囲を容易に飛び出して、浮動小数点数側では0や∞になってしまいます。

これを防ぐには、正の実数の対数を足し合わせればよい。

タグ： Julia言語統計

posted at 17:32:07

#統計確率論がらみの数値計算では「沢山の正の実数を掛け合わせたもの」を計算したくなることがよくあります。尤度函数

(a,b)↦p(X₁|a,b)…p(X_n|a,b)

はまさにそういう例になっている。しかし、これはコンピューター上で標準的に使われている浮動小数点数とは極めて相性が悪い。

タグ：統計

posted at 17:26:18

#統計 #Julia言語ベイズ統計の場合には被積分函数(本質的に尤度函数)をコンピューター上で適切に作ること自体に工夫が必要です。

logsumexpという函数の実装で注意するべきことを知っている人にはお馴染みの問題が生じる。

en.wikipedia.org/wiki/LogSumExp

タグ： Julia言語統計

posted at 17:21:51

#統計 #Julia言語要するに、多変数函数の数値積分を計算したいだけなのですが、それが結構大変だという話。

タグ： Julia言語統計

posted at 17:18:34

#Julia言語 #統計さらに試しに、MCMCで事後分布のサンプルを作らずに、数値積分パッケージ

github.com/stevengj/HCuba...

を使って、(a,b)に関する積分も実行してみましたが、この場合には2変数の積分で済んでいるので結構うまく行きます。

この情報は自分で実験したい人には参考になるかも。

タグ： Julia言語統計

posted at 17:16:52

#Julia言語 #統計 a,bに関する定数函数を除いて、よく呼ばれる対数尤度函数を簡素化した。

DynamicHMC.jl は対数尤度函数+対数事前分布を与えればNUTSでMCMCしてくれるだけの素朴な作りになっているので、こういうことを試すには都合がよいです。

github.com/tpapp/DynamicH...

タグ： Julia言語統計

posted at 17:16:51

#Julia言語のDynamicHMC.jlに混合正規分布

p(x|a,b)=(1-a)N(x)+aN(x-b), N(x)は標準正規分布

で作った尤度と事前分布φ(a,b)=exp(-b^2/20)/√(20π)の積の対数を与えてMCMCさせてみたのだが、尤度部分を

log((1-a)+a*exp(b*x-b^2/2))

に変えたら、サンプルサイズn=2^10で3割程度速くなった。 #統計

タグ： Julia言語統計

posted at 17:16:48

kagshun/EMANON @kagshuntravel

S (ﾂｲｰﾄはｽﾚｯﾄﾞ全体をご確認く @esumii

カバ。ライオンの5倍、サメの50倍人間を殺している獰猛な生き物ですが

このスピードを見てください。
そして何より驚くのは、泳いでるのではなく、水中の底を走っているという事実。 twitter.com/thescarynature...

タグ：

posted at 16:27:07

S (ﾂｲｰﾄはｽﾚｯﾄﾞ全体をご確認く @esumii

@shiba_masa なお、現行の入試が「暗記中心」というのも、少なくとも多くの国立大学については（重大な）事実誤認です。

タグ：

posted at 15:57:30

@shiba_masa 国大協は国立大学の「総意」として明確に反対意見を述べましたが、その発言を座長が遮ってまで強引に推し進めた、というのが事実です。審議会の議事録をご確認ください。web.archive.org/web/2015022415...

タグ：

posted at 15:56:09

栗原裕一郎 @y_kurihara

文芸誌が「お金」特集するっていって呼んで来たのが大澤真幸とか進歩しないねえ。

タグ：

posted at 15:44:36

natsuo @natsuokawai

駆け足でしたがベイズ学習に入門した話をしました。
発表に使った notebook はこちらです。#StudyAIchi

github.com/natsuokw/machi...

タグ： StudyAIchi

posted at 15:20:06

はい、みなさんが静かになるために空気抵抗を設定しました.jl pic.twitter.com/oOol5TdRVB

タグ：

posted at 14:52:13

naoki10007 @naoki10007

@konamih なんでしょう　楽しみですね

タグ：

posted at 14:45:51

広島ベイズ塾 @bayesjuku

【塾生の活動記録】Stanでポン・デ・リングを作ってみた bayesmax.sblo.jp/article/181987...

タグ：

posted at 14:31:15

ボールチェーンがカップから飛び出る様子を再現するためには、カップの底面に接している質点には重力が働かないようにして、初期条件を適切に(←これが問題)すればいいのかな？

手を動かしていないのでまだ感覚がつかめていない。

タグ：

posted at 14:14:16

#Julia言語 #Makie.jl

ソースコードの yossha7.mp4 にうけた！

twitter.com/kb_satou/statu...

タグ： Julia言語 Makie

posted at 14:11:20

計算に対して動画出力が長すぎる……
だが record 以外の動画出力方法がわからん。 pic.twitter.com/caBnyPovfh

タグ：

posted at 14:03:22

分岐

twitter.com/irobutsu/statu...

タグ：

posted at 14:00:50

1つ前のツイートで紹介したボールチェーンもしくはニュートンビーズを数値計算で再現してくれないものかしら。

twitter.com/kb_satou/statu...

タグ：

posted at 14:00:17

重力も入ってます、github にソースを置いたのでご確認いただければ。
twitter.com/ktonegaw/statu...

タグ：

posted at 13:59:51

ごまふあざらし(GomahuAzaras @MathSorcerer

github
github.com/KB-satou/Gotch...

タグ：

posted at 13:58:08

漢字で書いてある日本語の本をひながなにして
大人に渡すとナメてんのかって思われるのに，
数学書を数式を使わずに書いた本を書店に置くと
バカ売れする現象に名前をつけたい．

タグ：

posted at 13:54:59

Lukas Schlögl @LukasSchlogl

Just priceless. An estimated 20% of genetic research papers contain errors because Excel converted some gene names into calendar dates. twitter.com/EconAndrew/sta... pic.twitter.com/ts6sZbFa1C

タグ：

posted at 13:37:45

N=500, 紐っぽい pic.twitter.com/tMX8h5UBR4

タグ：

posted at 13:24:12

N重振り子3D.jl (調子に乗ってN=100)
ちゃんと発散せずに動く。ルンゲ=クッタならね。 pic.twitter.com/MUPmO4hf2E

タグ：

posted at 13:09:38

Andrew Whitby @EconAndrew

This is top shelf trolling, because thanks to Excel “1 in 5” genetics papers contain errors in gene names. www.sciencemag.org/news/2016/08/o... twitter.com/msexcel/status...

タグ：

posted at 12:35:00

非公開

タグ：

posted at xx:xx:xx

飯間浩明 @IIMA_Hiroaki

「『十回』の『十』は『ジュッ』でなく『ジッ』と読む、本来の字音『ジフ』の変化だからだ」という主張があります。歴史的にはそうですが、現在の常用漢字表の「十」には「ジッ」の備考欄に「『ジュッ』とも」と記されています。現在一般的な「ジュッ」がお墨付きを得た格好。テストでは両方マルです。

タグ：

posted at 11:48:39

#Julia言語から #Wolfram言語をMathLinkライブラリ経由で使う方法

まだ表示が数式に整形されないのですが、JuliaからWolfram言語を使う道が開かれたので、今後の発展が期待されると思います。誰かが表示の整形をなんとかしてくれるだけで、相当に実用的になると思う。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/4meE1jlXvh

タグ： Julia言語 Wolfram言語

posted at 11:15:19

かず@猫の下僕 @ns_pooh_

コバ先生にhpvワクチン絶対に娘にも自分にも打たないとリツイートしていた女性。

まだお子さんは5歳のようですので、まだ猶予があります。あと7〜10年くらい、、、
それまでに厚労省が推奨に舵きりをしていただけたらいいなと思います。

タグ：

posted at 11:12:51

ヤマンザ子 @ktykgucci

ルーン@( ｡･ω･)⊃━☆゜*・。 @rune_worth

悲しみの実戦詰将棋です(44金打てば悪いなりにもう少し長引くと思った) pic.twitter.com/yGq5dwMQwr

タグ：

posted at 10:42:56

いずれ来るとは思っていたけどやはり来たか。ウチにはまだ広告来てないけど。結局最後の稼ぎ場所は医療系に落ち着くんだよなぁ…。
NMRパイプテクターの原理を応用した謎の「血液還元装置」 www.gohongi-beauty.jp/blog/28817/ @kuwamitsuosamuさんから

タグ：

posted at 10:15:23

たまさか @TamasakaTomozo

娘さんから、自分に孫ができることをサプライズで知らされた男性のリアクションが、反則級に可愛い。
pic.twitter.com/etoTXxWXlZ

タグ：

posted at 09:56:14

おゆ#世界中の子どもたちが安穏に♨ @oniyuri_

#統計指数型分布族にのみ適用できる方法を見たら、「特別に簡単な場合にのみ適用できる方法であり、より複雑なモデルには適用できない」と判断するべき。

指数型分布族はかなり一般的な確率分布族であるというのはひどい誤解。

twitter.com/genkuroki/stat...

タグ：統計

posted at 09:47:09

@konamih わくわく♪ (*´ω｀*)

タグ：

posted at 08:30:56

twitter.com/konamih/status...

楽しみ❤️

タグ：

posted at 08:13:28

こなみひでお @konamih

山田太郎 ⋈（参議院議員・全国比例） @yamadataro43

今日も日本システム企画の熊野社長がカンカンになるようなネタを投下しますので，みなさんお楽しみに。火に油を注ぎ返してやるのだ。

タグ：

posted at 07:58:36

35人が亡くなり多くの負傷者も出た京都アニメーション放火事件。25億円を超える寄付金への税制優遇を決めました。当初から粘り強く働き掛けてきた結果です。これが京アニ事件だけでなく放火事件で苦しむ他の事案についても税制優遇の道を切り開く前例になればと考えています。www3.nhk.or.jp/news/html/2019...

タグ：

posted at 07:36:30