黒木玄 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月07日(土)
dragoner@2日目東イ14b @dragoner_JP
しかし、カズオ・イシグロの「心情的には日本人でいたかったけど、二重国籍を日本が許さなかったからイギリス国籍になった」という話が色々と悲しい → ノーベル文学賞 カズオ・イシグロが語った日本への思い、村上春樹のこと bunshun.jp/articles/-/4429
タグ：
posted at 00:16:05
#数楽 #JuliaLang 正規分布N(0,σ)の代わりに、ρ×t分布 ρt(ν) を使った線形回帰で遊んでみた。optim()函数の初期値依存の実験。
動画の作り方は次のリンク先に
nbviewer.jupyter.org/gist/genkuroki...
添付動画の収束先は最小二乗法と同じ。 pic.twitter.com/jbuIp4XKNz
posted at 04:01:59
#数楽 #JuliaLang 続き。自由度の大きなt分布は実質的に正規分布。
このツイートの添付動画の収束先は最小二乗法とは異なり、(最尤法としての)最小二乗法よりも尤度が大きくなる。動画では対数尤度の-1倍を計算。
optim函数はデフォルトでNelder-Meadを使う。 pic.twitter.com/cyowFlULim
posted at 04:04:57
#数楽 #JuliaLang パラメーター nu はt分布の自由度。 nu が100以上だとt分布はほぼ正規分布。ρは正規分布の標準偏差に対応するパラメーター。nuが大きなときには、ρは正規分布の標準偏差だと思ってよい。
posted at 04:07:54
#数楽 #JuliaLang 動画の動きが結構面白い。たぶん、Nelder-Meadの単体の動き方が回帰直線がふらふら動く様子として見えているのだと思いますが、本当にそうなのかどうかはわかりません。どうなんですかね？
posted at 04:12:51
#数楽 #JuliaLang 線形代数で解けるところを最適化函数optim()
julianlsolvers.github.io/Optim.jl/lates...
を使った最尤法で最小二乗法をやってみた。左右の「外れ値」のせいで傾きが0.2ではなく0.11に収束します。
Nelder-Mead法 pic.twitter.com/Fe94gwDWqK
posted at 11:29:01
#数楽 #JuliaLang ソースは
nbviewer.jupyter.org/gist/genkuroki...
にあります。
BFGS法 pic.twitter.com/KgtNzBukhU
posted at 11:30:43
#数楽 #JuliaLang
LBFGS法 pic.twitter.com/lYM7LuuOk3
posted at 11:31:34
#数楽 #JuliaLang
Conjugate Gradient法 pic.twitter.com/cakgFYOwUm
posted at 11:34:35
#数楽 #JuliaLang
Gradient Descent 法 pic.twitter.com/edgC7pws18
posted at 11:35:35
#数楽 #JuliaLang
Accelerated Gradient Descent 法 pic.twitter.com/Zk44DcXKYS
posted at 11:36:28
#数楽 #JuliaLang 動画的にはNelder-Mead法が一番面白いかな？
興味がある人向けリンク
bicycle1885.hatenablog.com/entry/2015/03/...
JuliaによるNelder-Meadアルゴリズムの実装
posted at 11:42:24
朝ごはんしながら、放送大学を見てたら、自然科学を学ぶ心構えみたいなことをやってて、数学の先生が「数学は自由」って言ってて、非現実的な事象も数学では組み立てられるって話なんだけど、やっぱりローカルルールを押し付けて教えてる超算数の話とか、どうしても思い出しちゃってねぇ…
タグ：
posted at 11:50:52
Dr. Chris Rackauckas @ChrisRackauckas
@genkuroki Would you be interested in contributing this to Optim.jl?
タグ：
posted at 12:07:23
面白すぎた。 -- ソ連の宇宙技術は最強過ぎたのだが、それを西側諸国が完全に理解したのはつい最近だった - ラヴェル船長の受難とソビエト及びロシアだけが保有する衛生管理技術 ncode.syosetu.com/n4343eh/1/
タグ：
posted at 13:51:25
マジか偏微分勉強しよう、という中学生も出てくる予感がしますが、数学に興味を持ってもらうという観点からはどちらでもOKですかね twitter.com/genkuroki/stat...
タグ：
posted at 15:05:37
内部留保課税に対して，正式な会計用語でないと言ってみても何の批判でもないし，二重課税だというのもピントがずれてる
でも，こういう微妙な反論になってしまいがちなくらい
内部留保課税わけわかんないんよ．何言ってるのか理解できないので芯を喰った批判もできないんじゃないかな？
タグ：
posted at 15:22:31
@E_Fessenden @sekibunnteisuu @OKWAVE 趣味の問題 減点せんでいい（英語やったらrtiangleやろ） #超理科
タグ： 超理科
posted at 16:06:58
@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
東工大の工大祭、ロボット技術研究会展示レポート。とても不思議な2枚歯の歯車。同じ向きに回転させることも出来ます。 pic.twitter.com/tr6hB2z3EX
タグ：
posted at 17:20:21
僕自身の経済状況がよくなった実感はまったくないけれども、学生の就職状況は明らかによくなっている。まあそういうことです。自分に直接の影響が及ばないからといって、経済状況が改善していないと考えるのは典型的な誤り
タグ：
posted at 18:21:43
自分の収入に関係ないけど、患者さんのコンプライアンスは改善しますね＞景気改善
どん底の頃は、インスリンを「節約」しちゃう人もいましたので…
（もちろん血糖管理は悪化しますが、お金がないという事で使用量減らしてた方がおられました…）
タグ：
posted at 18:25:24
#数楽 外れ値が出るような分布で残差を生成して作ったサンプルに最小二乗法をあてはめた結果。動画を面白くするために最尤法にしてoptim函数を使って解いた。初期値も動画を面白くするための設定にした。
サンプル生成のモデルは
y = 0.2 x + 0.1×TDist(0.8). pic.twitter.com/1eS3FvVd9q
タグ： 数楽
posted at 18:25:45
#数楽 動画はモデル y = a + bx + ρ×TDist(ν) に適用。外れ値を上手に処理しているように見える。νが大きなときこのモデルは標準偏差ρの通常の最小二乗法と同じ。ρとνが小さくなるときには、最尤法さんが「残差の裾野が重い」と判断していることを意味しています。 pic.twitter.com/y4CzztRyuw
タグ： 数楽
posted at 18:32:14
#数楽 #JuliaLang 最尤法の計算途中の動画作成のソースコードは以下の場所で公開。以下の前者にはまだ最新版が反映されていないが、じきに反映されるでしょう。
nbviewer.jupyter.org/gist/genkuroki...
gist.github.com/genkuroki/979c...
posted at 18:34:39
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
Dr. Chris Rackauckas @ChrisRackauckas
@genkuroki Ahh I see. github.com/JuliaNLSolvers... .
タグ：
posted at 23:06:14
