黒木玄 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 を見て下さい。
2017年10月08日(日)
@sekibunnteisuu いくらなんでもテントウムシにしてしまったら鳳さん（与謝野晶子の実兄）がかわいそうで見てられない。
#超算数 #モルグリコ
posted at 00:06:37
非公開
タグ：
posted at xx:xx:xx
#数楽 #JuliaLang
gist.github.com/genkuroki/d34e...
線形回帰のアニメーションの作り方
juliabox.com にアップロードして使える
出力をクリアしてアップロードできるサイズにおさえました。
posted at 06:25:30
#数楽 #JuliaLang
nbviewer.jupyter.org/gist/genkuroki...
線形回帰のアニメーションの作り方
このファイルは生成されたアニメーション付き
Julia言語からmatplotlib.animationを使ってGIFアニメを作る方法もわかります。
posted at 06:28:29
#数楽 #JuliaLang 以上の2つは以下の2つの続きです。
gist.github.com/genkuroki/bc35...
測定スキルの異なる7人の科学者たち←外れ値で遊んでみる話
gist.github.com/genkuroki/1e1e...
動物の体重と脳の重さの関係
posted at 06:31:02
中山成彬氏の「日教組の強いところは学力が低い」発言のファクトチェックを何度もした。最新作→ oku.edu.mie-u.ac.jp/~okumura/stat/... twitter.com/sspmi/status/9...
タグ：
posted at 08:32:31
おお、これを使ったら本当にJuliaのコードをコンパイルして実行ファイルができたぞ。。 #julialang JuliaComputing/static-julia: Static AOT Julia compilation github.com/JuliaComputing...
タグ： julialang
posted at 09:59:17
統合失調症を患った人を医療機関につなげる探偵会社。その社長・重川亮さんによると、「精神疾患のある人をカモにして成り立っている探偵会社は多く、業界の表沙汰にされたくない部分」だそうです。ライターの片岡健さんが取材しました。www.bengo4.com/other/n_6761/
タグ：
posted at 10:37:54
非公開
タグ：
posted at xx:xx:xx
維新公認候補となる長谷川豊、去年夏、死刑囚について「ゴミはごみ箱に捨てるべきだ」「死刑にするなら、もっと残酷に殺すべき」「小学校時代からその死刑執行シーンはみんなに見せた方がいい」なんてことも書いてたんですね。維新こそゴミなのでは。news.ameba.jp/entry/20160727...
タグ：
posted at 11:24:06
眼しか動かせないALS患者の榊さんがOriHime eyeを使い、視線入力で描かれた最新作がこちら。
おそらく、眼だけでここまで絵をかける人はおらず、人類未踏の領域。人の可能性を更新しているのは痺れる。 pic.twitter.com/y6HEG7cmC8
タグ：
posted at 15:18:59
非公開
タグ：
posted at xx:xx:xx
結果はなかなかおもしろい。きちっと検証できるツールを使ってファクトチェック。長野さんGJだね。
小池知事「核武装容認」サイトから削除？ 検証ツールで確認すると… - withnews（ウィズニュース） withnews.jp/article/f01710... #withnews
タグ： withnews
posted at 20:15:53
便所は大切＞ソ連の宇宙技術は最強過ぎたのだが、それを西側諸国が完全に理解したのはつい最近だった『ラヴェル船長の受難とソビエト及びロシアだけが保有する衛生管理技術』 ncode.syosetu.com/n4343eh/1/ #narou #narouN4343EH
タグ： narou narouN4343EH
posted at 22:31:06
#数楽 #JuliaLang q(x|ν)をt分布の確率密度函数とし、確率密度函数
p(x|ρ,ν) = q(x/ρ|ν)/ρ
のサンプルサイズが小さい場合の尤度函数の様子をプロットしてみました。
nbviewer.jupyter.org/gist/genkuroki...
続く
posted at 23:43:49
#数楽 #JuliaLang その確率密度函数は
(エラー) = ρ×T、 Tは自由度νのt分布に従う確率変数
というモデルの確率密度函数です。エラーのサンプルを与えて、ρとνの値を最尤法で計算しました。右半分がサンプル(縦軸がエラー、横軸は無意味な番号付け)。続く pic.twitter.com/A9V0c0F404
posted at 23:49:03
#数楽 続き。左半分の色付けは対数尤度函数の-1倍です。色が黒い部分ほど対数尤度の-1倍が小さいので、尤度は高い。最適化による最尤法の計算は△からスタートして□でストップしています。きちんと黒い部分に収束している。続く
pic.twitter.com/A9V0c0F404
タグ： 数楽
posted at 23:51:22
#数楽 #JuliaLang 続き。自由度nuの大きなt分布は実質的に正規分布と区別がつかないので、添付画像の例で最尤法は正規分布になっていると判断していることになります。すなわち、最尤法さんは「外れ値はない」と判断している。
pic.twitter.com/A9V0c0F404
posted at 23:53:53
#数楽 確かにサンプルの上下のばらつき方を見ると外れ値はないように見える。
t分布の自由度nuの値がどんどん大きくなって行っているということから、最適化函数さんが「このサンプルには外れ値がないよ！！！」と叫んでいるように見える。
pic.twitter.com/A9V0c0F404
タグ： 数楽
posted at 23:56:01
2017年10月09日(月)
#数楽 特別に数学的に厳密な証明をする訓練を受けていなくても、以上のような計算実験をすれば、モデルの尤度函数の様子や最尤法で使用する最適化函数の「気持ち」がわかるようになると思う。
ソースファイル(汚いです！)は次の場所で全公開
nbviewer.jupyter.org/gist/genkuroki...
タグ： 数楽
posted at 00:01:17
#JuliaLang の Optim.jl パッケージのoptimize函数の値
r=optimize(f,init,Optim.Options(store_trace=true,extended_trace=true))
の中身は dump(r) すればわかります。続く
タグ： JuliaLang
posted at 00:10:18
#JuliaLang 続き。トレースのオプションをオンにしておくと、
r.trace
に途中の計算の様子の記録が格納されています。その中身も dump(r.trace) すればわかります。
r.trace[10].metadata
に10番目の計算結果が入っている。
タグ： JuliaLang
posted at 00:12:28
#JuliaLang 続き。計算途中のパラメーターの値は、デフォルトのNelder-Mead法では
r.trace[i].metadata["centroid]
に入っており、他のsolverでは
t.trace[i].metadata["x"]
に入っています。
タグ： JuliaLang
posted at 00:16:00
#JuliaLang の使い方を覚えるコツは、
?hoge
でドキュメントを検索しまくったり、
dump(hoge)
で hoge の中身を覗いてみたりすることだと思います。あと
{at-mark}time hoge
でhogeの実行時間とメモリ使用量を常にチェック。
タグ： JuliaLang
posted at 00:19:33
現時点で、枝野と前原なら、枝野支持するし
立憲民主党と希望党なら、立憲民主党支持するけど
民主党が政権取る時に、少しだけ期待してた時の事を思い出す。
子ども手当の公約や、小沢鋭仁、馬淵澄夫、金子洋一に期待したんだ。
結局、裏切られた。民主党では財務官僚と戦えない。
タグ：
posted at 00:38:45
#数楽 #JuliaLang 上のt分布モデルを線形回帰に利用する場合は対数尤度函数の-1倍の値を大きなν(>10^8)と小さなρ(<0.1^4)で∞に設定しておくと、挙動が安定し易くなります。大きなνは正規分布を意味するので、そうしておいてもたぶん問題ない。
posted at 01:02:22
先生もくっつけちゃったから、わざわざ修正テープで修正。教員は忙しいと言いながら、こんなどうでもいい事に時間使ってるんだなぁ、と思ってしまう。くっつくかくっつかないかは「フォントの差」です。正誤には無関係。ほら、「夜」、くっついてる。 pic.twitter.com/lR0Cp2uGWf
タグ：
posted at 01:53:19
@genkuroki 『テコンドー韓流ダイエット 岡本依子のネリチャギビューティー』
とかいう物が実在するので、どちらの観測も当たっているのでは。
（予告編 www.youtube.com/watch?v=gpxFKi...）
タグ：
posted at 02:02:23
Distributing points on a sphere
www.maths.unsw.edu.au/about/distribu... pic.twitter.com/FTiOgMKUDi
タグ：
posted at 03:52:49
渡辺先生のベイズ統計本で、ようやくWAIC導出までの式が追えた。ミソは、特異点解消してδ超関数を漸近展開するところ。メリン変換後の極の最大位数mがlog 0型の漸近挙動を招くことで、自由エネルギーFの漸近展開にlog(log(n))なるクソキモい項が現れる。極の値-λは、Fにlog(n)の寄与を与える。
タグ：
posted at 05:11:09
渡辺先生のベイズ本4章を読み始めた頃は、「統計力学やったことあるし、まぁ簡単にいけるやろ」などとタカをくくっていたが、見積もり甘すぎた。WAICの導出は2010年らしいし、そんなに簡単なわけないよな。でも何事につけ「いけるやろ」という心境でトライするのは重要だと思う。
タグ：
posted at 05:34:16
#数楽 #JuliaLang
gist.github.com/genkuroki/469e...
y = bx + ρT(ν)、T(ν)は自由度νのt分布に従う
というモデル(パラメーターはb,ρ,ν)の尤度函数で真のミニマム以外にローカルミニマムを持つ例のJupyter notebook
posted at 11:35:21
#数楽 リンク先の動画がローカルミニマムで止まる例。外れ値を無視できずに最小二乗法解に落ち込んでいる。
twitter.com/genkuroki/stat...
タグ： 数楽
posted at 12:30:48
JuMPを使ったナップサック最適化問題をこれでコンパイルしたら、コンパイル込みで11秒かかっていたコードが、0.8秒で実行できた。。。すごいJulia最強ではないか。。 #julialang
タグ： julialang
posted at 14:25:37
#数楽 真の解と別にローカルミニマムがある小さな例
model: y = rho*TDist(nu)
sample: [0.01, 1.0]
初期条件 (rho,nu)=(1.0,1.0)だとローカルミニマムに落ちて行ってしまう。(0.9,1.0)にずらすと真の解へ。 pic.twitter.com/eUsc0w5tUm
タグ： 数楽
posted at 14:34:36
#数楽 #JuliaLang 一つ前のツイートのサイズ2のサンプルで最尤法の真の解とは異なるローカルミニマムがあるケースのプロットのソースコードは
nbviewer.jupyter.org/gist/genkuroki...
の下の方にあります。
posted at 14:41:50
#数楽 自由度νのt分布はν→∞で標準正規分布に収束します。
添付画像は確率密度函数の対数のグラフです。ν=100のt分布は±5σ程度まで標準正規分布に一致し、ν=1000なら±10σ程度まで、ν=10000なら±20σ程度まで一致します。上のローカルミニマムでは ν＝1億！ pic.twitter.com/FaDKIXCoDh
タグ： 数楽
posted at 15:03:06
先ほどの資料です。
deeplearn.js入門
docs.google.com/presentation/d... #devfest17
タグ： devfest17
posted at 15:03:17
非公開
タグ：
posted at xx:xx:xx
小学校の先生しか分からない苦労選手権
最優秀賞
問題は生徒対応より保護者対応。
金賞
高学年女子の陰口の恐ろしさ
入選
1＋1という名の哲学講義
男の子は大体瞬足をはいているので一緒に走るとコーナーで差をつけられる
キラキラネームが解読不可能 pic.twitter.com/6XafhP7PkV
タグ：
posted at 15:28:51
#数楽 #JuliaLang のDistribution.jlパッケージのソースコードを見て(ドキュメント
juliastats.github.io/Distributions....
からソースコードにリンクがはってあるので必然的に見ることになる)、ちょっと勉強になったこともあります。それは～続く
posted at 15:32:13
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#数楽 続き。よく出て来る正規分布とガンマ分布が自然に出て来る中心極限定理とは別の仕組みについては以下の2つのリンク先を参照。
正規分布↓
twitter.com/genkuroki/stat...
ガンマ分布↓
twitter.com/genkuroki/stat...
タグ： 数楽
posted at 15:57:31
非公開
タグ：
posted at xx:xx:xx
4×6と6×4は答えは一致するが、式の意味が違う。4が6個か6が4個か。その違いを説明出来ない人が多い。2年生の掛け算の導入で習ったはず。
式の順序の大切さは「行列」(旧数学Ｃ，現在は高校で扱わない)を習うまでイマイチ理解しにくい。そんな算数・数学の現状をどうにか出来ないものか。 twitter.com/bozu_108/statu...
タグ：
posted at 19:52:34
link: C言語でインクルードするだけで使えるNon-movingで正確なコピーGCを作った - Qiita:
qiita.com/wasabiz/items/...
タグ：
posted at 20:46:52
link: nyuichi/gc.h: Header file-only Non-moving & Precise GC for C (< 100LOC w/o comments):
github.com/nyuichi/gc.h
タグ：
posted at 20:48:55
UFO教授 (藤木文彦 Fumihiko @UFOprofessor
@genkuroki ちゃんとした説明のしてある本があったら私も教材に使いたいです。去年苦労したので。
タグ：
posted at 21:06:03
#数楽 #JuliaLang
下のリンク先のツイートのモデル(と同値な階層モデル)とサンプルでベイズ推定してみました。Julia言語でのMCMC適用の非常に小さな実例です。
nbviewer.jupyter.org/gist/genkuroki...
twitter.com/genkuroki/stat...
posted at 21:37:41
#数楽 #JuliaLang のMCMCパッケージのMamba.jlは複数のチェインを並列処理してくれます。添付画像は同時に作った8つのチェインの最後の一つです。事前分布は対数座標での正規分布。
詳しくはソースコードを参照。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/2yZpMMKr1z
posted at 21:40:38
#数楽 #JuliaLang なぜか事前分布として最初手抜きの一様分布を使っていたのだが、それを正規分布に変えたら、計算がかなり速くなった。Mamba.jlについては全然理解していません。
mambajl.readthedocs.io/en/latest/
posted at 21:46:03
魔法少女アニメ用音響エフェクトの自動生成に関するNicographの論文。Akinori Ito は私ではない。微妙に近い分野で同姓同名がいるとまぎらわしいなあ ieeexplore.ieee.org/abstract/docum...
タグ：
posted at 21:54:30
UFO教授 (藤木文彦 Fumihiko @UFOprofessor
@genkuroki R入門の本が余りにひどかったので、某社に、
書きましょうかと提案しているところです。
タグ：
posted at 22:31:27
置き場所が変わったみたいなので、ＵＲＬを更新する。
大学の数学を苦痛と捉える学生の実態把握と、その原因の探求
～小・中学校数学専攻 4 回生を対象としたアンケート調査を通して～
satsuki.ex.osaka-kyoiku.ac.jp/~j169302/miyaz...
タグ：
posted at 23:20:48
Bayesian HMM、Infinite HMM、無限木構造HMMによる教師なし品詞タグ推定の実装をすべて書き直し誤りを修正しました。iTHMMは実験中ですが、深さを固定した時にIHMMと同様の結果になるのでバグはなさそう。
github.com/musyoku/unsupe...
タグ：
posted at 23:28:58
