黒木玄 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年02月20日(木)
Parallelizing a Specialized Matrix Type in #julialang
mdavezac.github.io/blog/
タグ: julialang
posted at 00:01:19
「小6の息子が消しゴム買ってっていうから、黒いのはどうしたの?って言ったら「先生がダメだって」と言うので、そんな訳ない普通の消しゴムだしと思って、明日もう一回聞いてみな、って理由を聞かせたら、本当にダメだったの!理由は「白じゃないから」なんだって!」という話を聞きクラクラした🌀🌀
タグ:
posted at 00:09:49
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
#統計 【この辺の注意を明瞭に書いてある解説は珍しいです】
その珍しい例
島谷健一郎著
『ポアソン分布・ポアソン回帰・ポアソン過程』
統計スポットライト・シリーズ2
近代科学社 2017
www.amazon.co.jp/dp/B076LVLLSF
Kullback-Leibler情報量が多項分布の確率の対数の漸近挙動に出て来ることの説明。 twitter.com/genkuroki/stat... pic.twitter.com/N5gZ1K41Vg
タグ: 統計
posted at 00:19:49
#統計 添付画像は p.85より。pp.83-84に詳しい計算が書いてある。
genkuroki.github.io/documents/2016...
の第1節
www.jstage.jst.go.jp/article/butsur...
p.610の右側
と同じ計算。 pic.twitter.com/y1JiW3wp34
タグ: 統計
posted at 00:23:05
非公開
タグ:
posted at xx:xx:xx
#統計 私のノート
nbviewer.jupyter.org/gist/genkuroki...
では、予測分布で比較しているのでベイズと最尤法が一致する(他のモデルでは普通は不可能な)事前分布は improper な "Beta(0,0)" になっていますが、『社会科学のためのベイズ統計モデリング』ではMAP法と最尤法が一致する一様事前分布を使っていますね。
タグ: 統計
posted at 00:50:21
#統計 ベルヌーイ分布モデルの場合もそうなのですが、簡単な確率モデルの式にはガンマ函数(やベータ函数)がよく出て来る。それらの導函数を扱うと、ディガンマ函数やトリガンマ函数が出て来る。 pic.twitter.com/OThDAj3Icv
タグ: 統計
posted at 01:19:34
@thomas_berg1 Nå har du grunn til å installere #JuliaLang Thomas! Med Cxx.jl pakken får du C++ på kommanlinjen om du vil 😉 Grei til å teste ut STL thing f.eks. uten å kompilere et helt program. github.com/JuliaInterop/C...
タグ: JuliaLang
posted at 01:38:38
Grassmann.jl uses the OG (Original Grassmann) definition for complement duality, compatible with differential geometry and Hodge-DeRahm co/homology theories ~ grassmann.crucialflow.com pic.twitter.com/9nvkpi6wzT
タグ:
posted at 01:45:07
Plotly's interactive parallel coordinates plot is really nice to build some insights on the impact of hyper-parameters on the generalization performance. pic.twitter.com/ekOVDqIjwi
タグ:
posted at 03:23:07
Here is the full notebook with the interactive plots at the bottom if you want to play with them:
nbviewer.jupyter.org/github/ogrisel...
Gradient Boosting to predict the price of houses in Ames. pic.twitter.com/eKjuay2aJ3
タグ:
posted at 03:23:10
By the way in that same notebook you can also see SHAP (master branch) explaining the decision of sklearn.ensemble.HistGradientBoostingRegressor: pic.twitter.com/tTAl6QKsUa
タグ:
posted at 03:26:44
Masahiro Ono 小野 昌弘 @masahirono
ただの感想ですが、岩田先生のビデオは国際的に通用する論理で話をしていると思う。高山先生の記事は、国内向けの論理ばかりで疑問に思ったーこれはたとえば、英語に訳したら意味不明な点が明確になるかなと思う。あと、後者は前者の信用を落とすための余計な文章が多くて、ちょっと残念に思った。
タグ:
posted at 04:52:53
仰るとおりです。問題の本質に反論できず、「現場はがんばってるから仕方ない」としか言ってません。それでも感情に訴えて流れをつかみ、勝てて(?)しまっています。これこそが日本を停滞させている本質だと思うのですが。 twitter.com/masahirono/sta...
タグ:
posted at 06:46:53
高山氏の投稿を読んで、素人の私が最も感じたのは、この非常事態においても横の機関に気を遣いながら仕事する余裕があるならまだまだ大丈夫なのかなという漠然とした安心感を与えてるよなって。(褒めてない)
タグ:
posted at 07:35:53
非公開
タグ:
posted at xx:xx:xx
『数とは何か』を教師志望の人に勧める場合には、その本が有名なトンデモ本であるシュペングラー『西洋の没落』を繰り返し引用しまくっていて、文化相対主義に異様に肩入れしていることへの警告が必要。
その本のそういう所を受け入れた人が中高生に数学を教えるようになるのは悪夢。 twitter.com/kamo_hiroyasu/...
タグ:
posted at 08:07:48
この機会に高山義浩先生を知った人も多いのかもしれませんけど、高山先生こそ世界と日本の地域医療の現場を知り抜いた第一線の医師ですよ。しかも若月俊一の農村医学と地域医療の系譜の最良の部分を受け継いで、現代医療に接続されてる方だと思います。この本とか読むといい。
www.amazon.co.jp/dp/4260028197/...
タグ:
posted at 08:37:03
個人の感想ですが、10年くらい見てきた岩田先生は、リアルは知らんけど、ネット上では、感染症の分野ではみんな認めるが、ざっくり言うと医療クラスタから感染症以外の信頼はあんまされてないってかんじ。
タグ:
posted at 08:52:50
なんぼでもネットでの争いネタはあると思うけど、女医も育休中にヒマを探して論文書いたりできるっしょという論争でやり取りするのはほぼやめました
感染症の分野では尊敬するべきかなと思うので本も読むが
タグ:
posted at 08:57:17
岩永直子 Naoko Iwanaga『言 @nonbeepanda
岩田先生が動画を削除したことを追記しました。なお、動画の内容に一部反論された高山義浩先生にも取材を申し込みましたが、ご自身ではないほうがいいとして、別の先生をご紹介いただいているところです。
www.buzzfeed.com/jp/naokoiwanag... @nonbeepanda
タグ:
posted at 09:24:04
#統計 以下のリンク先のWBICの計算でも、ガンマ函数(分配函数Zの一種)の対数の導函数(digamma函数)が出て来ている。
特殊函数というくくりで見ると、全然整理されていないように見える汚い計算が、「統計力学でよくやる計算だ」と知っていれば、見通しを立て易くなる。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 09:44:14
岩田健太郎 K Iwata, MD, M @georgebest1969
高山先生がFBから情報公開しています。彼と公の場で議論するのはよくないし目標でもないので細かいことへの反論はしません。ただ大きな誤解あるといけないので一点だけ。
タグ:
posted at 10:19:49
岩田健太郎 K Iwata, MD, M @georgebest1969
高山先生は「実際はゾーニングはしっかり行われています。完全ではないにせよ」とお書きになっていますができてなかったのは事実でよって感染の危険を強く感じました。派遣前「クルーズ船の中の本部を外に出すようぜひ進言してほしい。私も何度も主張しているのですが」とおっしゃったのが高山先生です
タグ:
posted at 10:19:50
#統計 統計力学でよくやる計算とは添付画像のような計算のこと。分配函数 Z(β) の対数導函数で、ハミルトニアン H(w) の平均と分散が自然に出て来る。 twitter.com/genkuroki/stat... pic.twitter.com/sR2YQ4UY8S
タグ: 統計
posted at 10:25:58
#統計 続き。以下のリンク先の添付画像の汚い手書きのメモではまさにそういう計算を行なっている。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 10:33:50
今回の騒動でひとつ良かったのは、高山義浩先生を知ったこと。
彼の書いた、インフルエンザについてみんなが知っておいたほうがよい基礎知識。
やっぱり丁寧だし、読みやすくて素晴らしいな。
インフルエンザについて診察室で説明していること www.huffingtonpost.jp/yoshihiro-taka... @HuffPostJapanさんから
タグ:
posted at 10:36:24
このニュース,大学入試の共通試験に対する信頼性,credibilityを根本から揺るがす大事件なのだけれど,英語の共通試験に関してはベネッセがこれをしてもいいと文科省はお墨付きを出してましたもんね.
東大や京大の試験作成の責任者が対策問題集を出すようなもん.この視学官も,そのうち天下るのか twitter.com/kotani_eisuke/...
タグ:
posted at 11:08:30
シュペングラーの引用が非本質的な一部分に限るのであれば、大して問題にならないと判断するのですが(完璧な本など存在しない)、全体が影響を受けているように見えた。
私ならばその本を他人には絶対にすすめない。
タグ:
posted at 13:45:28
あと、大森荘蔵の時間三部作やそこに書いてあることを褒め称えているように見える本も他人には絶対にすすめない。
「権威」ある思想的なトンデモ系は本当に有害だと思う。
タグ:
posted at 13:48:28
【募集】
コロナでの経済的損失、大打撃すぎてただただ傍観してるだけじゃ辛すぎる……
「お客さんが激減してる!」という飲食店やホテル、旅館、サービス関係のみなさま、よければリプ欄でぜひ宣伝してもらえませんか…
イチオシの写真と、Googleマップ等住所情報もつけていただければ最高です…
タグ:
posted at 13:56:21
#統計 しかし、ベイズ自由エネルギーの近似計算としてのWBICについては原論文
www.jmlr.org/papers/volume1...
がわかりやすいと思います。
これを読めば、逆温度を巧妙に使える理屈もわかるし、コンピュータでの計算もできると思う。
タグ: 統計
posted at 13:57:25
黒木さんのSanovの定理の説明は、学部2年次くらいまでの数学で書かれていて、扱う内容に比してその意味で「易しい」というのが私の数学の先生の話。渡辺ベイズ本は、ルベーグ積分など、ガチの数学科の人がやる内容も載っているから、それは流して読むしかないんじゃないの?とのこと。
タグ:
posted at 13:59:48
じゃあまあ、自分は順調に数学の勉強を続けられたとして、2年後くらいにSanovの定理がわかれば嬉しいかな…。『社会科学のための~』をポチッたので、それも読んでみるつもりだが、どっちにしろ数学を進めないと手が出なそう。
タグ:
posted at 14:02:37
#統計 大学1年生向けの微積分の題材として、ガンマ函数函数は定番だと思うのですが、
「ガンマ函数の対数導函数達(高階を含む)にはディガンマ、トリガンマ、ポリガンマなどの名前が付いていて、コンピュータで使える特殊函数ライブラリで大抵の場合に使える」
という注意もした方が良さそうですね。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:05:16
#統計 訂正
❌ハミルトニアンという名の確率変数H(X)という名のモーメント母函数
⭕️ ハミルトニアンという名の確率変数H(X)のモーメント母函数
この手のミスが結構多い。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:36:02
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
書きました。はてなブログに投稿しました #はてなブログ
3.2.2:カテゴリ分布の学習と予測【緑ベイズ入門のノート】 - からっぽのしょこ
www.anarchive-beta.com/entry/2020/02/...
タグ: はてなブログ
posted at 18:00:25
#数楽 【黒木さんのSanovの定理の説明は、学部2年次くらいまでの数学で書かれていて、扱う内容に比してその意味で「易しい」】
こういうコメントはめちゃくちゃうれしい。
genkuroki.github.io/documents/2016...
Kullback-Leibler 情報量と Sanov の定理
の想定読者は、大学2〜3年生なのでドンピシャ! twitter.com/yosuke_hash/st...
タグ: 数楽
posted at 18:50:56
#数楽 予備知識が少なくてもすむように、Stirlingの公式の導出の仕方の一例の説明から始まる。
高校数学IIIの教科書を読んでいれば、不定積分∫e^{-x} x^n dxに繰り返し部分積分を適用してΓ(n+1)=n!を自分で証明できるはず。
結構、この「自分でできる」が数学では大事。
genkuroki.github.io/documents/2016... pic.twitter.com/YJggSCb6CE
タグ: 数楽
posted at 19:05:01
私の様なCPLEXとかGurobiソルバーが使えない「一般市民」は代替えソルバーはどんなものがあるか真剣に悩んでいます。#Julia言語
JuMP
www.juliaopt.org/JuMP.jl/dev/in...
タグ: Julia言語
posted at 19:06:56
float f,g;void setup(){size(720,720);}int F(float x,float y,float p,float q,int d){float s=noise(g++,f);fill(s*255);rect(x,y,p,q);return d++>8?1:d%2<1?F(x,y,s*p,q,d)+F(x+p,y,s*p-p,q,d):F(x,y,p,s*q,d)+F(x,y+q,p,s*q-q,d);}void draw(){F(g=0,0,720,720,1);f+=.01;}#つぶやきProcessing pic.twitter.com/GlU27t7WaL
タグ:
posted at 19:15:03
不謹慎だと怒られそうだけど、岩田医師と高山医師の組合せキャラのバランスが絶妙で、ドグラ・マグラにおける正木博士と若林博士みたいだなて思った。こういうの腐ってる界隈が好物なんじゃないかと思うしもう薄い本書きはじめてる人がいてもおかしくなくてそのばあい多分高山受
タグ:
posted at 20:54:37
#数楽 渡辺澄夫『ベイズ統計の理論と方法』の第4章の内容は、特異点解消定理さえ認めてしまえば、
genkuroki.github.io/documents/2016...
一般化されたLaplaceの方法
(ガンマ函数を使って頑張る計算)
または
nbviewer.jupyter.org/github/genkuro...
分配函数のゼータ函数
の方法でかなり初等化できます。後者の方がシンプル。
タグ: 数楽
posted at 21:11:07
【私が関係者に聞いた、岩田先生がダイヤモンドプリンセス号から下された事件の真相は、こちらの内容に近かったです】のように「私が関係者に聞いた~の真相」という言い方は非常に良くないと思う。このタイプの発言に頼ろうとした人達は恥を知るべき。
その「関係者」の名前を明示するべき。 twitter.com/yusuke_tsugawa...
タグ:
posted at 21:17:32
「知り合いから聞いた」タイプの情報はデマである可能性が高いとみなされるべき。
そういうタイプの言い方で社会的地位によって信用されそうな人が情報を流すことが許されると、我々の社会はデマに対してより脆弱になってしまうと思う。
タグ:
posted at 21:20:39
以下のリンク先の添付画像を見たら、岩田健太郎さんに対して「国民を不安と混乱に陥れる発言はどうぞ控えて下さい」とありますが、私のTLには不安と混乱に陥った国民は見当たらず、
* 岩田健太郎氏はすごい専門家である。
* 岩田健太郎氏は大迷惑な奴だ。
という話が流れまくっていて笑っている。 twitter.com/j_geiste/statu...
タグ:
posted at 21:35:11
#統計 事前分布φ(w)の分配函数の対数の-1倍
Ψ(β) = -log Z(β) = -log∫exp(-βH(w))φ(w)dw
とカノニカル分布 p(w|β)= exp(-βH(w))φ(w)/Z(β) に関する平均をE[ |β]と書くと、
Ψ(0)=0, Ψ'(β)=E[H(w)|β]
なので
Ψ(β)=∫_0^β E[H(w)|β]dβ.
この表示がWBICの出発点。続く
twitter.com/genkuroki/stat...
タグ: 統計
posted at 21:50:05
Dr. Chris Rackauckas @ChrisRackauckas
This is absolutely amazing, and indeed it's design auto composes with DifferentialEquations.jl #julialang. BlockBandedMatricrs already gives around a 5x speedup over standard sparse techniques for implicit PDE solving: this would only compound that advantage. twitter.com/ArturoErdely/s...
タグ: julialang
posted at 21:56:59
Dr. Chris Rackauckas @ChrisRackauckas
I plan to write some semilinear PDE tutorials using all of the new tooling soon.
タグ:
posted at 21:57:32
岩田さんの主張やなくて人柄ややり方を批判する方々がすごい勢いで湧いてきてるのを眺めていると、ああ日本は実に陰口といじめの国なのだなあとしみじみ感じる。正論が目立ったり秩序を乱すならば「あいつはダメだ」「やり方がダメだ」と横槍で攻撃が始まる。オレはこういう国じゃないようにしたい
タグ:
posted at 21:59:48
twitter.com/0211jaguar/sta... まさにこういうのがそれですね。いじめの国の物言いの典型だと思います twitter.com/0211jaguar/sta...
タグ:
posted at 22:06:10
twitter.com/kyoujikyouji/s... これもいじめの論理ですね。関係ないことを持ち出して人格を非難することに夢中になるのがいじめの国です。そりゃグラミー賞の人も出ていくよね… twitter.com/kyoujikyouji/s...
タグ:
posted at 22:08:46
「クリスマスシーズンの12月23日が天皇誕生日でなかったことや、台風・暖冬などで12月の消費の落ち込みが想定以上だった」
台風・暖冬は関係ないというデータについては→note.com/iida_yasuyuki/...
祝日云々は,季節調整データは曜日の影響を調整しているので筋違い
jp.reuters.com/article/japan-...
タグ:
posted at 22:21:15
読み始めた.まだ序盤だがついていける.計算過程については極めて丁寧なので,高校数学で十分追える(慣れている人からすると冗長かも).しかし一方で,不親切とまでは言えないが,初心者に対してはもう少し説明が必要と思われるところもある.
タグ:
posted at 22:38:01
#統計 『社会科学のためのベイズ統計モデリング』の読者は
須山敦志『ベイズ推論による機械学習』
www.amazon.co.jp/dp/B07L2V4H59
のセンスの良い計算の仕方も参考になるかもしれません。統計モデリングでは各種の確率分布をレゴブロックのように組み合わせて使うので、簡単な例の様子を見ることは大事。
タグ: 統計
posted at 22:41:13
#統計
github.com/sammy-suyama/B...
須山敦志『ベイズ推論による機械学習入門』のサンプルコード
正誤表もここにあります。
この本の第2章「基本的な確率分布」はレゴでの各種ブロックの解説だと思ってよい。『社会科学のためのベイズ統計モデリング』の第5章より詳しく、カラーの図も豊富です。
タグ: 統計
posted at 22:47:41
「確率モデル」の初出時(2章 p.25)の説明で,各サンプルがi.i.dの下でベルヌーイ分布に従うこと(Y_i ~Bernoulli(q))を例に挙げているが,「ベルヌーイ分布」を確率モデルと呼ぶと書いた直後に,確率モデルは実現値はこの確率分布から生成されたという「仮定」を表すとあったり・・・
タグ:
posted at 22:51:15
つまり,この例の場合に確率モデルが次のうちのどちらを指すのかちょっと混乱した
(a) 確率モデル = Bernoulli(q)
(b) 確率モデル = Y_i ~Bernoulli(q))(i.i.d)
モデルというからには,(b)という気がするが
タグ:
posted at 22:54:14
#統計 あびこ(牛)さんも紹介なさっていますが、
watanabe-www.math.dis.titech.ac.jp/users/fujiwara...
特異モデルにおけるベイズ検定と変化点発見への応用
藤原香織 渡辺澄夫
の話は「仮説検定とベイズ統計は水と油だ」のような印象論が__誤り__であることの証拠になっている。
間違った印象論を広めるような教育はまずい twitter.com/abiko_ushi/sta...
タグ: 統計
posted at 23:31:01
PseuDoctor@フォロワーさん以外 @_pseudoctor
続き)そもそも今回のCOVID-19に関してもメディアに出て来るのは医者ばっかり。専門家と一般市民の間の橋渡しをしてくれる筈の人達がどっかに居た様な気がするんだけど…
まぁ震災直後みたいにアカラサマな妨害を仕掛けてこないだけマシかもしれないけどね(某医療ジャーナリスト氏()を除く)。
タグ:
posted at 23:51:24