黒木玄 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月31日(月)
#統計 McShane-Gelman 2022 stat.columbia.edu/~gelman/resear... のpp.2-3のIllustrationの節のFig.1を3つの図に分解を再改訂。
帰無仮説「効果=0」のP値をnull P値と書き直した。 pic.twitter.com/7cjBaR9UWe
タグ： 統計
posted at 00:07:30
#統計 特にこの3枚の図は、null P値、点推定、区間推定の3点セットの解釈の仕方を知る上で特に有益だと思います。 twitter.com/genkuroki/stat... pic.twitter.com/We4aJCjcL8
タグ： 統計
posted at 00:13:58
#統計 P値側から定義された信頼区間の易しい例については
github.com/genkuroki/publ...
を参照。WaldとWilsonの信頼区間の2つを理解すればひとまずは十分だと思う。
タグ： 統計
posted at 00:17:42
#統計 あと悪しき信頼区間警察はベイズ版の信用区間警察でもあることが多い。非ベイズとベイズの違いを強調してはいけない場合に強調するのはやめて欲しいと思います。
その辺については以下のリンク先の長大なスレッドを参照。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:22:56
#統計 事前分布を使うベイズ統計であっても、nを大きくすれば、尤度函数の部分が支配的になるので、シンプルなモデルでは事前分布を使わない非ベイズの場合と本質的に同じになります(注意: 一応数学的には自明でない部分があるので注意！)。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:25:07
#統計 添付画像は、二項分布モデルの通常のP値函数(通常の信頼区間を与える)と一様事前分布のベイズ統計版のP値函数(ベイズ版信用区間を与える)の同時プロットです。
この程度の違いしかないのに、違うことを強調するのはひどくミスリーディングです。
github.com/genkuroki/publ... pic.twitter.com/4yBVRYmAti
タグ： 統計
posted at 00:31:13
#統計 P値函数やベイズ統計での事後分布と点推定と区間推定の関係。こういう感じで、非ベイズとベイズを同じように気軽に解釈して大丈夫。
どうしてこういう気軽さを崩そうとするのか？
おそらく、万人が有用なツールを使えるようになることより、自分の主義の布教の方が優先順位が高いのだろう。 twitter.com/genkuroki/stat... pic.twitter.com/N3JJmdk54f
タグ： 統計
posted at 00:45:06
びっくりしたところに赤で印をつけてみた。
これかなりまずいですね。 twitter.com/caazzy/status/... pic.twitter.com/Tx39HNlEk0
タグ：
posted at 01:02:25
akima@4th vaccinated @akima9936
www.mext.go.jp/content/202209...
↑
令和4年9月5日版はこちら
twitter.com/genkuroki/stat...
タグ：
posted at 01:13:27
akima@4th vaccinated @akima9936
黒木さんと同じく「うわっ」という反応しててなんかワロタ
twitter.com/akima9936/stat...
タグ：
posted at 01:14:44
宮澤大輔氏が息子をくも膜下出血で亡くした宮坂教授に対してワクチンで脳卒中起こす、息子に嘘ついたなと口さがない投稿をしてますが
くも膜下出血は超大規模調査で接種後の発症集積も見られず(図3),コロナ流行が限定的だった昨年の日本ではくも膜下出血死亡動態の異常増加は一ミリも見られません(図4 pic.twitter.com/ME3zz94TPG
タグ：
posted at 12:39:07
@Joe_medi3 @sekkai なぜかJoeさんは何が問題であるかに一切言及していない。
検定法Tと整合的な効果の慣習的指標がないという理由で指標θを使うこと自体は確かに問題ないです。
しかし、θの信頼区間が0をまたいでいるのに「効果無しを意味する仮説θ=0は棄却されなかった」と正直に言わないことには問題があります。
タグ：
posted at 13:23:54
@Joe_medi3 @sekkai 続き。
「効果無しを意味する仮説θ=0は棄却されなかった」と正直に結果報告に書かないだけではなく、効果の指標の統計的有意性を判定できない検定法TのP値が5%を少し切ったことを強調して「効果は統計的に有意」と強調することにも問題があります。
続く
タグ：
posted at 13:27:45
@Joe_medi3 @sekkai 続き
さらに、もしもθについて【信頼区間までは出さないことが多い】という慣習があって、θの点推定値だけを報告して、点推定の精度をわからなくすることが許されているとしたら、統計学入門レベルでアウトです。
普通の常識は「パラメータの推定では点推定値だけを報告してはいけない」です。
続く
タグ：
posted at 13:32:33
@Joe_medi3 @sekkai こういうことが許されているというのが結論でしょうか？
↓
検定法Tと整合的な効果の慣習的指標がないという理由で指標θを公式に採用した。しかし、θの信頼区間が0をまたいでしまったので、信頼区間は報告しないことにした。信頼区間は報告しなくてよいことになっている。😈
これは非常にまずい。
タグ：
posted at 13:37:19
@Joe_medi3 @sekkai あと、Joeさんのそのアカウントですが、2009年作成のアカウントなのに、この話題に関するリプライ以外の発言が見当たりません。
わざわざそういうアカウントを「これ専用」に使ってコメントする行為を私が不快に感じたことも付け加えておきます。 pic.twitter.com/X2krnV4Dx3
タグ：
posted at 13:41:39
@Joe_medi3 @sekkai 私の側からは、Joeさんは「これ専用」として採用したアカウントから、何が問題であるかに一切触れないように注意しながらコメントしているように見えています。
こういう話題ではそういう怪しげな行為は避けて欲しいと思います。 pic.twitter.com/ru8RQkePUc
タグ：
posted at 13:45:38
#統計 補足。データの数値xに区間[L(x),U(x)]を対応させる適当な条件を満たす任意の区間推定法が与えられているとき、仮説θ=θ₀の棄却域を{x|θ₀∉[L(x),U(x)]}と定義することによって逆に検定法が得られます。
だから、θの区間推定法と仮説θ=θ₀の検定法(θ₀は自由に動かせる)は一対一対応している。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:01:47
#統計 95%信頼区間の意味での95%を我々は確実さの尺度として使いたいのですが、実際には、信頼区間の構成に使った統計モデルの妥当性の程度によって、95%という数値自体が現実においては信頼できなくなる場合が結構あります。
特に統計モデルが複雑だとそうなり易い。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:40:52
#統計 しかも、シンプルなモデルとおとなしめの事前分布の場合には、信頼区間とベイズ版信用区間の間に実践的に意味のある違いがないことがよくあり、そうなりがちな数学的理由もあります。
信頼区間とベイズ版信用区間の違いがほぼない場合にまで95%の意味の違いを強調しても実践的には無意味です。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:46:20
#統計 信頼区間はP値函数から得られ、ベイズ版信用区間は事後分布から得られます。P値函数の値には
統計モデルのパラメータ値のデータの数値との相性の良さ(compatibility, 両立性)
という解釈があるのですが、ベイズ統計での事後分布も(統計モデルに事前分布を含めれば)同じように解釈できます。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:50:58
#統計 繰り返しになりますが、これって、「今まで普通に許されて来たことをいきなり否定されては困る」という話になっているんじゃないか？
別スレ
↓ twitter.com/genkuroki/stat...
タグ： 統計
posted at 16:00:32
#統計 関連スレッド
以下のスレッドでは、モデルがシンプルな場合のベイズファクターのP値への翻訳の方法も紹介しています。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 16:08:34
うちの小学校は式にも単位や助数詞を強いているんだけれど、この÷２のところの『つ』はさすがにヤバくない？（書かないと×でやり直し）
こういうのも ＃超算数 ってことでいいのかなw pic.twitter.com/nIq2b90I7l
タグ： 超算数
posted at 16:40:50
#統計 非常に素直で自然な素晴らしい疑問だと思いました。
一般的に通用する粗い結果(本質的に渡辺澄夫『ベイズ統計の理論と方法』第3章に書いてある正則モデルの漸近論)で満足できないなら、ケースバイケースで地道に手とコンピュータによる計算する必要があります。
地道さ、大事。 twitter.com/e0_rainy/statu...
タグ： 統計
posted at 17:29:33
#超算数 しばしば【算数の時間に先生が言った内容をどれくらい覚えているかをはかる】ことに価値があるという主張が見られるのですが、算数の学力と無関係な【先生が言った内容】にはどれだけの価値があるのでしょうか。
先生の好きな食べ物は？ 答え、カレー
というのとどう違うのでしょうか。 twitter.com/rio__tan/statu...
タグ： 超算数
posted at 17:43:00
#統計 そういう自然な疑問から得られる知識はワクチンの治験の結果の報告を読むときにも役に立ちます。
現代人的には結構重要な教養でしょう。
実際、以下の論文に目を通した人は意外なほど多いと思われます。推定にシンプルなベイズ統計の方法を使っている。続く
www.nejm.org/doi/10.1056/NE...
タグ： 統計
posted at 17:54:46
#統計 添付画像は www.nejm.org/doi/10.1056/NE... より。
Table 2(添付画像②)には95%ベイズ信用区間が書いてあります。
親切なことに非ベイズの通常の95%信頼区間もTable 3(添付画像③)に書いてあります。続く pic.twitter.com/bogz1WBqeb
タグ： 統計
posted at 17:54:49
#統計 続き。しかも、Table 3に書いてある95%信頼区間はClopper-Pearsonの信頼区間で保守的に幅が広くなりがちな信頼区間になっています。
しかし、そういう保守的な信頼区間を使っても、ベイズ信用区間から見て取れる結果と実践的に意味のある違いはないことも分かります。 pic.twitter.com/HAooFn5Aui
タグ： 統計
posted at 17:54:52
@genkuroki 理由は、たぶん出版のページ数制約のが強いと思います。
容量制限の緩いsupplementにグラフと生数字を載せるようになるといいなとは思います。
タグ：
posted at 18:10:12
#超算数 自分ちの子が
❌「ずつ」が出て来たらかけ算の問題なんだよ。
とか
❌「しき」では「ずつ」のついた数を先に書くんだよ。
とか言うようになっていたら、子どもが相当にまずい状況に陥っていることになるので、きちんと説明して助けてあげる必要があります。
算数教育は結構怖い。 twitter.com/kazumarugou/st...
タグ： 超算数
posted at 18:24:52
#羽生善治 九段の勝利
【第72期ALSOK杯 #王将戦 挑戦者決定リーグ】
●永瀬拓矢王座 VS ○羽生善治九段
第72期ALSOK杯王将戦挑戦者決定リーグ戦、永瀬拓矢王座VS羽生善治九段戦が10月31日（月）に行われ、羽生九段が94手で永瀬王座に勝ちました。
buff.ly/3DLp9sd
posted at 19:05:56
@kosugitti @KunisatoY Quarto内で任意の言語を動かすのは特殊なQuarto組み込み機能を除いてJupyterエンジンとKnitrエンジンが担っているのですが、Jupyterは1ファイル1カーネルの原則があるため1ファイル内に多言語を共存させるのはknitrエンジンを選択する必要があります。要するにRMarkdownと同じです。
タグ：
posted at 19:20:52
@kosugitti @KunisatoY そしてknitrがPythonを実行するにはreticulate、Juliaを実行するにはJuliaCallというRパッケージが必要です。（knitrのドキュメントに書いてあります）
逆に言えばこれらを入れるだけで動きます。例はこちら。
github.com/quarto-dev/qua...
タグ：
posted at 19:23:40
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
佐伯胖氏のスウキョウキョウ批判は鋭い。
最初の見開きと最後の見開きだけ読んでも十分、本質をえぐっている twitter.com/OokuboTact/sta...
タグ： 超算数
posted at 20:41:06
OokuboTact 大久保中二病中年 @OokuboTact
スウキョウキョウの理論には全く感心しないけれど、理系の専門家達が算数教育を改革しようと文部省と戦っていた時期があった。
今では同じ穴のムジナという印象しかない
＃超算数 twitter.com/OokuboTact/sta...
タグ： 超算数
posted at 20:43:25
#超算数 家庭教師派遣会社の虎の巻が【
こういう問題で「ずつ」という言葉が出たら
手っ取り早く「ずつ」の数が先〜！
解答パターンの決まった式として”考えないように教えなさい”】
と指導している場合があるのか！ twitter.com/kazumarugou/st...
タグ： 超算数
posted at 20:49:03
@bob3bob3 @kosugitti @KunisatoY knitrなのは同じなので記法も同じですよ。私がそのGitHubに挙げた例はコードブロックを跨いでPythonとJuliaがそれぞれxを保持していることを示しているだけでPythonとJuliaの間は跨いでいないです。
これがJupyterのマジックコマンドだとコードブロックを跨げないので同一言語でもxを保持しません。
タグ：
posted at 22:16:40
www.buzzfeed.com/jp/naokoiwanag...
【厚労省「小中高校担当とは連携しているのだが...」文科省「HPVワクチンの位置付けが変わったことを知らなかった」】
【位置付けが変わった】？？？
時期によらず、反ワクチンの言説を学生に吹き込めという内容の通知を大学に出す行為は相当にひどい。 twitter.com/bfjmedical/sta...
タグ：
posted at 22:27:33
#超算数 【学習指導要領の「場面を式に表したり，式を読み取ったりすること」の部分を常識的な解釈しか許さない別の表現に改定し、その改定理由も公式に公開する】という提案は結構有効だと思います。 twitter.com/genkuroki/stat...
タグ： 超算数
posted at 22:54:24
新しい記事がQiitaにアップされました！#Julia言語 #Julia日本語記事
qiita.com/MandoNarin/ite...👈
タグ： Julia日本語記事 Julia言語
posted at 22:57:23
#超算数 かけ算順序固定強制指導を批判している人まで、「かけ算順序指導はできない子を助けるための教育的配慮である」という主張をある程度受け入れていたりする。とても残念。
以下のリンク先で引用されている添付画像の資料を見て、できない子への教育的配慮に見えるとしたら、相当におかしい。 twitter.com/genkuroki/stat... pic.twitter.com/T6FXKSNbTm
タグ： 超算数
posted at 23:08:34
#超算数 添付画像の資料を見ればわかるように、かけ算順序固定強制指導は「子どもに間違いをただす」ために行われるようになったのであり、できない子への配慮から始まったのではない。 twitter.com/genkuroki/stat... pic.twitter.com/nHKlmTPgQ8
タグ： 超算数
posted at 23:49:34
