黒木玄 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 を見て下さい。
2018年02月05日(月)
２刷り以降は直っておりますが、『我々はなぜ我々だけなのか』(川端裕人著・海部陽介監修)１刷りにある間違えのリストです。#我なぜ twitter.com/Rsider/status/...
タグ： 我なぜ
posted at 23:58:42
非公開
タグ：
posted at xx:xx:xx
須山敦志 Suyama Atsushi @sammy_suyama
僕は深層生成モデルは使わない（使いたいと思う機会がない）ですが，normalizing flowは一般的な潜在変数モデルにも使えるのでHMCの代わりに使ってみたいですね．
タグ：
posted at 23:42:44
「お金いれたのに出てこない。おかしい。壊れてるぞ。」って言いながらガチャをガシガシ蹴ってる感じの財務省主計局次長のコメントの右側に「水道から錆び混じりの水が出てくる…なかなか直してもらえない…」っていうどこかの准教授のコメントが映える twitter.com/Shusaku_SASAKI...
タグ：
posted at 23:40:22
須山敦志 Suyama Atsushi @sammy_suyama
ガウス分布からのサンプルに非線形関数を何回か潜らせれば，徐々に複雑な確率分布からのサンプルに変化していきます．これを深層生成モデルなどの事後分布の近似に使うのがnormalizing flowと呼ばれる手法です．
arxiv.org/abs/1505.05770 pic.twitter.com/poGkoxmY2v
タグ：
posted at 23:39:12
非公開
タグ：
posted at xx:xx:xx
twitter.com/hiding_koukyo/...
#Julia言語
nbviewer.jupyter.org/gist/genkuroki...
カルマン渦
にベクトル場と "p" の動画を追加しました。
MAXSTEP = 10000 にしました。倍に増やしたので、その分だけ計算時間は余計にかかっています。計算135秒、以下の動画26秒でした。 pic.twitter.com/iuznbjQccy
タグ： Julia言語
posted at 22:31:59
@ahe100 #Julia言語 でGIF動画を作るもっと小さな例が以下の場所にあります。
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/gist/genkuroki...
現在試しているのは最後のやつの一番下の方にあるPlots+GRを使う方法。
タグ： Julia言語
posted at 22:27:11
@ahe100 #Julia言語 を使って数学GIF動画を気軽に作っています。
ソースファイルは
nbviewer.jupyter.org/gist/genkuroki...
で公開しています。ただし現在ではそこに書いてあるよりも高速な方法をメインで使いつつあります。
juliabox.com ならブラウザ経由でJulia言語を使えます。Julia言語について要検索。
タグ： Julia言語
posted at 22:23:58
#Julia言語 やっていることは動ける範囲に制限を設けたランダムウォークに過ぎません。資産の加法平均には1上限を設けるが、対数平均(expすれば相乗平均)に下限を設けないと、対数平均は-γ=-0.5772に収束します。対数平均には「平均効用」の意味を持たせられる。
pic.twitter.com/r3n968MjN4
タグ： Julia言語
posted at 22:09:36
#Julia言語 資産の分布のランダムウォークの動画作成のソースファイルへのリンクをはるのを忘れていた。
nbviewer.jupyter.org/gist/genkuroki...
です。平均に上限を設けるだけではなく、対数平均に下限を設けると分布はガンマ分布に収束します。 pic.twitter.com/0Hvm7mAmsj
タグ： Julia言語
posted at 22:06:55
@genkuroki ただ、例えばDistributionのパッケージや、（以前解説してくださった）VectorとMatrixの違い、PyPlotの呼び出し(Symbolなど)、言語としてPythonとの対応関係を探していたのです。ならそう書けばよかったのですが、この辺り言葉足らずだったと反省してます
タグ：
posted at 21:28:51
@genkuroki 失礼しました。どうも自分のほうが勝手に非難されたと感じてしまったようです。数学のロジックをすっ飛ばして手抜きして理解しようとしている・・・勝手な勘違いだったと思います。
タグ：
posted at 21:25:18
ようやく試してみようとしたが…昼間の職場の窓10JuliaProはガンガンエラー吐いたし，家の林檎本ProでNotebookで試したらGR入れてと言われて入れたけど，大量のエラー吐かれて挫折。やはり自分で確かめながらでないと使えんのか。 twitter.com/tsatie/status/...
タグ：
posted at 20:44:20
@MathSorcerer 須山さんのウェブサイトを読んで「これは面白い」と思った人は買って損無し。計算も丁寧に書いてあります。
第4章の近似推論の章の解説は特に良いです。
第5章も楽しい。
その良さは見れば多分すぐにわかる。
娯楽にもなる良い本。
タグ：
posted at 20:03:30
@Dsuke_KATO 「須山さんの本のサンプルコードを理解しようと思ったら、須山さんが解説している数学の内容を理解しなければいけないだろうし、そうするのが近道だ」と私は主張しています。それに反対であればそのように述べればよいと思う。私が言っていないことについて私は知りません。私は誰も非難していません。
タグ：
posted at 19:45:45
@Dsuke_KATO 嫌がらせのつもりはまったくありません。
私はどうすれば理解できるかについて本当の話を書いたつもりです。
実際、私の発言は完全にそういう内容になっていると思います。
タグ：
posted at 19:39:32
#Julia言語 ではπやeだけではなくEuler定数γも大事に扱われている点がうれしい。
Euler定数 γ の-1倍は「一人当たりの平均資産を1以下に制限にして非負の資産分布をランダムウォークさせたときの、資産の対数平均の収束先」の形で自然に出て来る。
mean(log(X)) が-γ=-0.5772…に収束している。 pic.twitter.com/r3n968MjN4
タグ： Julia言語
posted at 19:32:56
「買うかも」が半分買うとすると70%の人が買うので，1億2700万人いる日本では約8,890万部ほど売れることになるな (なおサンプリングバイアス)。
twitter.com/bicycle1885/st...
タグ：
posted at 19:26:31
@watarukun1005 #Julia言語 いえいえどういたしまして。
π 以外にも e とか γ も同じ扱いです。
π = 3.1415926535897...
e = 2.7182818284590...
γ = 0.5772156649015...
π,eだけではなく、γ (Euler's γ constant)も入っていることがちょっとうれしい。任意桁を自動的に計算してくれる。 pic.twitter.com/vaB5MmCD3y
タグ： Julia言語
posted at 19:19:41
須山さんの本を読んで、 #Julia言語 を勉強したい場合には、須山さんの本に書いてある擬似コードを参考にして、Julia言語で書いてみればよいと思います。書き方はJulia言語のドキュメントを検索すればわかります。
docs.julialang.org/en/stable/
そこのSearch docsの利用の仕方をまず覚えるべき。
タグ： Julia言語
posted at 18:24:51
「ベイズ推論＝MCMCをまわすこと」のように誤解している人には須山さんの本は非常によいと思う。
手計算で得た有益な公式や適切な近似を使うと圧倒的にシンプルで高速な方法で必要な結果を得られることが説明されています。手計算ができない人はそういう美味しい道具を自分のものにできない。
タグ：
posted at 18:22:58
twitter.com/Dsuke_KATO/sta...
須山さんの『ベイズ推論による機械学習』のサンプルコード (#Julia言語)
github.com/sammy-suyama/B...
を理解したかったら、Julia言語よりも、本に書かれている数学的内容を理解することに重点をおくべきだと思います。数学が理解できればコードもすぐに理解できると思う。
タグ：
posted at 18:19:50
非公開
タグ：
posted at xx:xx:xx
twitter.com/yuki0626chan/s...
#Julia言語 を試しに使ってみることは Jupyter notebook が適していると思う。 juliabox.com なら即使える。
Jupyter notebook で
@ which factorial(im)
@ which gamma(im+1.0)
@ which lgamma(im+1.0)
(@ の後の空白は除く)
を順番に実行すれば～続く
タグ： Julia言語
posted at 18:09:45
非公開
タグ：
posted at xx:xx:xx
dockerでjupyter インストールした際に
juliaって言語入ってたんだけどすごいね
数式はほぼ数学の書き方どうりに定義できる
f(x)=xみたいな感じで
タグ：
posted at 14:25:22
@genkuroki 6個のコインは
ooo oooooo
ooo
のように長方形や一列に並べられるが、
5個だと
ooooo
一列にしか並べられない。
こんな説明もたまにします。
タグ：
posted at 13:46:34
#Julia言語 に興味を持って、いきなり重い数値計算をやらせようとしてしまうタイプの人は、まず
nbviewer.jupyter.org/gist/genkuroki...
Julia言語で計算が遅くなった場合の解決法
を見ておいた方が無難。特に配列に関する1.4節の内容。可読性を下げずに高速化する書き方を紹介しているつもり。
タグ： Julia言語
posted at 12:36:46
財務省の言い分は一貫してはいる。
大学に出す金そのものは減っていない。使い方が悪いから成果が出ない。勝手に決めさせるとろくな使い方をしないから、使い方に口は出す。
普通に考えれは、変に口を出すからパフォーマンスが下がるのだ、という結論になりそうなものなのだが。
タグ：
posted at 11:55:45
非公開
タグ：
posted at xx:xx:xx
twitter.com/xt_yuki/status...
#Julia言語 DifferentialEquationsパッケージは現時点ではupdateしない方が無難。詳しくはこのツイートが繋がっているスレッドを参照。
私が全パッケージを入れ直す手間をかけた結果の結論は
Pkg.pin("DifferentialEquations",v"3.1.0")
タグ： Julia言語
posted at 09:45:49
非公開
タグ：
posted at xx:xx:xx
JuliaでPyPlot.jl (matplotlib) を使ったヒストグラムのアニメーションを効率的に作成する
JuliaでもPyPlot.jlによってmatplotlibが使えるそうです。この記事はJuliaでアニメーションを作る話がきっかけになっていますが、実質は「matplotlibのplt.histあ… ift.tt/2E0L1yL
タグ：
posted at 02:30:34
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Perfectionatic @perfectionatic
@fermatslibrary I did a #julialang implementation of this a while back. See perfectionatic.org/?p=399
タグ： julialang
posted at 01:48:43
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@genkuroki matplotlibのマニュアルに書いてあったヒストグラムのアニメーションを効率よく描く方も加えてQiitaににまとめました。 qiita.com/skotaro/items/...
タグ：
posted at 00:06:35
