黒木玄 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 を見て下さい。
2021年02月07日(日)
「Julia言語で入門するプログラミング」第7回を書きました。多重ディスパッチ、Union型、パラメトリック型、Holy Traitsパターンといったあたりの話題に触れています。難易度の高い文法事項で読むのが大変なので、ボリュームは少なくしています。 #Julia言語
muuuminsan.hatenablog.com/entry/2021/02/...
タグ: Julia言語
posted at 03:05:04
私のブログ、サイトに移行したくなってきてる。とりあえずスモールスタートだと無料のはてなブログからスタートしたけど。
やるならWordpressかなー、と思ってたけど、どうせなら JuliaのWebフレームワークを使ってみたい気がする。でも構築にそっちにハマって移行が進まない懸念が。
タグ:
posted at 08:09:29
Akira Suzuki | Engin @akr_suzuki
そのPythonもJuliaへの流れが既に起きています。 / 2021年 ITエンジニアが学びたいプログラミング言語、第1位はPython (TECH+) #NewsPicks news.mynavi.jp/article/202102...
タグ: NewsPicks
posted at 09:51:47
今日の発表の補助的資料ですが, 最適フィードバック制御モデル (optimal feedback control; OFC)をTodorov先生のMATLAB実装を参考にJulia実装した記事を作成しました.
compneuro-julia.github.io/chap13/13-3_op... pic.twitter.com/mA38gbof4t
タグ:
posted at 10:18:05
非公開
タグ:
posted at xx:xx:xx
仕事で静的型付けの言語をたくさん書いた後にJuliaに戻るとそのシンプルさと自由さに感動して改めて自分の好きな言語なんだなって確認できる一方、静的型付け言語(今回はScala)を書いてる時に使える素晴らしいIDE機能達を早くJuliaにも持ち込みたいなって気持ちになりますね(metals-vscode良かった)
タグ:
posted at 10:46:49
It's been pretty nice experience to write Scala with metals-vscode. I like it :)
scalameta.org/metals/docs/ed...
タグ:
posted at 10:49:17
母は手工芸が得意で、晩年は木彫りを楽しんでいました。せめてもの供養にと、母の作品の一部を皆さんにも見て頂こうと思います。(この木彫り作品の柄は母のオリジナルデザインです)こういう作品を多数作っていました。 pic.twitter.com/c3wc8jaZ4B
タグ:
posted at 11:14:26
#統計
以前、ガンマ分布の最尤推定を考えたことがあり、添付の方法を自力で思いついた(先人が発表している可能性は大いにあり)。簡易ニュートン法で更新式が簡単。
速さは結構いいと思う。 twitter.com/genkuroki/stat... pic.twitter.com/NG39pfotTR
タグ: 統計
posted at 12:00:29
ごまふあざらし(GomahuAzaras @MathSorcerer
#Julia言語
#JuliaLang
JSServe.jl を使って WGLMakie のグラフをブラウザに表示することができました. pic.twitter.com/fFy3g5gsso
posted at 14:10:48
ごまふあざらし(GomahuAzaras @MathSorcerer
jl.dev という #Julia言語 開発オンラインでやった成果です.
github.com/terasakisatosh...
タグ: Julia言語
posted at 14:10:48
#Julia言語
nbviewer.jupyter.org/gist/genkuroki...
@ timeと@ code_warntypeの使い方
を大幅に拡充。
@ time や @ code_warntype に余力函数の検査を行わずに、使い回しをするコードを大量に書いた後で、大きな速度的な劣化を招く書き方をしていることに気付くと、かなりの時間を無駄にすることになります。
タグ: Julia言語
posted at 14:37:14
#Julia言語 あと「Juliaでは静的型付けもできるので速い」というような「静的型付け必須カルト」的な誤解をしてしまうことも、@ code_warntype を見ていれば防げます。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/pxBfyNVfAW
タグ: Julia言語
posted at 14:46:50
最近模索中の統計モデリングの使い方は、汎用フィルタあるいは特徴量ジェネレータとしての使い方。例えば、時系列を複雑な要素分解して、ローパスフィルタでは得られないような推定した値を別の機械学習の入力に使う。階層モデリングして色々な〇〇差を除去した個人ごとの値を入力に使う、など。
タグ:
posted at 15:01:42
ほんとは全部丸ごと統計モデリングにぶち込めるのが理想。しかし、尤度の綱引きが解析者に意図通りにいかないこともあったり、速度的に限界があったり、出力部分で性能の良いブラックボックスなモデルを使う方が良い場面もあったりするので。
タグ:
posted at 15:05:46
坂口安吾の随筆読んでたら「友達がメチルで死んでから体に気を遣って飲むようになった、仕事中はヒロポンで5日くらい徹夜するから泥酔せんと眠れないので書き終えたらウイスキー2本飲んで気絶してる」とか書いててマジ文壇のパワーを見せつけられた、ストゼロで退廃気取ってるネット人類とは格が違う
タグ:
posted at 15:06:38
#数楽 個人的な意見では、三角函数の加法定理は中学校数学の幾何の範囲内で理解できる非常に面白い話なので、特に時間をかけて教えることは極めて教育的だと思います。(三角函数の加法定理にはラジアンの知識は不要なので中学校レベル)
直角三角形の並べ方を変えても本質的に同じ公式が得られる。 twitter.com/genkuroki/stat... pic.twitter.com/5iQlScpD9g
タグ: 数楽
posted at 15:20:39
My #julialang package JET.jl is gradually growing and now its source code has >3000 lines of Julia code. I'm still happy that Julia's powerful AST macro allows me to abstract away the duplicated works and helps me maintain such a certain amount of code :)
タグ: julialang
posted at 15:39:18
Currently within JET.jl's src directory, I have 10 macro definitions. You can have a look at some of them here:
github.com/aviatesk/JET.j...
For example, the `at-withmixedhash` macro I added today effectively abstracts the following kind of pattern:
github.com/aviatesk/JET.j... pic.twitter.com/BFYQw5U5pF
タグ:
posted at 15:39:19
I found some of those macros might be useful generally, and I actually upstreamed into Julia Base; `at-invokelatest` and `at-invoke` macros will be available as of v1.6. I'm happy if you find it interesting.
github.com/JuliaLang/juli...
github.com/JuliaLang/juli...
タグ:
posted at 15:39:19
ネットの記事で見た渡辺先生のベイズ統計の本、Kindle 版がなかったから久々に紙の本で買ったけど、パラパラ覗いたらイメージしてたより読みやすそうだった。
タグ:
posted at 15:44:57
統計学の本に関しては、数学的にきっちり書いてある本の方が好みだな。。。変にゆるく説明してある本だと、細かいロジックが気になってかえって悩まされることが多い。😅
タグ:
posted at 15:46:53
前にブログのネタにしたことあるけど、平均値の線形性
E [a X + b Y ] = a E [X ] + b E [Y ]
とかも、ちゃんと証明するにはフビニの定理とか必要なはずなんだけど、雑にしか説明してない本が多い。
タグ:
posted at 15:54:10
RC1!!! 🎉 Julia v1.6.0-rc1 is now available - Announcements - JuliaLang discourse.julialang.org/t/julia-v1-6-0...
タグ:
posted at 16:07:28
E[aX+bY]=∫_Ω (aX(ω)+bY(ω))dμ(ω) がwell-definedな時点で積分の線形性から出る話なので、Fubiniの定理は無関係です。
定義域が違うXとYの定義域の直積に適切な測度を入れる話はFubiniの定理に関係しますが、期待値汎函数の線形性はXとYが独立でなくても成立しているのでやはり違う話題になります。 twitter.com/r1ms31dk/statu...
タグ:
posted at 16:10:11
#数楽 それはそうだとしても、渡辺澄雄さんが書いた本はどれも読み易いと私も思います。
個人的にはベイズ統計に関する解説では最も(考え方が)易しく解説してあると思う。
難しいこと__も__書いてあることは確かですが、易しく読めることも確かだと思う。 twitter.com/r1ms31dk/statu...
タグ: 数楽
posted at 16:12:53
久しぶりに3次元イジングモデルをだらだらといじって遊んでいます。モンテカルロシミュレーションの「適切な手順」をチュートリアルにしておきたいわけですが、まあとりあえず計算が動くだけで面白い。50^3なんて、僕が博士課程で研究してたときには大変だったのに、今やMacbookでできるもんな
タグ:
posted at 16:16:54
非公開
タグ:
posted at xx:xx:xx
#数楽 渡辺澄雄『ベイズ統計の理論と方法』の第4章での分配函数のゼータ函数の使い方については、私の解説
nbviewer.jupyter.org/github/genkuro...
が現存する最も易しい解説だと思います。
分配函数の漸近挙動については大学1年レベルの直接計算も可能です。次のリンク先を参照。
genkuroki.github.io/documents/2016...
タグ: 数楽
posted at 16:24:08
適切な手順としては、まず初期緩和に必要な時間を見積もっておく必要があり、これにはそこそこの試行錯誤が必要。非平衡緩和の緩和時間と平衡での緩和時間は同じなんだけど、念を入れるなら平衡緩和時間も出しておく。これで必要なステップ数がだいたい分かる。あとは独立なランを10回くらいやる
タグ:
posted at 16:26:14
富谷(助教);監修 シン仮面ライダー @TomiyaAkio
当たり判定【プログラミング講座】 youtu.be/DXajm0sHfp0 via @YouTube
タグ:
posted at 16:26:19
#数楽 特異モデルではなく、指数型分布族の特に易しい正則モデルの場合のWAICの公式を作っておいても理解が捗ります。
Bernoulli分布モデルの場合↓
nbviewer.jupyter.org/gist/genkuroki...
正規分布モデルの場合↓
nbviewer.jupyter.org/gist/genkuroki...
自由エネルギー→LOOCV→WBIC→WAICの順に公式の導出が面倒になる。
タグ: 数楽
posted at 16:30:50
#数楽 最もシンプルな特異モデルである混合正規分布モデルの場合の尤度函数の形については
nbviewer.jupyter.org/github/genkuro...
尤度函数のプロット
を参照。
平坦事前分布の事後分布∝尤度函数
なので尤度函数の形を沢山知っておかないと、ベイズ統計については理解不可能になります。 twitter.com/genkuroki/stat...
タグ: 数楽
posted at 16:37:35
#統計 最も易しい特異モデルの混合正規分布モデルと正規分布モデルのWAICによる比較については
nbviewer.jupyter.org/gist/genkuroki...
1次元の混合正規分布モデルのWAIC twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:43:02
非公開
タグ:
posted at xx:xx:xx
#統計 そういう色の付け方をしたプロットの例↓
汎化誤差=予測誤差+定数です。
WAICの側は原点を各サンプルごとに適切にずらしてあります。
適切なずらし方は渡辺澄雄『ベイズ統計の理論と方法』のp.119, 定理15をみれば分かる。
こういうプロットが私以外の人から沢山出て来るとうれしいです。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 17:37:40
Julia nativeの機械学習ライブラリあるんだ。そのうちjuliaからPythonに逆輸入とかあるだろうか twitter.com/poyothon/statu...
タグ:
posted at 17:38:26
八代市立太田郷小学校のサイト、どうも時々更新されているような気がする。
【※このかけ算順序問題は、大学の先生たちなどの間でもいまだ論争があるそうです】
この記述、以前は無かったと思う。
es.higo.ed.jp/otago/blogs/bl... #掛算
タグ: 掛算
posted at 17:40:55
#Julia言語 入れた。
一度、v1.6以上を使うとv1.5以前には戻れない感じ。
Julia v1.6.0-rc1 is now available - Announcements - JuliaLang discourse.julialang.org/t/julia-v1-6-0...
タグ: Julia言語
posted at 18:40:23
HAYASHI Tomohiro @SokoranoKumasan
うっわ。リテラシー低いな。
十年前で現実への認識止まってんじゃん。 twitter.com/hirox246/statu...
タグ:
posted at 18:45:23
#統計
rpubs.com/hoxo_m/255994
で紹介されているガンマ分布の最尤法の数値計算法は、素直にかつシンプルにみんな知っているニュートン法だとみなす方が分かりやすいです。
添付画像のノートを見てください。
零点の近くでグラフが直線でよく近似されているとニュートン法は爆速になります。 twitter.com/shoei05/status... pic.twitter.com/1IXSPOqxem
タグ: 統計
posted at 19:08:25
#数楽 添付画像1の状況でニュートン法を使うよりも、添付画像2の状況でニュートン法を使う方が圧倒的に収束が速いことは、ニュートン法を知っている人にとっては明らかだと思う。
添付画像1の状況で零点を近似計算したいなら、"x"をかけて添付画像2の状況に変換してから、ニュートン法を使うべき。 pic.twitter.com/Sk9lh2HMTW
タグ: 数楽
posted at 19:15:52
#数楽 微少修正版
f(α) = 0 に直接的にNewton法を適用せずに、αf(α) = 0 にNewton法を適用するだけの話。これはガンマ分布の最尤法に限らず広く使えるテクニック。
ガンマ分布の最尤法における工夫されたNewton法 pic.twitter.com/uSMQjHYZzk
タグ: 数楽
posted at 19:25:27
寄稿しました。お読みいただければ幸いです。
囲碁プロ棋士の賞金ランキング トップ10棋士の意外な「使い道」とは|NEWSポストセブン www.news-postseven.com/archives/20210... #newspostseven
タグ: newspostseven
posted at 19:29:27
#数楽 比較のための画像を作り直した。
添付画像1と添付画像2のどちらがニュートン法にとってありがたいか?(笑)
ψ(a) = (d/da)log Γ(a) はディガンマ函数と呼ばれる基本特殊函数で、コンピュータの基本特殊函数ライブラリで計算できます。導函数 ψ'(a) はトリガンマ函数と呼ばれている。 pic.twitter.com/R2rbmSOyv3
タグ: 数楽
posted at 19:35:00
#統計 去年の2月の終わり頃に出た新型コロナの潜伏期間の推定の論文
pubmed.ncbi.nlm.nih.gov/32046819/
の結果の再現を、
nbviewer.jupyter.org/gist/genkuroki...
でWAICの実装をするついでにやりました。
添付画像は原論文にはない予測分布のプロット。
ここでもフィッティングの対象にガンマ分布も使われています。 twitter.com/genkuroki/stat... pic.twitter.com/frghVA7j3S
タグ: 統計
posted at 20:14:05
Julia v1.6.0-rc1 is now available
discourse.julialang.org/t/julia-v1-6-0...
#juliaLang
タグ: juliaLang
posted at 20:49:33
非公開
タグ:
posted at xx:xx:xx
なんか延々とはてなブログに快適に投稿するための環境整備をしている
はてなブログに投稿しました #はてなブログ
Markdownの数式をKaTeXで高速で表示できるようにした - あおいろメモ
solkul.hatenablog.com/entry/2021/02/...
タグ: はてなブログ
posted at 21:17:56
Finally getting around to fixing my least favorite thing about JSON3.jl: github.com/quinnj/JSON3.j...
#julialang
タグ: julialang
posted at 21:31:47
超ざっくりとシャコパンチの4棒リンク機構モデル(McHenry+2012)を割り箸と輪ゴムで。ゴムで止めている位置が二種類あります。動画の前半と後半で位置が変わると入出力関係が変わります。てこと同じことだなということがこれで体感できます。後半で、より小さい運動入力で大きい出力がでます。 pic.twitter.com/shlnglWWOh
タグ:
posted at 21:52:04
たとえば、McHenry+2012 JEBの図1から:
jeb.biologists.org/content/215/7/...
a-b-c-d が輪ゴムの部分になっています。
シャコの種類が変わるとこの形が変わりますのでシャコごとに部位を計測するとこでどのような入出力関係になるかが推測できます。この論文ではさらに「拳」にかかる水の抵抗を推測しています。 pic.twitter.com/w3gYmjI8mq
タグ:
posted at 22:08:33
訂正:
❌@ time や @ code_warntype に余力函数の検査を行わずに
⭕️@ time や @ code_warntype による函数の検査を行わずに twitter.com/genkuroki/stat...
タグ:
posted at 22:19:38
OokuboTact 大久保中二病中年 @OokuboTact
#超算数
小学校での「0は偶数」の教え方。
丁寧に教えている。
意外だ! pic.twitter.com/jzAaJ4JmQv
タグ: 超算数
posted at 22:38:19
@k46347506 インストールせずにというか移行はできます。ユーザのホームディレクトリ以下に .julia/environments/v1.5/Project.toml というファイルがあるのでそれを .julia/environments/v1.6 にコピーして、Julia1.6初回起動後にPkgREPLでupすれば大体いけるはず
タグ:
posted at 22:45:28
OokuboTact 大久保中二病中年 @OokuboTact
#超算数
「倍数」に関して、「0」の扱いの説明が素っ気ない。 pic.twitter.com/0555oz4Ypo
タグ: 超算数
posted at 22:46:28
@KB_satou @genkuroki Unicode 3 switched the code point definitions and since then it caused confusion in fonts
タグ:
posted at 22:57:36
反ワクチンから「自分は接種したら証明を公開するのか?」と訊かれての反応。
反ワクチンの人は
「推進派は、自分は受けないで他人に勧めている」
という理解なのか。物が歪んで見えるのだな。 twitter.com/MIKITO_777/sta...
タグ:
posted at 23:15:52
"成年年齢引き下げで未成年者取消権が18歳まで引き下がる1年前にこの改悪。高齢者の消費者被害がこれでもかというくらいに蔓延しているときにこういう仕打ち。
身近な消費者被害を防いだり救ったりできなくなります"
現在消費者庁が進めている「特定商取引法契約書面デジタル化」、許されません。 twitter.com/syoyuri/status...
タグ:
posted at 23:19:55
非公開
タグ:
posted at xx:xx:xx
【静的型づけ言語みたく変数に型を書けばいいというわけではなくて,賢いjuliaのコンパイラが推論しやすいように書く】
まさにそれ。
函数の引数に型を書いても #Julia言語 だと全く速くなりません。
Juliaが函数の引数の型を他のモノの型に適切に伝搬させるように書くと速くなる。 twitter.com/sf_trainee/sta...
タグ: Julia言語
posted at 23:25:15
ジャパンライフ問題で天下りが指摘されたので、そう見られても当然ですよね・・・
ただ、マルチ商法や訪問販売等を規制する特商法の今回の改正に、「契約書デジタル化」が入ったのは自民党の消費者問題調査会からのようです。
www.nikkei.com/article/DGXZQO... twitter.com/hyourokudou_6/...
タグ:
posted at 23:32:43
ごまふあざらし(GomahuAzaras @MathSorcerer
AbstractAlgebraで計算させたいな。
#JuliaLang twitter.com/ilovegalois/st...
タグ: JuliaLang
posted at 23:38:48