黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2022年10月22日(土)
うーん…同じ入力に対してpythonの方は2回ずつしかスパイクが立ってないのに、juliaの方は3回ずつスパイクが立ってる…脱分極の速度が違うみたいに見えるなー。そのあたりをヒントにコードを見直すか。。
#ゼロから作るSpikingNeuralNetworks
#スパイキングニューラルネットワーク pic.twitter.com/9xahnNtPa8
タグ: スパイキングニューラルネットワーク ゼロから作るSpikingNeuralNetworks
posted at 00:14:21
#統計 以下では「教科書を読め!」という話にしていますが、手元に学部レベルの入門的教科書がある人は、その入門的教科書では、
信頼区間は推定の章
に
P値と検定については仮説検定の章
に解説が載っていて、
それらの関係について何も説明していない
ことに気付く可能性が高い。続く twitter.com/genkuroki/stat...
タグ: 統計
posted at 03:40:53
#統計 昔からあるそういう傾向のせいで、大学で統計学の講義を担当している人でさえ、区間推定と仮説検定を全然別の話だと信じている場合がかなりあるように思えます。
個人的な意見では、これは信頼区間の実践的な解釈を高等教育を受けても全然できなくなる原因の1つ。続く twitter.com/genkuroki/stat...
タグ: 統計
posted at 03:41:02
#統計 一方、タイトルが「数理統計学」になっている有名な教科書には、信頼区間(もしくは信頼領域への一般化)と検定の双対性(表裏一体性)が標準メニューとして解説されています。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 03:41:05
#統計 以上では、検定側の定式化を実践的統計ユーザーにとって馴染み深いP値で行ったせいで、細部の論理的厳密さが完璧ではなくなっているのですが、検定側の定式化を棄却領域で行えばもっと楽に検定法と信頼区間の計算法の自然な一対一対応が得られます。 pic.twitter.com/lRw5jXTScs
タグ: 統計
posted at 03:41:29
#統計 繰り返し強調していることですが、添付画像のように、検定法と信頼区間を表裏一体のものだと理解することは、単に数学的にシンプルな話であるだけではありません。
検定(P値)と信頼区間の双対性は、P値と信頼区間の実践的な使い方を理解するために役に立つのです。続く pic.twitter.com/n2cXlOy63Q
タグ: 統計
posted at 03:48:45
#統計 P値と信頼区間の表裏一体性を理解していれば、データの数値「n回中k回成功」が与えられたとき、成功確率pをP値pvalue(k|p)に対応させる函数のグラフと、信頼区間を同時プロットしたグラフをイメージできるようになります。
添付画像では点推定の結果も図示しています。 pic.twitter.com/caedmXCBTA
タグ: 統計
posted at 03:48:47
#統計 こういうグラフをイメージできた人は、信頼区間の解釈はP値の解釈に帰着できることも当然だと感じるでしょう。続く pic.twitter.com/sJrsHXET2r
タグ: 統計
posted at 03:50:27
#統計 ここから先に進むためには、次の短い論文がお勧めです。
journals.sagepub.com/doi/10.1177/02...
Discuss practical importance of results based on interval estimates and p-value functions, not only on point estimates and null p-values
Valentin Amrhein and Sander Greenland
2022 pic.twitter.com/43rBzEMsVG
タグ: 統計
posted at 03:54:34
#統計 その論文の共著者の2人は、800人以上の科学者達が統計的有意性に反対したという記事
www.nature.com/articles/d4158...
の3人の共著者のうちの2人です。
この記事の重要キーワードはcompatible, compatibilityです。
上の論文を読めばcompatibilityによるP値と信頼区間の解釈について学べます。
タグ: 統計
posted at 04:00:58
イベルメクチン vs プラセボ JAMA. October 21, 2022.
「軽~中等度のCOVID-19外来患者において、イベルメクチン治療は、プラセボと比較して、回復までの時間を有意に改善しなかった。よって、軽~中等度のCOVID-19患者におけるイベルメクチンの使用を支持しない」 twitter.com/JAMA_current/s...
タグ:
posted at 04:01:17
#統計 以上で紹介したような、統計学の基本的な使用法に関する現代的な議論を読んだ人達は、このスレッドで詳しく説明したP値と信頼区間の表裏一体性に関する直観的な理解が議論の理解に非常に役に立つことに気付くはずです。
journals.sagepub.com/doi/10.1177/02... pic.twitter.com/o5UzF42WVA
タグ: 統計
posted at 04:04:56
#統計 この話を理解できた人にとって面白く感じられることは、一見して全然関係なさそうに思える区間推定法の設定と検定法の設定が、実は
同じことを別の言葉で述べているだけ
であったことです。
そういう話題について「推定と検定は違う」とどうして言えてしまうのか?😅 pic.twitter.com/U0ctuYAKMK
タグ: 統計
posted at 04:18:39
そもそもこの図で説明されないと分からない,というのも何だかな,と思わないでも無い程,自明で,この図で説明されても分からないとしたら検定とかやめた方がいいと思う. twitter.com/genkuroki/stat...
タグ:
posted at 08:08:40
#統計 二項分布モデルのP値と信頼区間の構成結果の具体例については
nbviewer.org/github/genkuro...
比率のP値函数達
に非常に詳しい解説があります。
P値函数のグラフは概ね「とんがり帽子」型もしくはその変種の形になり、信頼区間はそのグラフの高さαの切断で得られる線分になります。 twitter.com/genkuroki/stat... pic.twitter.com/asFNoAUZAZ
タグ: 統計
posted at 11:15:33
#統計 そこでは、4種のP値と信頼区間の組み合わせが定義実装されているのですが、それらはnを大きくするとほぼ一致するようになります。
#Julia言語
↓
nbviewer.org/github/genkuro... pic.twitter.com/v75iJ6gZTQ
posted at 11:18:01
#統計 nbviewer.org/github/genkuro... では ベイズ統計版の信用区間も扱われています。
信用区間も区間推定の一種とみなせるので、自然にP値函数が対応しています。
続く pic.twitter.com/CqQsxE1yMA
タグ: 統計
posted at 11:23:04
#統計 添付画像
①一様事前分布の事後分布とHDI版の95%信用区間
②③④一様事前分布の事後分布とHDI版の95%信用区間に対応するP値函数は近似的にWilsonのP値函数によく一致する。
nbviewer.org/github/genkuro... pic.twitter.com/Q1obUlnKyq
タグ: 統計
posted at 11:23:06
#統計 統計学教育の世界には、
検定と信頼区間の双対性に一切触れない教え方のせいで、それらの実践的な解釈に仕方の理解が難しくなっている
という問題の他に、
非ベイズとベイズ統計の主義思想哲学の違いのみを強調して、実際にどれだけ違うかを適切に説明しようとしない
という問題もあります。 pic.twitter.com/Bq0KmnbFZS
タグ: 統計
posted at 11:27:19
#統計 添付画像の例は二項分布モデルの場合には、通常のWilsonのP値函数と一様事前分布のベイズ統計のP値函数がよく一致していることを示しています。
これは二項分布モデル特有の性質ではなく、シンプルなモデルの多くで成立しています。
この場合には信頼区間とベイズ信用区間はよく一致します。 pic.twitter.com/Vx3LPsXPdd
タグ: 統計
posted at 11:32:15
#統計 今年の5月にイベルメクチンの新型コロナへの効き目が疑わしいという研究の1つ www.nejm.org/doi/full/10.10... が出たのですが、ベイズ信用区間を使っています。
案の定、ベイズ統計を使っているから主義思想哲学が違い、解釈も変えなければいけないというコメントをしている人達が結構いました。続く
タグ: 統計
posted at 11:42:58
#統計 しかし、そこで使われている2つの二項分布の積モデルの一様事前分布のベイズ統計は、数学的には、非ベイズとほぼ同じ結果を与える場合になっています。(3×4=4×3に似ている(笑))
だから、主義思想哲学が違うのように言っても無意味な場合になっています。
詳しくは
↓ twitter.com/genkuroki/stat...
タグ: 統計
posted at 11:42:58
#統計 関連
「{信頼,信用}区間警察」行為は、シンプルなモデルのシンプルなベイズ統計の結果が非ベイズの場合とよく一致していて、互換性がある、という数学的事実を使うことを妨害しているので注意してほしいと思います。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:09:25
#統計 特に、ベイズ信用区間については
❌95%信頼区間には真の値が含まれる確率は95%だと考えてはいけないが、ベイズ統計での95%信用区間であれば真の値が含まれる確率は95%だと考えてよい
とする過剰広告をしても恥じないが、社会的地位は立派な人達が大量にいるので注意した方が良いです。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:16:45
#統計 統計学で主義は無用であることについても、渡辺澄夫さんに従うのが正解だと思います。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:29:18
#統計 Greenlandさんも、Modern Epidemiology, Fourth EditonのChap.23で渡辺澄夫さんと同じように "incorrectly" と説明しています。
It is often said (incorrectly) that “parameters are treated as fixed by the frequentist but as random by the Bayesian.” twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:37:26
#統計 Greenlandさんも「統計学における本当の話を正直にしてくれる人」の1人で、その第23には以下のリンク先添付画像のようなことも書いてあり、非常に面白い。
非ベイズでもベイズでも統計モデルの選択が重要なのですが、モデルの多くが正当化されずに使われている点に注意した方がよい。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:46:59
#統計 非ベイズとベイズの関係については、
xcelab.net/rm/statistical...
Statistical Rethinking
もお勧めです。
これは統計学の世界に伝統的に蔓延る悪しき慣習からの離脱に非常に役に立つ本です。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 12:56:32
#統計 これも関係ある話題。
論文に書いてあったワクチンのefficacyの信頼区間を数値的に再現しようとする試み。
ほぼ成功しているが、一部に微小な数値の違いが生じている。
私の側の誤りなのか、論文側が微小に間違っているのか、今でもわからない。
論文の結論にこの食い違いは関係しません。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 13:11:57
ピンカーせんせの「チーズケーキ」のくだりは,『心の仕組み』下巻(NHKブックス版)の,このあたりですかね. pic.twitter.com/nwjJEOrrOc
タグ:
posted at 14:40:36
#統計 ベータプライム分布(第2種ベータ分布) en.wikipedia.org/wiki/Beta_prim... もベータ分布と同じくらい大事。
ベータ分布は比率pの分布とみなせるが、ベータプライム分布はオッズ比u=p/(1-p)の分布とみなされる。
X~Gamma(a,θ), Y~Gamma(b,θ) (独立)のとき、
X/(X+Y)~Beta(a,b),
X/Y~BetaPrime(a,b). twitter.com/nekokurumihime...
タグ: 統計
posted at 14:51:17
#統計 私も引用ツイートに感想に同感です。
しかも以下のような事実もある。
添付画像はLehmann, Testing Statistical Hypotheses, 1959 (第1版)より。これは仮説検定については鉄板の有名な教科書なので、検定と信頼区間の表裏一体性は遅くとも1959年には常識になっていたことが分かります。 twitter.com/quesokis/statu... pic.twitter.com/irvios9Ouh
タグ: 統計
posted at 16:20:23
#統計 添付画像はLehmann, Testing Statistical Hypotheses, 1959 (第1版)より。
信頼区間には2つの解釈がある。
1. 通常の区間推定
2. 仮説θ=a (aは具体的な数値)の型の検定法の族
こういうことが、信頼区間の実践的に役に立つ解釈の基礎になります。 pic.twitter.com/04a7t5XYUp
タグ: 統計
posted at 16:25:59
#統計 「双対性」という用語について。
「双対性」は数学ではよく使われる用語でdualityの翻訳語で、独特の肯定的なニュアンスがあります。
以下のリンク先の添付画像中にもdualという単語が "a dual interpretation" の形式で登場しています。# twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:40:12
#統計 現代的には、Lehmannさんが "confidence statements" と書いている部分は、"compatibility statements" と積極的に言い変えた方が自信過剰的な解釈を抑制できるので好ましいと言えると思います。
journals.sagepub.com/doi/10.1177/02... pic.twitter.com/sNskZ5cfv1
タグ: 統計
posted at 16:45:17
#統計 尤度函数の側から見ると、P値函数も事後分布も尤度函数の(互いに異なる)正規化。
データの数値xとパラメータ値θについてP値函数の値pvalue(x|θ)はデータの数値xとモデル+パラメータ値θのcompatibilityの指標の1つ。
事後分布の密度函数φ(θ|x)の値も同様にcompatibilityの指標とみなせます。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 17:03:39
#統計 以上で解説した事後分布についての見方は
journals.sagepub.com/doi/10.1177/02...
を読んで理解できた人が、ベイズ統計の事後分布を使うときに役に立つと思います。
非ベイズのP値函数とベイズの事後分布は似たようなものだと考えるとよい。
タグ: 統計
posted at 17:10:39
#統計 以下のリンク先のように、非ベイズのP値函数とベイズの事後分布の間の近似的対応関係を作る(コンピュータ上で実装する)こともできます。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 17:14:56
おお、分かり易い。
#Julia言語 でのgeneric functionについても同じくらい易しく丁寧な解説があるとよいと思いました! twitter.com/y__mattu/statu...
タグ: Julia言語
posted at 17:39:33
Q:今までで一番苦労したことは?
A:ゲーム攻略AIの作成、AIもプログラミングもよく知らない中、独りで本とネットの知識だけで作り上げるのは大変だった
Q:学習の流れは?
A:c→Python→c++→Java→html&css→Dart→TypeScript→Julia
Q:好きな言語は?
A:DartとJulia
タグ:
posted at 18:53:32
今日の結婚式の後「最近zabadakを知って聴き始めて、カラオケでも歌っていて、まさかご本人と会えるとは」と20代女子から話しかけられました。「吉良さんがいたzabadakを知らないのですが」というファンの方が出現してきています。吉良くん、誇らしいだろうなあ
タグ:
posted at 19:33:12
小児科学会の救急市民公開フォーラムはいよいよ明日です!(申込は終了)
Twitterをしていなかったら接点がなかったであろう坂本先生@oshietedoctorやほむほむ先生@ped_allergyと公開講座をできるのは感無量です。
できるだけ分かりやすく話しますので、楽しんでいただければと思います。 pic.twitter.com/pxULnc78kv
タグ:
posted at 20:07:18
明日の私のパートは「風邪と風邪薬」のお話です。
岩崎書店@IWASAKISHOTENさんの担当者様にお願いをして、「ウイルスのはなし」のキャラクターもスライドに登場します😊
amzn.to/3sgOBiK
タグ:
posted at 20:11:21
風邪と風邪薬の話は、こちらの本の内容を中心にお話しします。
医学書で一般の方々には難しいと思いますので、当日はもっとずっと分かりやすくお話しできるよう努力します。
システマティック・レビューとメタ解析で読み解く 小児のかぜの薬のエビデンス
amzn.to/3TGks8k
タグ:
posted at 20:13:33
教員で初任の方が休職・退職してしまうケースが過去と比較して明らかに増えている。
たらればの話をしても仕方ないが、教頭にゆとりがあれば対応できたかも、救えたかもしれないというケースがある。
教頭が忙しすぎる。もう少し校内に目をむけるゆとりが欲しい。
タグ:
posted at 20:44:50
【正方形を見せて「この形は?」という問題に、正方形、長方形、平行四辺形、菱形、四角形のどれで答えても正解になってしまう。】
どれで答えても正解にしない先生は、児童が馬鹿になる方向に誘導しており、虐待していると思います。 #超算数 twitter.com/panterave/stat... pic.twitter.com/OIiYl02sqP
タグ: 超算数
posted at 23:05:39