黒木玄 Gen Kuroki(@genkuroki)/2022年09月11日

ペンローズタイルはde Bruijnによって5次元超立方格子の2次元への射影によって得られることが示された、その論文がこれっぽい
pure.tue.nl/ws/files/43441...

タグ：

posted at 01:48:34

三重積 @triprod1829

実装している方がいる
qiita.com/cometscome_phy...

タグ：

posted at 01:51:09

#統計個人的な意見

統計学入門の教科書によく書いてある比率に関するWaldの信頼区間については、以下の添付画像のようなひどく粗い正規分布近似になっているので、リンク先のWilson版(スコア検定版)との比較で「うげっ」となって欲しいところ。

正規近似の様子が動画のように心の中で見えると楽。 twitter.com/genkuroki/stat... pic.twitter.com/MGxQh4gnsY

タグ：統計

posted at 09:55:45

#統計以上ではWald版の信頼区間の粗さを強調しましたが、「実装の単純さ」との関係でWaldの方法も便利な場合は多いので、Waldの方法も否定的に見ない方がよいです。

二項分布モデル程度なら「Wilsonの方を使え」と言い易いですが、もっと複雑な場合だと言い難い場合もあります。続く

タグ：統計

posted at 10:01:22

#統計続き。例えば、2×2の分割表の場合に、Wilson版の信頼区間を実装してくれているRのパッケージを見つけることに私は失敗しました。あるならば教えて欲しいです。

Wald版は疫学関係のパッケージに見つかった。

タグ：統計

posted at 10:01:22

古いPCを再利用--お勧め軽量Linuxディストリビューション5選 - ZDNet Japan japan.zdnet.com/article/351925... @zdnet_japanより

タグ：

posted at 10:10:41

#統計ガンマ分布とPoisson分布の関係は、確率密度函数のレベルでがんばろうとせずに、累積分布函数のレベルで比較した方が易しいです。

個人的な意見では、ガンマ分布とPoisson分布の離散時間版である二項分布、負の二項分布、ベータ分布の関係の連続時間極限で理解しておくとよいと思います。 pic.twitter.com/3qbwXcNQ5l

タグ：統計

posted at 10:42:15

#統計和で定義されるPoisson分布の累積分布函数をガンマ分布の累積分布函数=正則化不完全ガンマ函数(基本特殊函数ライブラリで効率的に実装済み)で書く公式は、Poisson分布の累積分布函数の効率的な実装を与えます。

これは二項分布とベータ分布の場合と全く同様。

タグ：統計

posted at 10:42:16

#統計 Poisson分布の累積分布函数をガンマ分布の累積分布函数で書く公式は、

　Poisson分布モデルでの片側検定のP値

の

　improper共役事前分布のベイズ統計

でのexactな解釈を与えます。

P値とベイズ統計は「水と油」ではないです。 pic.twitter.com/Vks5wGpog3

タグ：統計

posted at 10:42:19

鰹節猫吉 @sunchanuiguru

#超算数
当然のごとく擁護する馬鹿が出てくるでしょう。
教員が考えている正解は
①1200÷4
②3✕4✕25
のどちらかであると予想されますが、

①か②か、授業でやったほうを書かなければならない

という馬鹿者が出現すると予想していました。 twitter.com/sekibunnteisuu...

タグ：超算数

posted at 11:03:21

#統計再投稿

【論文に記載されたt検定の内容を理解するために必ずしも確率分布の概念は必要ではない】という有害な極論に名を連ねてつつ、t検定について間違った考え方を教え続けているわけです。

こういう人達が統計学の誤用を再生産しているのです。

高等教育の問題の1つがここに見えています。 twitter.com/genkuroki/stat...

タグ：統計

posted at 11:17:37

#統計復習

drive.google.com/file/d/1Qkx5dB...

には【論文に記載されたt検定の内容を理解するために必ずしも確率分布の概念は必要ではない】と書いてあった。

【～日本心理学会～心理統計法標準カリキュラム小委員会
　三浦麻子(委員長)～】続く twitter.com/genkuroki/stat... pic.twitter.com/ljIpZGPtkB

タグ：統計

posted at 11:24:01

#統計心理統計標準カリキュラム小委員会委員長の三浦麻子さん曰く、

【確率分布を教えないと統計神がお怒りになるなら背教しますと言いたくなる】

これ、【論文に記載されたt検定の内容を理解するために必ずしも確率分布の概念は必要ではない】との関係的にまずすぎ。

twitter.com/asarin/status/... pic.twitter.com/JpcAeBEoCo

タグ：統計

posted at 11:31:21

#統計しかも、私のような正当一直線であるがゆえに有益な批判者😊をブロックして恥じない人がそう言っているわけ。次世代の将来が心配だ。

タグ：統計

posted at 11:37:29

#統計難しいことを教えようとすると色々大変になることは、教える仕事をしている人であればみんな知っていて、精神的ダメージをみんな受けている。

人間相手の仕事はそういう精神的ダメージが本当につらい。

でも、そのことが原因で不正直な極論にジャンプしてしまうのは非常にまずいと思いました。 twitter.com/genkuroki/stat...

タグ：統計

posted at 11:37:30

#統計率直に言うと、【論文に記載されたt検定の内容を理解するために必ずしも確率分布の概念は必要ではない】のような不正直な極論に走る人達は、確率分布について十分に理解しているとは思われず、理解度が低いことが原因でどこをどのように誤魔化して教えるのが適切かも全然分かっていないのでは？

タグ：統計

posted at 11:42:22

#統計すでに書いたことですが、t検定で、まず重要なのは、t分布を使った補正ではなく、中心極限定理による正規分布近似の方です。

だから、正規分布からt分布を出す面倒な計算は、実践的にt検定を使う人にとって必須ではないです。

しかし、中心極限定理による近似がいつ有効であるかの知識は重要。

タグ：統計

posted at 11:51:18

#統計そして、昔ながらの伝統的な統計学入門の教育でなされたあの悪しき検定法選択もフローチャートの一部分としてt検定を位置付けているような理解度の低い大学教員は、t検定において中心極限定理の理解が重要であることを何も理解していない可能性が高いです。

タグ：統計

posted at 11:51:21

#統計あのつまらない検定法選択フローチャートのイメージで、それに合わせて、確率分布について解説されても、つまらない話になることは確実で、ヘタをすれば数学が得意な人ほど嫌になっている可能性さえあります。

大学生に「くもわ」「はじき」を教えるのと同じことをしちゃダメです。

タグ：統計

posted at 11:51:22

#統計【教室の空気が濁って淀む】ことを嘆くときには、自分自身の理解度の低さについても嘆くべきだと思います。

数学を教えていると、毎回自分自身の数学の理解度の低さを把握させられます。これほんとの話。

そういうことに気付いたら、勉強し直して、結果をみんなで共有するとよいと思う。

タグ：統計

posted at 11:56:00

ねとらぼ @itm_nlab

お馬さん～入ってますか～？

女の子「トントン、入ってますか～」→馬「入ってますよ～」　女の子と馬の交流がほほえましい【米】
nlab.itmedia.co.jp/nl/articles/22...
@itm_nlabzoo pic.twitter.com/f5MP2b1tF2

タグ：

posted at 12:05:00

#統計

Poisson分布のcdf = 1 - ガンマ分布のcdf

に関する説明の訂正版。青字の部分が訂正箇所です。

こういう公式はよく細かい所をよく間違うのでコンピュータで幾つか例を計算して確認した方が良いです。 twitter.com/genkuroki/stat... pic.twitter.com/hNU9uj7spD

タグ：統計

posted at 12:09:18

#統計 #Julia言語

JuliaのDistributions.jlでは、
パラメータλτのPoisson分布は Poisson(λ*τ) と書け、
パラメータα, 1/τのガンマ分布は Gamma(α, 1/τ) と書ける。

cdf(Poisson(λ*τ), α-1) = ccdf(Gamma(α, 1/τ), λ) の確認
↓
github.com/genkuroki/publ... pic.twitter.com/DutiE9Ckis

タグ： Julia言語統計

posted at 12:40:25

Hideki Kawahara: WAS @hidekikawahara

#統計

 #Julia言語のDistributions.jlパッケージの記号法は、相当に教育的で分かり易い。例：

　cdf(Poisson(λ*τ), α-1)

は

　期待値λτのポアソン分布の
　累積分布函数(cumulative distribution function, cdf)の
　α-1での値

と素直に読める。ほとんどこのまま数式として使いたいレベル。

タグ： Julia言語統計

posted at 12:44:14

Julia v1.8.1が出ていた。

#julialang
julialang.org/downloads/

タグ： julialang

posted at 13:19:17

【教室の空気が濁って淀む】に続く言葉は「講義が下手くそでごめんなさい」であって欲しかった。

【確率分布を教えないと統計神がお怒りになるなら背教します】

【論文に記載されたt検定の内容を理解するために必ずしも確率分布の概念は必要ではない】

はまずすぎ。

drive.google.com/file/d/1Qkx5dB... twitter.com/genkuroki/stat... pic.twitter.com/2YrQ533LMM

タグ：

posted at 13:51:55

Hiroyasu Kamo @kamo_hiroyasu

@a_saitoh 「正方形は長方形でない」も素でそう思い込んでいる小学校教員が多数いるだろうと疑っています。

タグ：

posted at 13:58:36

もなか @monamour555

知念実希人　物語り @MIKITO_777

これまた主語がでかいなぁ。
C を出してくるエンジニアは、善管注意義務を果たしている有能者。
もう少しレベルを落として、言われたとこしかしない(できない)エンジニアとお付き合いすれば、円満解決。大概はコストも安いし弾数も多い。 twitter.com/tmkwsn/status/...

タグ：

posted at 14:03:15

これ……2年前かと思ったら、今日のツイートだぜ……

うそだろ……まだイソジン諦めてなかったのかよ…… twitter.com/hiroyoshimura/...

タグ：

posted at 14:29:38

#統計「確率分布曼荼羅」の類を統計学入門者に勧める傾向も酷いと思います。

「こんなものを勧めて理解できると思っているのか？」と「確率分布曼荼羅」を勧める人を非難する傾向がないのは非常に残念なことだと思います。

悪しき「確認分布曼荼羅」とはこれ
↓
www.math.wm.edu/~leemis/chart/... pic.twitter.com/QTTALgE2Kv

タグ：統計

posted at 14:34:12

#統計こういう「曼荼羅」を見て喜ぶ人達がいることも知っているし、面白がるのもありだと思います。

しかし、これで勉強することを勧めたり、こういう雑過ぎていかにもダメそうなイメージで統計学では使う確率分布について理解することを勧めるのは常軌を逸しています。

www.math.wm.edu/~leemis/chart/... pic.twitter.com/h09H91Wk3w

タグ：統計

posted at 14:38:49

#統計立派な経歴を持つ社会的な信頼度が高そうな人達が、人間が何か難しいことを理解することについて普通に常識的に考えていれば常軌を逸したことを言っていることが明瞭な場合があるので注意が必要です。

後の世代は「あれは酷かった。見習わないようにしよう」となって欲しいです。

タグ：統計

posted at 14:41:43

#統計悪しき「検定法選択のフローチャート」とはこれ。

添付画像は Statistical Rethinking Chap.1 xcelab.net/rmpubs/sr2/sta... より。

Statistical Rethinkingは個人的な意見では普通に常識的な考え方に沿って書かれていて非常によいです。この手のフローチャートも厳しく批判している。 twitter.com/genkuroki/stat... pic.twitter.com/gcupaN4AmQ

タグ：統計

posted at 14:48:53

#Julia言語

自分のパソコンに入れた。パッケージのprecompileちう。

Julia v1.8.1 has been released - Announcements - JuliaLang discourse.julialang.org/t/julia-v1-8-1...

タグ： Julia言語

posted at 15:06:56

#Julia言語 precompileおわた。

404 dependencies successfully precompiled in 489 seconds

タグ： Julia言語

posted at 15:08:34

#Julia言語 Jupyterでもv1.8.1を使えるようにした。

github.com/genkuroki/publ... pic.twitter.com/w82yOTA7Hy

タグ： Julia言語

posted at 15:13:11

Hidehiko MURAO @HidehikoMURAO

Missing…😓 どこへ行ってしまったのでしょう…。
amzn.to/3xbLtb6

タグ：

posted at 16:04:57

Ishida the Brain Dam @tbs_i

昔（MacのCPUがモトローラだった頃）SASが出してたjmpという統計ソフトがこのフローチャートを全面的に取り入れた製品だった。このフローにあるような質問に回答していくといきなりp値が出るという。
どんな分析手法が選択されたかはマニュアル紐解かないと分からなかったはず。 twitter.com/genkuroki/stat...

タグ：

posted at 16:13:37

#統計 P値函数の全体までは行かなくても、「差は無い」という型の仮説のP値だけではなく、信頼区間も報告させることには次のツイートに書くメリットがあると思います。続く

タグ：統計

posted at 16:17:48

#統計続き

2×2の分割表の独立性検定

や

順位和を使うノンパラメトリック検定

でも、信頼区間も報告しなければならないようになれば、それらの検定のユーザーに

違いの指標となるモデルのパラメータが何であるか

を強制的に決めさせることができる。続く

タグ：統計

posted at 16:17:49

#統計

問題1：2×2の独立性のχ²検定を使ったときに、モデルのどのパラメータに関する信頼区間をどのように計算して報告すれば良いのか？

問題2：同様の問題をノンパラメトリック検定について考え、解答例について説明せよ。

こういう問題に答えることが強制されるのは良いことです。続く

タグ：統計

posted at 16:17:50

#統計なぜならば、我々が真に知りたいことは、2つの群が単に「違う」ということそのものではなく、2つの群が「どのようにどれだけ違うか」のはずだからです。

こういう科学的に常識的な思考をできれば、信頼区間も報告させることにもメリットがあることがわかる。

タグ：統計

posted at 16:17:50

#統計問題1について

2×2も分割表の独立性の検定で、「どれだけ独立でないか」を測るためのモデルのパラメータで最も扱い易いものは「オッズ比」です。

常に使えるとは限らないが、その他に「比率の比」や「比率の差」も「どれだけ違うか」に関する指標パラメータとして使えます。

続く twitter.com/genkuroki/stat...

タグ：統計

posted at 16:26:33

#統計続き。しかし、非常に困ったことに、私の探し方がまずかったのかもしれませんが、2×2の分割表の独立性のχ²検定と整合的な、それらの信頼区間を計算してくれるRのパッケージが見つからない。

見つけている人がいたら教えてください。

続く

タグ：統計

posted at 16:26:33

#統計 2×2の分割表の独立性検定には，Fisher検定もあって、Rのfisher.testはオッズ比の信頼区間も表示してくれます。

しかし、fisher.testが表示するP値はSterne型で、信頼区間はClopper-Pearson型なので、それらに整合性はありません。続く

タグ：統計

posted at 16:26:34

#統計続き。報告する数値の間に整合性がないのは、ちょっとよろしくないと思うので、fisher.testのその仕様は極めて不可解です。

その問題はexact2x2パッケージを入れれば解決します。

私も普段は #Julia言語の話しかしていないが、Rについても色々調べています。

タグ： Julia言語統計

posted at 16:29:03

#統計 Rにおいて、

fisher.test が表示するP値と信頼区間には整合性がない。
P値 < 5% なのに、95%信頼区間が1を含む！😱

exact2x2はfisher.testのP値と整合的な信頼区間を表示してくれる。

github.com/genkuroki/publ... pic.twitter.com/RHU9uHQ0dI

タグ：統計

posted at 17:02:24

#統計 Rにおいて、

fisher.test が表示するP値と信頼区間には整合性がない。
P値 < 5% なのに、95%信頼区間が1を含む！😱

exact2x2でtsmethod="central"を使うと、fisher.testの信頼区間と整合的なP値を表示してくれる。

github.com/genkuroki/publ... pic.twitter.com/LBKmj0Q01O

タグ：統計

posted at 17:02:26

#統計 RにおけるepiR::epi.2by2は以下を表示してくれる。

* Wald版のリスク比の信頼区間
* Wald版のオッズ比の信頼区間
* Yatesの連続補正されたχ²検定のP値
* Sterne版のFisher検定のP値

この場合はP値と信頼区間には色々整合性がない。

整合性は気にしない文化？

github.com/genkuroki/publ... pic.twitter.com/oExTADLJoZ

タグ：統計

posted at 17:08:00

#統計個人的な意見は、χ²検定では、Yatesの連続性補正は一切使うべきではないです。

χ²検定でYatesの連続性補正を使うと、Clopper-Pearson版(片側確率の2倍)版のFisher検定をよく近似するようになって、Fisher検定以上にP値が無用に高くなります。

最初からFisher検定を使った方がまし。

タグ：統計

posted at 17:14:59

#統計小サンプルではχ²検定を使ってはいけないという説がまるで標準的な考え方のように流布していますが、その説も間違っています。

コクランルールが正しいと思っている人は、教科書や査読論文の権威に頼らずに、自分で計算して証拠を示そうとすれば誤りに気付くと思います。 twitter.com/genkuroki/stat...

タグ：統計

posted at 17:18:23

#統計補正無しのχ²検定は結構悪くないです。

一般に正規分布近似を使う検定では、Wald法は誤差が大きくなりますが、スコアを使う方法(Wilsonの信頼区間やPearsonのχ²検定の類)は良い性質を持つことが多い感じです。(これは色々な場合を自分で実装して計算してみた経験則。)

タグ：統計

posted at 17:25:19

#統計 Wald法も適切な座標で正規分布近似している場合は全然悪くないです。

タグ：統計

posted at 17:25:19

#統計離散分布モデルで確率を正確に計算できる場合に、そのようにする検定を好むバイアスが結構あるように思うのですが、離散分布モデルでは確率をどんなに正確に計算しても、現れるP値の集合の離散性が原因の誤差が残るので、正規分布近似を使う方法よりも必ずしも優れているわけではないです。

タグ：統計

posted at 17:28:34

#統計そして、その優れているとは限らない側に寄せる補正が、Yatesの連続性補正なので、その必要性に疑問を持つことは自然で、実際、色々計算すると、「だれがこんなものを広めたんだ！手続きを無用に複雑にして、より状況を悪くしているように見える！」となることが多いです。

タグ：統計

posted at 17:31:57

積分定数 @sekibunnteisuu

#超算数
こんなのもあった。 pic.twitter.com/zbqi19OFQB

タグ：超算数

posted at 18:12:45

積分定数 @sekibunnteisuu

#超算数
28,30の意見にほぼ全面的に同意する。

もともとのテーマ、なぜ角度の計算式に「°」を付けるのか、自体は、私も面白いものとは思わないし、どーでもいいが、議論したい人はすればいいと思う。 pic.twitter.com/wmQ8TIcCl5

タグ：超算数

posted at 18:18:04

積分定数 @sekibunnteisuu

#超算数
しかし、曲がりなりにも「数学科　教師　講師」を掲げるコミュの参加者から

「360°＝0°」

＞中学生で「×」の省略を習っているのであれば、答えに「×」があること自体ダメなんじゃないかと思いますけどね。

などの発言が飛び出すのって、冷静に考えたら驚くべきこと（もう驚かないけど）

タグ：超算数

posted at 18:20:03

Julia 使って共同研究したいなあ

タグ：

posted at 19:19:38

準古典Eilenberger理論でRiccati方程式を解く必要があるのだけれど、Julia のパッケージにルンゲクッタがあるので一瞬でコードが書けた。先人のFORTRANコードを解読していた15年前位とは雲泥の差である

タグ：

posted at 19:24:47

しかもルンゲクッタの段や次数を変えたりして精度比較するのも簡単で、びっくり

タグ：

posted at 19:26:51

Home · DifferentialEquations.jl
さまざまな微分方程式を解いてくれる。
Julia のパッケージだけどPythonからも呼べるらしい。 docs.sciml.ai/stable/

タグ：

posted at 19:35:07

清水　団 Dan Shimizu @dannchu

M1 macで Julia1.8.1+Wolfram Engine13.1.0｜清水団 zenn.dev/dannchu/articl... #zenn

タグ： zenn

posted at 19:42:43

清水　団 Dan Shimizu @dannchu

グレブナー基底を調べてみました。 twitter.com/y_minoda/statu... pic.twitter.com/EHkUuQmd2s

タグ：

posted at 19:50:18

#統計信頼区間の報告を要求すると、以下のリンク先のフローチャートが常に使えるとは限らなくなる利点もある？

でも、P値1個とか、区間の両端の値2個とかに、情報を要約して良いという方針だと、簡単に対応されちゃいそうだな。 twitter.com/genkuroki/stat...

タグ：統計

posted at 20:06:08

統計たん @stattan

わかりみしかない。 twitter.com/saitoatsushi/s...

タグ：

posted at 20:19:35

非公開

タグ：

posted at xx:xx:xx

吉藤オリィ@分身 @origamicat

国際的な賞を受賞しオーストリアまで招かれるも、授賞式直前にコロナ陽性になってすぐ近くのホテルからロボットに代行してもらった話｜吉藤オリィ @origamicat #note note.com/ory/n/n5bd196a...

この悔しさはもはや笑い話にしよう

タグ： note

posted at 20:31:39

Hiroo Yamagata @hiyori13

これってプログラミング言語の話のネタだと思って読んでいたら、どうもマジらしくてびっくり。 twitter.com/tmkwsn/status/...

タグ：

posted at 21:04:36

Hiroyasu Kamo @kamo_hiroyasu

Julia日本語記事まとめ @julia_kizi

@ShaqCarson @a_saitoh 「人によっては直感に反することだからきちんと学ぶ必要がある」では納得できませんか？

タグ：

posted at 21:05:03

新しい記事がZennにアップされました！#Julia言語　 #Julia日本語記事
 zenn.dev/dannchu/articl...👈

タグ： Julia日本語記事 Julia言語

posted at 21:07:23

SYC（SY氏） @esuwaishi

Homebrewにも入ったのでアップデートした。
#julia言語 pic.twitter.com/g6EQcfNU6w

タグ： julia言語

posted at 21:10:12

#超算数 #言娯ミヒャエル・シュティーフェル『新算術』(1544)にまたお世話になる。81rページに「名数のかけ算について」という記述archive.org/details/bub_gb...。乗数が名数になることを暗黙に承認している。×記号以前の記述だが参考になる。Abstractus autem numerus, quemlibet denominatum numerum ... twitter.com/temmusu_n/stat... pic.twitter.com/27niENirrB

タグ：言娯超算数

posted at 22:05:26

#超算数 #言娯 multiplicare potest : ut 12 flor. per 10 multiplicati, faciunt 120 flo. Et sic semper recipit summa producta, denominationem, multiplicati uel multiplicantis. Sed 12 flo. per 10 fl. non multiplicantur, nec per 10 gros. &c,
抽象数は名数を乗ずることができる。

タグ：言娯超算数

posted at 22:05:27

#超算数 #言娯 12フロリンに10をかけると120フロリンになるように。積全体が被乗数もしくは乗数の単位(denominatio)を受け取る。しかし12フロリンに10フロリンなどを乗ずることはできない。

このツイの第一の文は12フロリンに10をかけても10に12フロリンをかけても同じといっているようにみえます。

タグ：言娯超算数

posted at 22:05:27

#超算数 #言娯 81r,vは包含除、等分除の話のようだ。

タグ：言娯超算数

posted at 22:12:44

Hiroo Yamagata @hiyori13