黒木玄 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日(日)
#数楽 確かにサンプルの上下のばらつき方を見ると外れ値はないように見える。
t分布の自由度nuの値がどんどん大きくなって行っているということから、最適化函数さんが「このサンプルには外れ値がないよ！！！」と叫んでいるように見える。
pic.twitter.com/A9V0c0F404
タグ： 数楽
posted at 23:56:01
#数楽 #JuliaLang 続き。自由度nuの大きなt分布は実質的に正規分布と区別がつかないので、添付画像の例で最尤法は正規分布になっていると判断していることになります。すなわち、最尤法さんは「外れ値はない」と判断している。
pic.twitter.com/A9V0c0F404
posted at 23:53:53
#数楽 続き。左半分の色付けは対数尤度函数の-1倍です。色が黒い部分ほど対数尤度の-1倍が小さいので、尤度は高い。最適化による最尤法の計算は△からスタートして□でストップしています。きちんと黒い部分に収束している。続く
pic.twitter.com/A9V0c0F404
タグ： 数楽
posted at 23:51:22
#数楽 #JuliaLang その確率密度函数は
(エラー) = ρ×T、 Tは自由度νのt分布に従う確率変数
というモデルの確率密度函数です。エラーのサンプルを与えて、ρとνの値を最尤法で計算しました。右半分がサンプル(縦軸がエラー、横軸は無意味な番号付け)。続く pic.twitter.com/A9V0c0F404
posted at 23:49:03
#数楽 #JuliaLang q(x|ν)をt分布の確率密度函数とし、確率密度函数
p(x|ρ,ν) = q(x/ρ|ν)/ρ
のサンプルサイズが小さい場合の尤度函数の様子をプロットしてみました。
nbviewer.jupyter.org/gist/genkuroki...
続く
posted at 23:43:49
便所は大切＞ソ連の宇宙技術は最強過ぎたのだが、それを西側諸国が完全に理解したのはつい最近だった『ラヴェル船長の受難とソビエト及びロシアだけが保有する衛生管理技術』 ncode.syosetu.com/n4343eh/1/ #narou #narouN4343EH
タグ： narou narouN4343EH
posted at 22:31:06
結果はなかなかおもしろい。きちっと検証できるツールを使ってファクトチェック。長野さんGJだね。
小池知事「核武装容認」サイトから削除？ 検証ツールで確認すると… - withnews（ウィズニュース） withnews.jp/article/f01710... #withnews
タグ： withnews
posted at 20:15:53
非公開
タグ：
posted at xx:xx:xx
眼しか動かせないALS患者の榊さんがOriHime eyeを使い、視線入力で描かれた最新作がこちら。
おそらく、眼だけでここまで絵をかける人はおらず、人類未踏の領域。人の可能性を更新しているのは痺れる。 pic.twitter.com/y6HEG7cmC8
タグ：
posted at 15:18:59
維新公認候補となる長谷川豊、去年夏、死刑囚について「ゴミはごみ箱に捨てるべきだ」「死刑にするなら、もっと残酷に殺すべき」「小学校時代からその死刑執行シーンはみんなに見せた方がいい」なんてことも書いてたんですね。維新こそゴミなのでは。news.ameba.jp/entry/20160727...
タグ：
posted at 11:24:06
非公開
タグ：
posted at xx:xx:xx
統合失調症を患った人を医療機関につなげる探偵会社。その社長・重川亮さんによると、「精神疾患のある人をカモにして成り立っている探偵会社は多く、業界の表沙汰にされたくない部分」だそうです。ライターの片岡健さんが取材しました。www.bengo4.com/other/n_6761/
タグ：
posted at 10:37:54
おお、これを使ったら本当にJuliaのコードをコンパイルして実行ファイルができたぞ。。 #julialang JuliaComputing/static-julia: Static AOT Julia compilation github.com/JuliaComputing...
タグ： julialang
posted at 09:59:17
中山成彬氏の「日教組の強いところは学力が低い」発言のファクトチェックを何度もした。最新作→ oku.edu.mie-u.ac.jp/~okumura/stat/... twitter.com/sspmi/status/9...
タグ：
posted at 08:32:31
#数楽 #JuliaLang 以上の2つは以下の2つの続きです。
gist.github.com/genkuroki/bc35...
測定スキルの異なる7人の科学者たち←外れ値で遊んでみる話
gist.github.com/genkuroki/1e1e...
動物の体重と脳の重さの関係
posted at 06:31:02
#数楽 #JuliaLang
nbviewer.jupyter.org/gist/genkuroki...
線形回帰のアニメーションの作り方
このファイルは生成されたアニメーション付き
Julia言語からmatplotlib.animationを使ってGIFアニメを作る方法もわかります。
posted at 06:28:29
#数楽 #JuliaLang
gist.github.com/genkuroki/d34e...
線形回帰のアニメーションの作り方
juliabox.com にアップロードして使える
出力をクリアしてアップロードできるサイズにおさえました。
posted at 06:25:30
非公開
タグ：
posted at xx:xx:xx
@sekibunnteisuu いくらなんでもテントウムシにしてしまったら鳳さん（与謝野晶子の実兄）がかわいそうで見てられない。
#超算数 #モルグリコ
posted at 00:06:37
2017年10月07日(土)
Dr. Chris Rackauckas @ChrisRackauckas
@genkuroki Ahh I see. github.com/JuliaNLSolvers... .
タグ：
posted at 23:06:14
Juliana LePine makes sculptures of celebrities, but hey, take a look at the amazing level of detials and anatomy youtu.be/L_vQW3xQEhc pic.twitter.com/jpNH5Bk4Mq
タグ：
posted at 21:02:15
#数楽 #JuliaLang 最尤法の計算途中の動画作成のソースコードは以下の場所で公開。以下の前者にはまだ最新版が反映されていないが、じきに反映されるでしょう。
nbviewer.jupyter.org/gist/genkuroki...
gist.github.com/genkuroki/979c...
posted at 18:34:39
#数楽 動画はモデル y = a + bx + ρ×TDist(ν) に適用。外れ値を上手に処理しているように見える。νが大きなときこのモデルは標準偏差ρの通常の最小二乗法と同じ。ρとνが小さくなるときには、最尤法さんが「残差の裾野が重い」と判断していることを意味しています。 pic.twitter.com/y4CzztRyuw
タグ： 数楽
posted at 18:32:14
#数楽 外れ値が出るような分布で残差を生成して作ったサンプルに最小二乗法をあてはめた結果。動画を面白くするために最尤法にしてoptim函数を使って解いた。初期値も動画を面白くするための設定にした。
サンプル生成のモデルは
y = 0.2 x + 0.1×TDist(0.8). pic.twitter.com/1eS3FvVd9q
タグ： 数楽
posted at 18:25:45
自分の収入に関係ないけど、患者さんのコンプライアンスは改善しますね＞景気改善
どん底の頃は、インスリンを「節約」しちゃう人もいましたので…
（もちろん血糖管理は悪化しますが、お金がないという事で使用量減らしてた方がおられました…）
タグ：
posted at 18:25:24
僕自身の経済状況がよくなった実感はまったくないけれども、学生の就職状況は明らかによくなっている。まあそういうことです。自分に直接の影響が及ばないからといって、経済状況が改善していないと考えるのは典型的な誤り
タグ：
posted at 18:21:43
東工大の工大祭、ロボット技術研究会展示レポート。とても不思議な2枚歯の歯車。同じ向きに回転させることも出来ます。 pic.twitter.com/tr6hB2z3EX
タグ：
posted at 17:20:21
@ChrisRackauckas The author of ANMS.jl already sent a pull request to the Optim.jl package.
See github.com/JuliaNLSolvers... and github.com/JuliaNLSolvers...
タグ：
posted at 16:44:05
@E_Fessenden @sekibunnteisuu @OKWAVE 趣味の問題 減点せんでいい（英語やったらrtiangleやろ） #超理科
タグ： 超理科
posted at 16:06:58
内部留保課税に対して，正式な会計用語でないと言ってみても何の批判でもないし，二重課税だというのもピントがずれてる
でも，こういう微妙な反論になってしまいがちなくらい
内部留保課税わけわかんないんよ．何言ってるのか理解できないので芯を喰った批判もできないんじゃないかな？
タグ：
posted at 15:22:31
マジか偏微分勉強しよう、という中学生も出てくる予感がしますが、数学に興味を持ってもらうという観点からはどちらでもOKですかね twitter.com/genkuroki/stat...
タグ：
posted at 15:05:37
面白すぎた。 -- ソ連の宇宙技術は最強過ぎたのだが、それを西側諸国が完全に理解したのはつい最近だった - ラヴェル船長の受難とソビエト及びロシアだけが保有する衛生管理技術 ncode.syosetu.com/n4343eh/1/
タグ：
posted at 13:51:25
Dr. Chris Rackauckas @ChrisRackauckas
@genkuroki Would you be interested in contributing this to Optim.jl?
タグ：
posted at 12:07:23
朝ごはんしながら、放送大学を見てたら、自然科学を学ぶ心構えみたいなことをやってて、数学の先生が「数学は自由」って言ってて、非現実的な事象も数学では組み立てられるって話なんだけど、やっぱりローカルルールを押し付けて教えてる超算数の話とか、どうしても思い出しちゃってねぇ…
タグ：
posted at 11:50:52
#数楽 #JuliaLang 動画的にはNelder-Mead法が一番面白いかな？
興味がある人向けリンク
bicycle1885.hatenablog.com/entry/2015/03/...
JuliaによるNelder-Meadアルゴリズムの実装
posted at 11:42:24
#数楽 #JuliaLang
Accelerated Gradient Descent 法 pic.twitter.com/Zk44DcXKYS
posted at 11:36:28
#数楽 #JuliaLang
Gradient Descent 法 pic.twitter.com/edgC7pws18
posted at 11:35:35
#数楽 #JuliaLang
Conjugate Gradient法 pic.twitter.com/cakgFYOwUm
posted at 11:34:35
#数楽 #JuliaLang
LBFGS法 pic.twitter.com/lYM7LuuOk3
posted at 11:31:34
#数楽 #JuliaLang ソースは
nbviewer.jupyter.org/gist/genkuroki...
にあります。
BFGS法 pic.twitter.com/KgtNzBukhU
posted at 11:30:43
#数楽 #JuliaLang 線形代数で解けるところを最適化函数optim()
julianlsolvers.github.io/Optim.jl/lates...
を使った最尤法で最小二乗法をやってみた。左右の「外れ値」のせいで傾きが0.2ではなく0.11に収束します。
Nelder-Mead法 pic.twitter.com/Fe94gwDWqK
posted at 11:29:01
#数楽 #JuliaLang 動画の動きが結構面白い。たぶん、Nelder-Meadの単体の動き方が回帰直線がふらふら動く様子として見えているのだと思いますが、本当にそうなのかどうかはわかりません。どうなんですかね？
posted at 04:12:51
#数楽 #JuliaLang パラメーター nu はt分布の自由度。 nu が100以上だとt分布はほぼ正規分布。ρは正規分布の標準偏差に対応するパラメーター。nuが大きなときには、ρは正規分布の標準偏差だと思ってよい。
posted at 04:07:54
#数楽 #JuliaLang 続き。自由度の大きなt分布は実質的に正規分布。
このツイートの添付動画の収束先は最小二乗法とは異なり、(最尤法としての)最小二乗法よりも尤度が大きくなる。動画では対数尤度の-1倍を計算。
optim函数はデフォルトでNelder-Meadを使う。 pic.twitter.com/cyowFlULim
posted at 04:04:57
#数楽 #JuliaLang 正規分布N(0,σ)の代わりに、ρ×t分布 ρt(ν) を使った線形回帰で遊んでみた。optim()函数の初期値依存の実験。
動画の作り方は次のリンク先に
nbviewer.jupyter.org/gist/genkuroki...
添付動画の収束先は最小二乗法と同じ。 pic.twitter.com/jbuIp4XKNz
posted at 04:01:59
dragoner@2日目東イ14b @dragoner_JP
しかし、カズオ・イシグロの「心情的には日本人でいたかったけど、二重国籍を日本が許さなかったからイギリス国籍になった」という話が色々と悲しい → ノーベル文学賞 カズオ・イシグロが語った日本への思い、村上春樹のこと bunshun.jp/articles/-/4429
タグ：
posted at 00:16:05
