黒木玄 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 を見て下さい。
2020年12月25日(金)
訂正： NetaUtils.jl じゃなくて、MetaUtils.jl です。
@ teval (:call, :sin, (:call, :/, π, 6))
→ 0.49999999999999994
に注目。tupleで書かれた式をevalできている。 #Julia言語
nbviewer.jupyter.org/github/genkuro... twitter.com/genkuroki/stat... pic.twitter.com/xGgBgIrhYr
タグ： Julia言語
posted at 23:57:18
尤度の正確な定義は「モデル内確率分布で測った観測データと同じ数値が生成される確率またはその密度」です。この定義およびそこから導かれる数学的性質を超える使い方は単純に杜撰な使い方とみなされます。
尤度は「もっともらしさ」でも「証拠」(の度合い)でもありません。
タグ：
posted at 23:25:03
添付画像は
web.archive.org/web/2018052405...
からの引用。松王さんはソーバーさんの本の翻訳者。
統計学における「尤度」を「証拠」とみなす行為は馬鹿げています。尤度は決してもっともらしさでない。
これを読んで大変な悪印象を持ちました。
どうしてこういう言説を周囲の人達は批判しないのか？ pic.twitter.com/5coTFZ2oPJ
タグ：
posted at 23:22:29
ソーバー『科学と証拠 統計の哲学入門』(松王政浩訳)が、尤度函数という情報を大量に捨てる要約の仕方に「主義」の力で特別な価値を与える、とか、尤度原理の類で仮説検定を批判する、というような論外な行為をしていないならば購入してもよいと思っているので、読んだ人がいれば教えて下さい。
タグ：
posted at 23:14:50
あと私は読んでいないソーバー『科学と証拠 統計の哲学入門』(松王政浩訳)の「尤度主義」の節に尤度原理にあたることが多分書いてあると思います。
尤度函数という情報を大量に捨てる要約の仕方に「主義」の力で特別な価値を与えている本に見えたので読んでいません。
www.unp.or.jp/ISBN/ISBN978-4...
タグ：
posted at 23:12:21
生活保護は現時点で、経験したことのない人が想像する以上に厳格です。本当に、還付金ですら受け取れないし、たった1円でも口座に振り込まれていたら次の月は給付がストップになる。そのくらいもうすでに厳しいものです。それでも、最低ラインの生きる権利であり、助かる人がいます。
タグ：
posted at 23:10:28
「尤度原理」の解説は特殊な文献にしか載っていないので、普通に手堅く統計学を勉強して来た人達の多くは知らないと思う。
載っているのは、豊田『瀕死の統計学を救え！』とか大塚『統計学を哲学する』とか、美添「統計と経済の間で」 www.terrapub.co.jp/journals/jjssj... など。
タグ：
posted at 23:03:35
私自身20代で非常に困窮し、癌になった母を生活保護で助けてもらった期間があります。幸い『戦場のコックたち』が直木賞の候補になり、ヒットしたおかげで生活保護を抜け元の生活に戻せたのですが、そんなラッキーなこと滅多にないでしょう。せめて、誰もが生活を送れることを。#生活保護は権利です
タグ： 生活保護は権利です
posted at 23:02:40
岩永さん自身の記事にはいいものが多い(全部ではない)けれども、「マスメディアは権力である」という自覚が足りないから、朝日新聞ではなく匿名医師を非難するのですよ。まずはメディア関係者自身が「権力を持つものとしての謙虚さ」を身につけなければ、信頼されない
タグ：
posted at 22:49:55
ごまふあざらし(GomahuAzaras @MathSorcerer
ちなみにその自分のコードは
Julia 標準ライブラリの REPL モジュールなどから拝借しました． twitter.com/abap34/status/...
タグ：
posted at 22:48:30
メディアに携わるものとしての反省が足りない。それを言うくらいならもっとメディア批判(非難)を展開するべき。例えば、火元が朝日新聞であることはもう世界中に知られているのだから、朝日新聞に取材して経緯を明らかにするなどを試みればいい。
メディアは単に身内に甘い twitter.com/nonbeepanda/st...
タグ：
posted at 22:45:32
#統計 主観確率のアドホックな公理からの導出とか、結果的に巧妙な詭弁になってしまっている「尤度原理」の「証明」を基礎に置いてしまった20世紀の主観主義ベイジアン達は、「あれもこれも主義に従わなくても適切に使えば合理的に役に立つ」という穏健で合理的な立場に立てなくなった。
日本の事例↓ twitter.com/genkuroki/stat...
タグ： 統計
posted at 22:34:22
#統計
phil-stat-wars.com/2020/12/03/jan...
January 7,～**note time modification and date change
Putting the Brakes on the Breakthrough,
or “How I used simple logic to uncover a flaw in a controversial 60-year old ‘theorem’ in statistical foundations”
Deborah G. Mayo
タグ： 統計
posted at 22:27:48
非公開
タグ：
posted at xx:xx:xx
A particle system with perceptually uniform colormaps
github.com/keijiro/Pugrad pic.twitter.com/IRzWOH3xiD
タグ：
posted at 22:18:32
#Julia言語
Juliaでのパイプは
github.com/jkrumbiegel/Ch...
Chain.jl
を使えばかなり便利だと思います。
マクロで何とかできるなら、アイデアを出せば誰か実現してくれるかも。 twitter.com/aivriahrgftt/s...
タグ： Julia言語
posted at 22:01:28
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
@bampaku ＞ 流行語大賞は「指導要領(意味：指導要領解説のこと)」ですかね。
「指導要領に書いてあります！」
＃超算数 ＃書いてない
＃学習指導要領解説は学習指導要領ではありません
タグ： 学習指導要領解説は学習指導要領ではありません 書いてない 超算数
posted at 21:22:48
無症状の時点で最も高い感染性を持つウイルスなので、弱毒化への淘汰圧はあまり高くないと思うのですよね、そもそも twitter.com/shinyamatsuura...
タグ：
posted at 21:21:13
matlabを手に入れます
↓
無料ライブラリが不足してpythonに行きます
↓
速度が足りなくなってc/c++に帰化します
↓
Juliaの誘惑 twitter.com/actinouran/sta...
タグ：
posted at 20:53:12
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
最強の将棋ソフト「やねうら王」からクリスマスプレゼント ～詰将棋500万問が無償公開／3手詰めから11手詰めまでの詰将棋をそれぞれ100万問ずつ、SFEN形式で forest.watch.impress.co.jp/docs/news/1297... pic.twitter.com/cn8bNFgk4a
タグ：
posted at 17:37:42
『社会科学のためのベイズ統計モデリング』、第7章まで読了。
第7章の最後に、「本章で解説したモデル評価の方法は、あくまで「ある限定的な観点からみた場合の指標」であることです。これらの指標で優れているモデルが、あらゆる観点でよいモデルであることを意味しません。」と記述されていた。
タグ：
posted at 17:01:25
「従ったフリをしておけ、でいいんじゃない？」から「それは意外に害が大きいんじゃない？」というのが、ここ1～2年で私が考えを改めた点ではありますな。「堂々とバツをもらって来い」が良いと思います。
タグ：
posted at 16:55:18
教員のなり手がいなくなる主たる原因は、
・給料が安すぎる。
・待遇が悪すぎる。
・求められることが常軌を逸している。
・職場環境が劣悪すぎる。
・プライベートを犠牲にせざるを得ない労働環境。
です。労働者として最低限守られるべき事が守られていないから。ただそれだけです。
タグ：
posted at 15:40:00
#統計 伝統的な主観主義ベイジアン(SavageやらBirnbaumが神様らしい)の偏狭で攻撃的な態度と、以下のリンク先スレッドで紹介したGelman-Shalizi(2013)がベイズ統計においても仮説検定❗️やグラフによる視覚的確認などのあらゆる手段を使う合理的な態度を比較すれば、どちらがまともであるかは一目瞭然。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:29:01
#統計 不合理な教条主義者でもなければ当然とみなせない前提を当然とした上で、一方的で非常に攻撃的な態度になっています。
「～という考え方もある」「同一の数学的道具を違う使い方をしてもよい」「あれもこれも有用である」のような穏健な態度を取っていないのでひどく有害。 pic.twitter.com/YzxZ832bhJ
タグ： 統計
posted at 15:18:58
現代のトップアスリートにとって、「数学的・物理学的アプローチでさらに上を目指す」というのはもはや常識。 twitter.com/vecchio_ciao/s...
タグ：
posted at 14:50:49
OKUMURA, Akira（奥村 曉） @AkiraOkumura
VS Code に Jupyter の plugin 入れて課題をやってる学生とかもいて 2020 年というのは素晴らしいですね。
タグ：
posted at 14:29:35
中村 良幸 (Nakamura Yosh @nakayoshix
JavaScriptからTypeScriptは自然な流れだし、Pythonから入ってJuliaへ流れるのも悪くないかなと。（人によってはGoへ流れるんだと思いますが、科学技術計算の文脈においては…という話です）
というわけで、「脱初心者した段階で他の言語に流出する」のはむしろ良いことじゃないかと思った次第です。
タグ：
posted at 14:20:57
掛算順序問題との類似
* 掛算順序が逆なら誤りだと教えている人達は論理的な考え方を教えていると信じていたりする。
* 主観主義ベイジアンは主観主義ベイズ統計学は公理から導かれる基礎を持つ論理的に整合的な理論だが、頻度主義統計学はそうではないと信じている。
😱😱😱😱😱😱
タグ：
posted at 13:53:10
カニは皆さん食べたことあればあのスジが腱です。外骨格とつながってます。並列に、斜めにたくさん配列することで力の効率という観点で我々なんかよりエレガントです。 twitter.com/katzkagaya/sta...
タグ：
posted at 13:36:23
#統計 美添泰人さんが主観主義ベイジアンの立場でベイズ統計について語っている文献を見つけたので紹介。
www.terrapub.co.jp/journals/jjssj...
経済と統計の間で
添付画像に引用した部分の「標本理論」を「頻度主義」に置き換えれば、主義によって何が定数になるかが変わるといういつものダメな話になります。 pic.twitter.com/TMokXoSgZY
タグ： 統計
posted at 13:35:04
生きているときに、、、なんてとても想像したくない。危険すぎる。我々脊椎動物よりも機械効率めちゃくちゃ高い設計になってますから。ぎっしりと並列化もしている。ボディビルがんばっても行き場のない筋肉が盛り上がるだけで並列化できないですが、これはそうではない。
twitter.com/K0TSU_K0TSU/st...
タグ：
posted at 13:17:55
『「空気が変わり始めてからものを言い始めた匿名の医療者がメディアを口汚なく批判」と前提条件』
HPVワクチンについて
『空気が変わり始める前から』
『実名』
という『前提条件』ならメディア批判する権利がある？
メディア側で医療者を『区分け』するんですか？
twitter.com/nonbeepanda/st...
タグ：
posted at 13:12:47
@jazzmas62 @maseguchi 以下のリンク先のスレッドでは、Gelman-Shalizin (2013)をネタに(情報量規準にこだわらない)「なんでもあり」のモデルのチェックについて詳しく説明しています↓ twitter.com/genkuroki/stat...
タグ：
posted at 13:06:01
@jazzmas62 @maseguchi 実際には、私以上に、情報量規準だけで判断することにどのようなリスクがあるかを、数値計算とグラフを使って解説をしている人を見たことがない。
真の予測誤差とAICやWAICの(差の)逆相関は私が繰り返し解説している話題の定番。(「逆相関」は書き間違いではない。漸近的に相関係数が-1になる！)
タグ：
posted at 13:02:44
まあ『口汚い批判』というのは、こういう事かと↓
『マッチョめ、滅びろ。男性ホルモンを名前にするアカウント、フォローしている人ブロックします』
『大嫌いだ。筋トレしてる奴、健康に気を使っているやつ』
twitter.com/io302/status/1...
タグ：
posted at 13:00:51
『安全地帯から正義の石を投げるのは楽しいでしょうね』とメディア側の人間が、医療者を蔑むのであれば、さらに溝は深まるばかり。
我々医療者は現場でHPVワクチン対応の矢面に立っており、『安全地帯』にいると思ったことは一度もない。
タグ：
posted at 12:55:48
@jazzmas62 @maseguchi 私がまるで「情報量基準だけで見ろ」と言っているかのように誤解されるような発言をするお馬鹿さんを、きちんとお馬鹿さん扱いしてくれる人達がいると個人的には非常に助かります。
私自身がいちいちその手のお馬鹿さんの相手をするのは正直辛すぎ。
タグ：
posted at 12:54:34
『しかし、空気が変わり始めてからものを言い始めた匿名の医療者がメディアを口汚なく批判するのも、安全地帯から正義の石を投げるのは楽しいでしょうね』
空気が変わり始める前から、HPVワクチンのツイ発信してた匿名の医療者から言わせてもらうと、メディアも医療も罪は重い。だから叩かれてる訳で twitter.com/nonbeepanda/st...
タグ：
posted at 12:54:15
@jazzmas62 @maseguchi 【黒木さんは情報量基準だけで見ろと言ってるわけではなくて～】
私が「情報量基準だけで見ろ」と言っているかのように誤解されるような発言をした人は、私に対してきちんと謝罪して欲しいと思います。
タグ：
posted at 12:50:17
#統計 関連
#Julia言語 の確率分布のパッケージ Distributions.jl では、最尤法が実装されていない場合に最尤法のメソッドを追加する方法↓ twitter.com/genkuroki/stat...
posted at 12:40:23
@genkuroki 字数に収まらないので一部だけですが。
@maseguchi: データ分析の専門家であっても、自分の分野の範囲の統計学しか知らないことがほとんどだと思うので、情報量基準示すだけでモデルの妥当性を説得される状況なんて見たことなくても、
タグ：
posted at 12:39:52
@jazzmas62 私がしている統計学の話題では、実際には仮説検定や信頼区間のような真に入門的な話題が多いです。
数値実験でもグラフを描いての視覚的確認を重視。実践的にもそうするのが当たり前。
どこをどう見れば「情報量基準だけ」とか言われることになるんだ？
タグ：
posted at 12:31:42
@jazzmas62 【黒木さんは情報量基準だけで見ろと言ってるわけではなくて～】
え？そういうことを言っているお馬鹿さん(達？)がいるの？
Gelman-Shalizi (2013)で、ベイズ統計においても、仮説検定やグラフを描いて確認することが大事(これは大昔から言われている常識)だと強調している話も紹介しているのに！
タグ：
posted at 12:29:25
黒木さんは情報量基準だけで見ろと言ってるわけではなくて情報量基準も使ってもっとシビアにモデルの妥当性を検証しろって言ってるんじゃないかな。モデル自体は当てはまりの良さそうなのを過去の知見と検証目的に従ってとりあえずは作るだろうけど。
タグ：
posted at 12:17:58
「いま数学で必要なのは三角関数でも微積分でもなく、統計です。」って、依存関係とかガン無視してね？
【文系・理系の壁】駿台・石原賢一さんに聞く 文系のサクセスロードが見えない理由とは(朝日新聞EduA)
#Yahooニュース
news.yahoo.co.jp/articles/46f2c...
タグ： Yahooニュース
posted at 10:25:42
news.yahoo.co.jp/articles/46f2c...
この石原賢一氏は何を言ってるんだろうか？
逆ですよね。三角関数や微積分を学べば統計や論理だけでなく数学を使うあらゆる分野へ進む足がかりが得られるわけで。 pic.twitter.com/k2iQ72mIhe
タグ：
posted at 09:40:53
#統計 以下の2つの結果は、データとモデルの尤度函数への要約で失われる情報を事前分布で補う方法を与えています。
nbviewer.jupyter.org/gist/genkuroki...
止め方ごとに最大期待予測誤差を最小にする事前分布を決める計算
nbviewer.jupyter.org/gist/genkuroki...
止め方ごとに違うP値をベイズ統計の方法でぴったり実現する方法
タグ： 統計
posted at 09:36:36
#統計 統計学の哲学者のMayoさんのウェブサイトの
errorstatistics.com/2017/04/01/er-...
には、p-hackingに類似の
Bayes hacking
の節があります。
ベイズ統計の事後分布で測った確率は通常のP値と漸近的に一致するものが結構あり(場合によってはぴったり一致)、p-hackingと同じことをできます。 pic.twitter.com/f9AixbkM7k
タグ： 統計
posted at 09:14:20
#統計 この「尤度原理」をまともなもの扱いしている人達の問題はかなり深刻で、このスレの上の方でも紹介したように、「尤度原理」を使ってp-hackingと同じようなことをやってよいことにする人まで出て来てしまいました。(豊田『瀕死本』の尤度原理の節) twitter.com/genkuroki/stat...
タグ： 統計
posted at 08:40:14
#統計 赤池さんは論文で推定の場合のみを扱っていますが、仮説検定(ただし片側)の場合に「データの取得法が違うという事前情報を使って異なる事前分布を採用してよい」という考え方を適用したのが、以下のリンク先の話です。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 08:34:53
#統計 Beta(a,a)型の事前分布で最大期待予測誤差を最小化するものを数値計算で求めると、試行回数nを固定した場合にはa=0.5のJeffreys事前分布が得られ、成功回数x(私に記号ではk)を固定した場合でx=kが小さい場合にはa<0.5の別の事前分布が得られます(x=kを大きくするとa=0.5に近付く)。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 08:29:06
#統計 1つ前のツイートに引用した赤池さんの論文の画像の3枚目の上半分
赤池さんの論文ではこの「止め方の違い」の件にこれ以上詳しく説明していないのですが、その続きにやるべき計算をしたのが、既出の
nbviewer.jupyter.org/gist/genkuroki...
止め方ごとに最大期待予測誤差を最小にする事前分布を決める計算 pic.twitter.com/vx9njR07Bk
タグ： 統計
posted at 08:24:22
#統計 添付画像は
www.uv.es/~bernardo/vale...
↓
www.uv.es/bernardo/Valen... (←歴史的資料)
↓
Likelihood and the Bayes procedure
Hirotugu Akaike
1980
より。
* Savage批判
* 事前分布はデータの取得法と無関係に決めるべきだという考え方の否定
* ベルヌイ分布モデルの例(赤字コメントを参照！) twitter.com/genkuroki/stat... pic.twitter.com/h6R3GjRlC0
タグ： 統計
posted at 08:19:48
#数楽 厳密な証明にはもちろんさらなる議論が必要です。
厳密な証明は自力でやると楽しいです。
そのためのヒントになる数値計算の結果が以下のリンク先にあります。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 07:41:27
【12/25 7時更新】1週の累計感染者数が2人/人口10万人 以上の都道府県は40
東京 神奈川 広島 大阪 京都 高知 埼玉 愛知 福岡 兵庫 千葉 岡山 北海道 奈良 沖縄 岐阜 群馬 長崎 宮城 熊本 栃木 福島 滋賀 香川 石川 山梨 静岡 茨城 宮崎 山形 三重 富山 佐賀 鹿児島 山口 大分 長野 新潟 島根 和歌山 pic.twitter.com/Gyu2xhV1gA
タグ：
posted at 06:48:18
New post: Sets, chains and rules - part II - www.juliabloggers.com/sets-chains-an... #julialang pic.twitter.com/qX15Pe1MUn
タグ： julialang
posted at 02:52:07
Sets, chains and rules – part I www.juliabloggers.com/sets-chains-an... #juliabloggers
タグ： juliabloggers
posted at 02:28:59
非公開
タグ：
posted at xx:xx:xx
shift_jisエンコードのウェブページにHTTP.jlでログインしてソースを取得できた
これまでPHPで日々取得していたSBI証券の口座情報をJUliaに移植してSQLite3にエクスポート pic.twitter.com/GjZEOaAigP
タグ：
posted at 01:13:31
非公開
タグ：
posted at xx:xx:xx
#統計 データn=12, k=3が得られたとき、試行回数nを固定した場合と成功回数kを固定した場合では、仮説「成功確率は0.5以上」の片側検定のP値はそれぞれ7.3%と3.3%と異なります。
oku.edu.mie-u.ac.jp/~okumura/stat/...
「尤度原理」の信者は「これは仮説検定の欠陥である」と言うのです！(わら)😱😱😱
タグ： 統計
posted at 01:02:29
#統計 試行回数nを固定した場合と成功回数kを固定した場合では、最大期待予測誤差(リスクの指標)を最小にするBeta(a, a)型の事前分布が(特にkが小さいときに)異なることも示せます。
データの取得法によってリスクを最小化する事前分布は変えるべきなのです。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:57:39
#統計 例えば、試行回数nを固定した場合と成功回数kを固定した場合において、事前分布をそれぞれBeta(1,0), Beta(0,0)にすれば、事後分布で計算した仮説「成功確率はp以上である」の確率は対応する場合の片側検定のP値にぴったり一致します。続く twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:57:38
非公開
タグ：
posted at xx:xx:xx
次やるとしたら、KONTAKT KSP、UVI Script、u-he Hive、Logic Scripter、HISEみたいな学術研究の対象になりにくいベンダー言語とか、OpenMusicのような楽譜・シーケンス言語、あるいはPython、Rust、Go、R、Juliaなど汎用的な言語のオーディオ機能あたりかも。やるかわかんないけど。
（おわり）
タグ：
posted at 00:24:52
#統計 赤池弘次さんが事前分布を使うといいことができることを嬉しそうに語っているシーンの紹介↓
私はこういう話が非常に楽しいと感じます。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:24:30
ごまふあざらし(GomahuAzaras @MathSorcerer
Julia 1.6の個人的に気になる新機能 qiita.com/ho-oto/items/e... #Qiita
タグ： Qiita
posted at 00:17:31
#統計 続き。「複数の仮説の真偽に関する事前の主観的信念の度合いがデータに基く信念の更新によって変わって行く」という見方をして良い場合もあるでしょうが、Gelman-Shalizi (2013)では実践的なベイズ統計の使い方とは全然違うことを指摘されています。続く
www.stat.columbia.edu/~gelman/resear...
の第1節を参照 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:15:53
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
