黒木玄 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 を見て下さい。
2020年07月02日(木)
@martin_trapp @PyTorch @soumithchintala @TuringLang @Hong_Ge_ Yup! It has also identified the gaps in the #julialang compiler - so that we can achieve the same or higher performance with native #julialang going forward. For small problems, CUDA.jl is within a factor of 2, and for larger ones, we need to fix GC.
タグ: julialang
posted at 00:33:19
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
Dr. Chris Rackauckas @ChrisRackauckas
@Raygun0001 @aureliengeron @hazemessamm If you're looking for some benchmarks related to this, look at DiffEqBenchmarks.jl github.com/SciML/DiffEqBe... . Generally there's more than an order of magnitude performance advantage for #julialang here: benchmarks.sciml.ai/html/MultiLang... .
タグ: julialang
posted at 01:56:59
Dr. Chris Rackauckas @ChrisRackauckas
@Raygun0001 @aureliengeron @hazemessamm and we have some benchmarks that are specific to high precision energy conservation, i.e. benchmarks.sciml.ai/html/Dynamical... and benchmarks.sciml.ai/html/Dynamical...
タグ:
posted at 01:57:46
One thing about #julialang for Statistics/Machine Learning/AI that doesn't get enough mention is how multiple dispatch is such a natural formulation for high performance numerical linear algebra code.
タグ: julialang
posted at 04:53:41
ごまふあざらし(GomahuAzaras @MathSorcerer
@Viral_B_Shah I could build Julia 1.x for Arm32 bit system using jlcross.
github.com/Julia-Embedded...
タグ:
posted at 06:07:51
ごまふあざらし(GomahuAzaras @MathSorcerer
@Viral_B_Shah Here is my summary what can Julia on RPi series do.
discourse.julialang.org/t/have-a-try-j...
タグ:
posted at 06:10:41
Nirmal Suthar @nirmal_p_suthar
I am delighted to share my blog post demonstrating 3D mesh reconstruction task. This blog covers my progress in GSoC project this year. Give it a read :) #julialang #gsoc
nextjournal.com/nirmal-suthar/...
posted at 06:21:37
@kamo_hiroyasu @tegosshy @1027stesc twitter.com/sekibunnteisuu... 教科書が、恒等式と方程式をこんな具合に区別している。
|a|=a は、
aの範囲を0以上の実数とするなら、「恒等式であって、方程式ではない」
すべての実数の範囲とするなら、「方程式であって、恒等式ではない」
なんとも不合理。
タグ:
posted at 07:12:28
#julialang plotting syntax feels great;
but then when I look at the actual content.....feels bad pic.twitter.com/NwWTMF3ZnL
タグ: julialang
posted at 09:11:27
結局、西村大臣はこういう無責任メディアにつけ入る隙を与えたということ。「感謝している」と言いながら、実際にやっていることは、専門家たちを感染症対策の矢面に立たせ非常に危険な場所に登らせるようなことをさせておきながら突然梯子を外したに等しい。
dot.asahi.com/wa/20200630000...
タグ:
posted at 12:38:00
#統計 MCMCで「そのような」擬似乱数列を作る様子の例。
nbviewer.jupyter.org/gist/genkuroki...
自前実装のGibbsサンプリング。
twitter.com/genkuroki/stat... pic.twitter.com/gyezVKRCrp
タグ: 統計
posted at 13:27:36
#統計 みんな大好きなNUTSによる擬似乱数の生成
nbviewer.jupyter.org/gist/genkuroki...
ちょっと古めの #Julia言語 Turing.jl を使っています。現在の版で動かすには微小な手直しが必要。 pic.twitter.com/Nou9QHJiw2
posted at 13:27:46
2日の夕刊文化面は「掛け算の順序問題」について、Twitter上で積極的に議論を重ねておられる「積分定数」氏こと、西沢宏明さんに、その問題点についてご寄稿いただきました。(典) pic.twitter.com/o9U1Q5S75Z
タグ:
posted at 14:17:08
「掛け算の順序問題」とは、自然数の掛け算は交換法則が成り立つため、5×4でも4×5でも答えは同じになるにもかかわらず、学び始めに「式に順序がある」という指導法が広く行われ、5×4は正答でも4×5は誤答とされるなどのケースがあり、その是非について長年、論争が行われている問題です。(典)
タグ:
posted at 14:19:06
@iida01 @bicycle1885 Radon-Nikodym微分は密度函数p(x),q(x)で確率測度が書けているケースではp(x)/q(x)になるので、「Radon-Nikodym微分」と言う必要はないです。
その辺の事情が明瞭でないと、測度論的なRadon-Nikodym微分が本質的に必要かどうかは不明だとするべきだと思います。
twitter.com/genkuroki/stat...
タグ:
posted at 16:44:41
@iida01 @bicycle1885 いずれにせよ、密度函数で確率測度が書けている場合で密度函数の比が必要になる場合を知らない人が、いきなり測度論的なRadon-Nikodym微分について勉強させられることは、24÷6のような計算さえしたことがない人に整数論の勉強をさせるようなもので、良くないと思います。
タグ:
posted at 16:47:44
Stan3, backendにtensorflow probabilityが選べるようになるんだね、びびった。たしかに最尤推定はStanよりtensorflowの方が良さそうと思ってる。MCMCはStanの方が性能がよさげ。変分ベイズは分からない。
github.com/stan-dev/stanc...
タグ:
posted at 17:09:47
OokuboTact 大久保中二病中年 @OokuboTact
#超算数
周回遅れの人のためのQ&A集が必要だと思う。
もぐら叩きは面倒
twitter.com/sekibunnteisuu...
タグ: 超算数
posted at 18:01:45
OKUMURA, Akira(奥村 曉) @AkiraOkumura
Zoom での講義は iPad で PDF を表示している状態を Zoom の AirPlay の共有機能を使うのが良いかなと思ったのだが、肝心の iPad 側でプレゼンに適した PDF viewer が見つからない。余計なツールバーは表示されず、指でページ遷移、Apple Pencil で書き込み、が固定できるようなもの。何かないですか?
タグ:
posted at 18:45:44
よくご存じですね。8月1日から研究員一同を引き連れて京都大学大学院医学研究科(社会健康医学系専攻)へ異動します。昨年11月に決まっていたので北大では周りの教授が「あいつ本当に引っ越すのか」ってずっと心配してくれてました。現在、引っ越し日程と流行分析のはざまであせあせしています。 twitter.com/anesthyuriko/s...
タグ:
posted at 18:50:44
Just a small #JuliaLang version of one of the plots by @standupmaths in his new video:
www.youtube.com/watch?v=ghxQA3...
Code: gist.github.com/Wikunia/dd53f8...
P.S. I just like \xi and I don't know whether there is a standard for (1-sqrt(5)/2) Yes I could use -1/\phi :D pic.twitter.com/HFLgiz586o
タグ: JuliaLang
posted at 19:02:15
OKUMURA, Akira(奥村 曉) @AkiraOkumura
@Hal_Tasaki @psr1509 @Hal_Tasaki ありがとうございます。この機能ですかね。is4-ssl.mzstatic.com/image/thumb/Pu...
タグ:
posted at 19:08:41
@opensourcesblog @standupmaths It's weird to see special Unicode characters like Greek or a sqrt symbol on real practical code 😅
タグ:
posted at 19:11:37
@AkiraOkumura @psr1509 そうですね。ぼくは最初は知らずにツールバーを表示してたんですが・・
あと、残像が残る赤いポインターもなかなかよいです。
タグ:
posted at 19:13:08
@opensourcesblog @standupmaths Do you write this Unicode characters by yourself? If so, how? Googling them and then copypasting into the code? 😳
タグ:
posted at 19:17:40
@Sh1bumi @standupmaths Just type the LaTeX code and press tab. So i.e. \sqrt[TAB] works in the REPL and supported editors like @code with the #JuliaLang extension
タグ: JuliaLang
posted at 19:19:19
@Hal_Tasaki @AkiraOkumura 私は偶然(?)Zoomで録画するとツールバーが収録されない設定(?)でした。気にしたことがなかったです。
タグ:
posted at 19:22:27
@Sh1bumi @standupmaths @code Funfact: If you don't know the name but somehow found it only you can copy it into the REPL help mode to see how to write it.
#julialang #unicode pic.twitter.com/eo254PoYKT
posted at 19:25:04
OKUMURA, Akira(奥村 曉) @AkiraOkumura
@psr1509 @Hal_Tasaki 学生の iPad に入っている GoodNotes を触らせてもらって(もうコロナ気にしてない…)もよく分からなかったんですが、Mac 側にミラーリングをしないとプレゼンテーションモードにならないということですか?
タグ:
posted at 19:30:04
新型コロナ患者を受け入れた医療機関が赤字になったり、医療者の給料がカットされたりするのはあまりにも理不尽なので、国が補償するべきだと思いますよ。
拍手だの感謝の言葉だのよりお金のほうが圧倒的にいいはずです
タグ:
posted at 19:31:16
#掛算 の順序問題は、「学び始め」ではないのです。
交換法則は2年生でかけ算を導入後すぐに指導されますが、「かけ算の順序を指導しないと3年生以降のわり算の理解に支障が出る」などという「算数教育の理論」が出てきて、3年生になってからも順序指導が行われるのです。
#超算数 twitter.com/tokyobunkabu/s...
posted at 19:33:50
@AkiraOkumura @psr1509 いいえ。ツールバーにある四角の真ん中から上向の矢印が生えてるやつをタップすると「共有」というメニューが開きます。この下三つにプレゼンテーションモードの選択があります。これで一番したの「ページ全体をミラーリング」を選択します。
タグ:
posted at 19:36:53
@tokyobunkabu 掛算順序問題はいつも大論争になりますが、エンジニアや科学者など、日常的に数学を使う大人の大部分が小学校での教え方に否定的です。
twitter.com/s_kajita/statu...
タグ:
posted at 19:38:17
#統計 添付画像1の最下段の注意を参照。
しかし、多項分布におけるPearsonのχ²統計量
Σ(観測量 - 期待値)²/期待値
はPoisson分布の直積の条件付き確率分布で多項分布を作れることから来ていることも理解しておいた方がお得。KL情報量が自然に出て来る。
分割表のχ²検定ではWilksの定理が重要。 twitter.com/wata0303jug/st... pic.twitter.com/18SqchuM0M
タグ: 統計
posted at 19:38:46
OKUMURA, Akira(奥村 曉) @AkiraOkumura
@Hal_Tasaki @psr1509 そこを選択して、Zoom や Mac はどういう操作をするのでしょうか。iPad 側で Zoom をホストするのですか?
タグ:
posted at 19:41:46
@AkiraOkumura @psr1509 そういう流儀もあるようですが、ぼくは
Mac で Zoom を走らせる
Zoom の share screen で iPad/iPhone via airplay を選ぶ
iPad 側で(ツールバーを引っ張り下ろしたりして)画面共有を許可
というふうにしています。
線でも繋げるのですが、なんかよくわからずもっぱら air play です。
タグ:
posted at 19:45:03
自分の中では結構画期的な発見。
新型コロナウイルスの再生算数を求める問題、発表・発症日から感染日を逆算するために、西浦先生含めてみんな逆畳み込みをやっているけど、この方法では数字が暴れてスムージング等が必要になる。
しかし、そういう人工的な操作はあまり入れたくない。どうすべきか。
タグ:
posted at 20:06:09
そもそも、難しいことを考えずとも、素直に患者ごとの発表・発症日から「感染から発表・発症までの日数分布に基づいて発生させた乱数」を引き算するだけで良いんじゃないか?と思い、やってみたのがこれ。
東京都の感染者数データや描画を含め、Mathematicaでこれだけ。 pic.twitter.com/u7hlObSjNG
タグ:
posted at 20:06:11
算出時点で未報告の患者を無視しているので最後が下がっているが、あっさり上手くいってびっくり。患者数6405人×モンテカルロ数分の乱数を発生させる必要があるが、独立なのでMCMC等を使う必要がなく、50000モンテカルロ並列化無しでも40秒程で完了。
タグ:
posted at 20:06:11
Rtは、Serial Intervalに応じて発生させた乱数を感染日に足して日ごとのデータに戻し、感染者数をそれで割るだけ。ここまで一切スムージング等の操作なし。
1つ1つのRt時系列は多少暴れているので、スムーズにするような感染日組合せを求めるのも良いかもしれない。その場合にはMCMCが必要になる。 pic.twitter.com/DnSiSST3GD
タグ:
posted at 20:06:13
未報告患者分を補正したものも作ってみた。日ごとの報告率で感染・2次感染者数を割る。ラスト3日分はカットしたが、2週間分ぐらいは信頼できないかも。
東京都は3月上旬にRt値のピークを迎え、4月上旬に1を切り、連休中に最低値、緊急事態解除につれて上昇し1.1~1.3ぐらいで安定している感じだなあ。 pic.twitter.com/nXk9q8WGd3
タグ:
posted at 20:06:15
#統計 二項分布の視覚化
nbviewer.jupyter.org/gist/genkuroki...
二項分布が正規分布でよく近似されている様子が見えます。 pic.twitter.com/ZmRUVTIPH5
タグ: 統計
posted at 20:21:20
#統計 二項分布Binomial(n, p)におけるピアソンのχ²統計量
(k - np)²/(np) + ((n-k) - n(1-p))²/(n(1-p))
の分布の累積分布函数の視覚化。自由度1のχ²分布で近似されています。
さくっと視覚化で確認できれば不安がなくなり、理解も進みます。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/FswfgwIZZj
タグ: 統計
posted at 20:24:02
@OokuboTact #超算数 とりあえず目を通して欲しい懐かしいサイト
かけ算の順序問題 - Wikipedia
ja.wikipedia.org/wiki/%E3%81%8B...
「掛け算順序固定」問題 - アットウィキw.atwiki.jp/multiplication/
算数の教科書とその指導書の問題点genkuroki.web.fc2.com/sansu/
タグ: 超算数
posted at 20:33:40
OKUMURA, Akira(奥村 曉) @AkiraOkumura
@psr1509 @Hal_Tasaki そうやれば、iPad 側ではツールバーが見えているのに、Zoom の接続先相手には PDF の部分しか見えない状態になる、と。
タグ:
posted at 20:36:24
😷拡散希望😷
西浦教授の応援イラストを描きたい素人絵描きです。
数学よわよわ民なのですが、背景にかっこよさげな数式を書き込みたいのです。現状、適当にベイズ理論のを拾ってきて書き込みましたが意味不明で恥ずかしい。
数学、統計ガチ勢の方、これ書いとけ!というお知恵をお貸しください🙏 pic.twitter.com/sZ3p8Sa6wX
タグ:
posted at 20:38:32
OKUMURA, Akira(奥村 曉) @AkiraOkumura
@Hal_Tasaki @psr1509 これは試してみないとそういう動作になると直感的に想像できないですね。980 円ですが、明日また学生の iPad 借りて試してみて、良さそうだったら自分の iPad にも導入します。ありがとうございます。
タグ:
posted at 20:40:59
さっき初めて計算した.速習本なので仕方ないが条件付き期待値に行く前に注意とか演習で書いてあると嬉しい. twitter.com/genkuroki/stat...
タグ:
posted at 21:04:32
#統計 Kullback-Leibler情報量の定義を
D(q||p) = ∫ q(x) log(q(x)/p(x)) dx
はRadon-Nikodym微分を使って書くと、
D(ν||μ) = ∫ (dν/dμ)(x) log((dν/dμ)(x)) dμ(x)
のように書けます。
しかし、これを知らなくても、例えば渡辺澄夫『ベイズ統計の理論と方法』を読むのに困らない。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 21:23:33
@AkiraOkumura @psr1509 そうですね。やはり試してみるのは大事だと思います。
あと、Good Notes は慣れるとペンやマーカーやポインターを駆使して手書きでバンバン書けるのですが、慣れないうちはペンのつもりで消しゴムで消してしまったりということが頻繁におきます。手書きをするなら修練が必要だと思います。
タグ:
posted at 21:34:56
非公開
タグ:
posted at xx:xx:xx
#超算数 教えている現場の先生がこう言っている。予想した通りである。「引き算の意味」として「残りはいくつ」と教えれば、そりゃ「求差は難しい」となるだろう。
遠山啓や日数教、筑波大学付属小算数部などの、算数教育の権威筋が言う「求差は求残より難しい」はマッチポンプ。 twitter.com/ccc_mihiro/sta...
タグ: 超算数
posted at 22:31:22
@tokyobunkabu 小学校で教える(かけられる数)×(かける数)のルールは、中学で代数計算が始まると明らかに無意味で邪魔なので、多くの人は無意味さに気が付くか、単に忘れてしまうのだと思います。
タグ:
posted at 22:43:03