黒木玄 Gen Kuroki(@genkuroki)/2017年10月07日

黒木玄 Gen Kuroki

@genkuroki

いいね数 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 を見て下さい。

Favolog ホーム » @genkuroki » 2017年10月07日

« 次の日| 前の日 »

並び順 : 新→古 | 古→新

2017年10月07日(土)

Dr. Chris Rackauckas @ChrisRackauckas

17年10月7日

@genkuroki Ahh I see. github.com/JuliaNLSolvers... .

タグ：

posted at 23:06:14

Massimo @Rainmaker1973

17年10月7日

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

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽

t分布×定数のモデル：私の正体がわかるかね？

数楽探偵：おまえの正体は正規分布の分散が逆ガンマ分布に従う階層モデルの一味だな！　測定値ごとに分散に「個性」があるモデルだ！

t分布×定数のモデル：よくぞ見抜いたな！それではさらばだ！

数楽探偵：いかないで！😅

タグ：数楽

posted at 18:47:36

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽

ユーザー：なあなあ、最小二乗法さん、あっちのやつは外れ値だから、無視して、きちんとフィッティングしてよ。

最小二乗法：そんなこと言われても困るよ。ぼくは「みんな等質だ」と思ってがんばってフィッティングしているんだから。外れ値のような「個性」を考えたいならまた別の……略

タグ：数楽

posted at 18:41:48

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽「最適化」の計算の途中の様子を動画にすると結構面白い。動画を見ると

こうかな？ちょっと違う。こっちの方がいいかな？これでどうだ？……

のような気分で計算しているというな「気持ち」を感じてしまう。

タグ：数楽

posted at 18:37:15

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 最尤法の計算途中の動画作成のソースコードは以下の場所で公開。以下の前者にはまだ最新版が反映されていないが、じきに反映されるでしょう。

nbviewer.jupyter.org/gist/genkuroki...

gist.github.com/genkuroki/979c...

タグ： JuliaLang 数楽

posted at 18:34:39

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽動画はモデル y = a + bx + ρ×TDist(ν) に適用。外れ値を上手に処理しているように見える。νが大きなときこのモデルは標準偏差ρの通常の最小二乗法と同じ。ρとνが小さくなるときには、最尤法さんが「残差の裾野が重い」と判断していることを意味しています。 pic.twitter.com/y4CzztRyuw

タグ：数楽

posted at 18:32:14

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽外れ値が出るような分布で残差を生成して作ったサンプルに最小二乗法をあてはめた結果。動画を面白くするために最尤法にしてoptim函数を使って解いた。初期値も動画を面白くするための設定にした。

サンプル生成のモデルは

y = 0.2 x + 0.1×TDist(0.8). pic.twitter.com/1eS3FvVd9q

タグ：数楽

posted at 18:25:45

KGN @KGN_works

17年10月7日

自分の収入に関係ないけど、患者さんのコンプライアンスは改善しますね＞景気改善

どん底の頃は、インスリンを「節約」しちゃう人もいましたので…
（もちろん血糖管理は悪化しますが、お金がないという事で使用量減らしてた方がおられました…）

タグ：

posted at 18:25:24

あ〜る菊池誠(反緊縮)公式 @kikumaco

17年10月7日

僕自身の経済状況がよくなった実感はまったくないけれども、学生の就職状況は明らかによくなっている。まあそういうことです。自分に直接の影響が及ばないからといって、経済状況が改善していないと考えるのは典型的な誤り

タグ：

posted at 18:21:43

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽平方完成について理解していれば最小二乗法を初等的に自力で思い付くことが可能。基本的な考え方の組み合わせ。

昔の偉い人が考えたことは、現代のちょっと数学を知っている人たちが自明に思い付くアイデアに過ぎなかったりする。これこそが数学の抽象化による平易化と数学教育普及の力。

タグ：数楽

posted at 18:19:32

KenKenMkIISR @KenKenMkIISR

17年10月7日

東工大の工大祭、ロボット技術研究会展示レポート。とても不思議な2枚歯の歯車。同じ向きに回転させることも出来ます。 pic.twitter.com/tr6hB2z3EX

タグ：

posted at 17:20:21

黒木玄 Gen Kuroki @genkuroki

17年10月7日

@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

A級３班国民 @kankichi573

17年10月7日

@E_Fessenden @sekibunnteisuu @OKWAVE 趣味の問題　減点せんでいい（英語やったらrtiangleやろ） #超理科

タグ：超理科

posted at 16:06:58

飯田泰之 @iida_yasuyuki

17年10月7日

内部留保課税に対して，正式な会計用語でないと言ってみても何の批判でもないし，二重課税だというのもピントがずれてる

でも，こういう微妙な反論になってしまいがちなくらい

内部留保課税わけわかんないんよ．何言ってるのか理解できないので芯を喰った批判もできないんじゃないかな？

タグ：

posted at 15:22:31

elcondor @elcondor

17年10月7日

マジか偏微分勉強しよう、という中学生も出てくる予感がしますが、数学に興味を持ってもらうという観点からはどちらでもOKですかね twitter.com/genkuroki/stat...

タグ：

posted at 15:05:37

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽平方完成が難しいと感じる中学生には、「最小二乗法を高級な偏微分を使って導いている理系の大学生のお姉さんやお兄さんたちは、もしかしたら、平方完成は難しいという理由で避けているのかもしれない」なんて話をすれば、「お姉さんやお兄さんよりも上に行こう」と思ってくれたりしないかな？

タグ：数楽

posted at 13:53:52

xmisao / ぺけみさお @xmisao

17年10月7日

面白すぎた。 -- ソ連の宇宙技術は最強過ぎたのだが、それを西側諸国が完全に理解したのはつい最近だった - ラヴェル船長の受難とソビエト及びロシアだけが保有する衛生管理技術 ncode.syosetu.com/n4343eh/1/

タグ：

posted at 13:51:25

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽最小二乗法の公式の導出は、中学校レベルの平方完成で可能なのに、わざわざ偏微分を使うのは簡単なことを難しく解いている感じがしないでもない。しかし、平方完成の計算は「難しいと感じられるような計算」だということなのかもしれない。

タグ：数楽

posted at 13:50:30

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽たぶん、平方完成は中学校で習うもっとも複雑な計算のテクニックだと思います。「こんなことをやって何の役に立つのか」と思う人は多いでしょうが、最小二乗法くらい知っておいてばちは当たらないと思う。普通に教えられている数学で役に立たないものは皆無。

タグ：数楽

posted at 13:45:31

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽最小二乗法の解説では、f(x) = x^2-2x を最小にするxを「f'(x)=2x-2=0 なので x=1」と求める方法の一般化が説明されていることが多いと思う。「f(x)=(x-1)^2-1 なので x=1」の一般化も同時に説明して欲しいところだと思います。

タグ：数楽

posted at 13:43:35

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽最小二乗法の解の公式を、偏導函数を0とおいて解くことによって得ている解説は多いが、中学校レベルの平方完成

ax^2 + bx = a(x-bx/(2a))^2 - b^2/(4a)

の一般化で初等的に出すこともできる。

タグ：数楽

posted at 13:41:24

Dr. Chris Rackauckas @ChrisRackauckas

17年10月7日

@genkuroki Would you be interested in contributing this to Optim.jl?

タグ：

posted at 12:07:23

らじうむ小山_PPPMP @Ra_koyama

17年10月7日

朝ごはんしながら、放送大学を見てたら、自然科学を学ぶ心構えみたいなことをやってて、数学の先生が「数学は自由」って言ってて、非現実的な事象も数学では組み立てられるって話なんだけど、やっぱりローカルルールを押し付けて教えてる超算数の話とか、どうしても思い出しちゃってねぇ…

タグ：

posted at 11:50:52

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 動画的にはNelder-Mead法が一番面白いかな？

興味がある人向けリンク

bicycle1885.hatenablog.com/entry/2015/03/...
JuliaによるNelder-Meadアルゴリズムの実装

タグ： JuliaLang 数楽

posted at 11:42:24

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang

Accelerated Gradient Descent 法 pic.twitter.com/Zk44DcXKYS

タグ： JuliaLang 数楽

posted at 11:36:28

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang

Gradient Descent 法 pic.twitter.com/edgC7pws18

タグ： JuliaLang 数楽

posted at 11:35:35

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang

Conjugate Gradient法 pic.twitter.com/cakgFYOwUm

タグ： JuliaLang 数楽

posted at 11:34:35

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang

LBFGS法 pic.twitter.com/lYM7LuuOk3

タグ： JuliaLang 数楽

posted at 11:31:34

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang ソースは

nbviewer.jupyter.org/gist/genkuroki...

にあります。

BFGS法 pic.twitter.com/KgtNzBukhU

タグ： JuliaLang 数楽

posted at 11:30:43

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 線形代数で解けるところを最適化函数optim()

julianlsolvers.github.io/Optim.jl/lates...

を使った最尤法で最小二乗法をやってみた。左右の「外れ値」のせいで傾きが0.2ではなく0.11に収束します。

Nelder-Mead法 pic.twitter.com/Fe94gwDWqK

タグ： JuliaLang 数楽

posted at 11:29:01

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 動画の動きが結構面白い。たぶん、Nelder-Meadの単体の動き方が回帰直線がふらふら動く様子として見えているのだと思いますが、本当にそうなのかどうかはわかりません。どうなんですかね？

タグ： JuliaLang 数楽

posted at 04:12:51

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang パラメーター nu はt分布の自由度。 nu が100以上だとt分布はほぼ正規分布。ρは正規分布の標準偏差に対応するパラメーター。nuが大きなときには、ρは正規分布の標準偏差だと思ってよい。

タグ： JuliaLang 数楽

posted at 04:07:54

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 続き。自由度の大きなt分布は実質的に正規分布。

このツイートの添付動画の収束先は最小二乗法とは異なり、(最尤法としての)最小二乗法よりも尤度が大きくなる。動画では対数尤度の-1倍を計算。

optim函数はデフォルトでNelder-Meadを使う。 pic.twitter.com/cyowFlULim

タグ： JuliaLang 数楽

posted at 04:04:57

黒木玄 Gen Kuroki @genkuroki

17年10月7日

#数楽 #JuliaLang 正規分布N(0,σ)の代わりに、ρ×t分布 ρt(ν) を使った線形回帰で遊んでみた。optim()函数の初期値依存の実験。

動画の作り方は次のリンク先に
nbviewer.jupyter.org/gist/genkuroki...

添付動画の収束先は最小二乗法と同じ。 pic.twitter.com/jbuIp4XKNz

タグ： JuliaLang 数楽

posted at 04:01:59