黒木玄 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 を見て下さい。
2021年07月27日(火)
MU ranakaSei JI | AW @SeijMura
JuliaのMCMCサンプラーか。
単一の言語でできるってのは強みになりそうな気もするなー。
github.com/TuringLang/Adv...
タグ：
posted at 00:02:19
So, bad news. I've been out with some sort of respiratory infection (not COVID, got tested) for most of the past two weeks and I still have no voice, so I'm cancelling my workshop tomorrow.
タグ：
posted at 00:27:28
However, I did open up the Diffractor source code for people who want to take a peek. I didn't really get around to cleaning it all up yet, but it should be usable.
タグ：
posted at 00:27:28
As for the workshop, the plan is to record the lecture parts of it and get it out separately and then have a live session at a later date.
タグ：
posted at 00:27:28
I'm quite disappointed, I had a bit of a spectacle planned, but unfortunately it's just not feasible at this point.
タグ：
posted at 00:27:28
Kareem Carr | Data S @kareem_carr
the first time you work with real data pic.twitter.com/QjoonHNglV
タグ：
posted at 00:48:09
Hi folks - unfortunately the workshop titled "Diffractor: Next-Gen AD for Julia" that was due to take place tomorrow is now cancelled.
The other workshop "Introduction to metaprogramming in Julia" will proceed at the usual time (2pm UTC)
タグ：
posted at 02:24:42
#統計 2x2の分割表
a b
c d
の独立性の帰無仮説を満たすとは限らない確率分布モデルの階層
4つのポアソン分布の直積
↓
↓a+b+c+d=一定
↓
四項分布
↓
↓a+c=一定
↓
2つの二項分布の直積
↓
↓a+b=一定
↓
Fisherの非心超幾何分布
独立性の帰無仮説を満たす場合には最後が超幾何分布になる。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 03:05:09
私は以前はPython matplotlibをJuliaのライブラリとして使って、数学アニを大量に作っていたのですが、現在はそれよりもお手軽で速い #Julia言語 Plots.jlを使っています。
数学アニメ作成にJuliaは非常に優れています。
RもJuliaのライブラリとして使えます(添付画像)。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/norimitsunishi... pic.twitter.com/87LXSS6rvn
タグ： Julia言語
posted at 03:24:35
Buitengebieden @buitengebieden
Some cases may never be solved.. 😅 pic.twitter.com/hLeIMcvfPq
タグ：
posted at 04:10:52
Vector fields, trajectories and ordinary differential equations
bit.ly/3mCDEnP (from Chaos: A Mathematical Adventure)
#math #science #iteachmath #mtbos #visualization #elearning #chaos pic.twitter.com/IxdgBZoWqk
タグ： chaos elearning iteachmath math mtbos science visualization
posted at 07:01:15
KURORO@Pythonレスキュー隊 @KURORO59126227
@9R0iX おはようございます。
> Juliaに書き直したら25倍くらい
おー、素晴らしい！それほど性能差が出るのですね〜
重い処理系 + 簡単なスクリプト類は、juliaへ移行してもいいかもと感じました。
貴重な体験の共有、ありがとうございます。🙇♂️
タグ：
posted at 08:31:18
関数ぶちこんだら，凸関数っぽいかどうか勝手に判定してくれるモジュール，Juliaにないかな？
ヘシアンを自分で計算しないとだめですかね．パラメータ数結構多くてめんどくさい．．
タグ：
posted at 11:58:30
お、これは分かりやすい疑似科学。このグラフ、一見すると「ワクチンを打った直後に死亡が多い（ワクチンで死んでる可能性が高い）」ように見えますが、実際は「接種から日数が経ってるほど『ワクチンは無関係』と判定されて報告されなくなる」というグラフですね。
twitter.com/NatsuPerson/st...
タグ：
posted at 12:14:18
モンテカルロ法に関する集中講義用に作った講義ノートのうちコラムだけ集めたPDFです。モンテカルロ法に限らず物理系の数値計算について私が考えていることをまとめてみました。興味がある方はぜひ読んでみてください。コメント等も頂けると幸いです。
www.dropbox.com/s/a26c00i0n0xq...
タグ：
posted at 12:35:09
@otapediatrician 一応、元ツイートの人は「単にワクチン接種後に死んだのではなく、ワクチンのせいで死んだのだ」って主張するためにこのグラフを作ったようですが、「報告するかどうかの時点で既に選別されている」という点を全く考えてないせいで、勝手に数字に騙されてるんですよね。
タグ：
posted at 12:39:50
#統計 以上の4つのモデルに関する独立性のカイ二乗検定を出す計算(最初の3つはWilksの定理を使い、最後の超幾何分布の場合はスターリングの公式を使う)の詳しい説明は
nbviewer.jupyter.org/gist/genkuroki...
r×cの分割表におけるPearsonのχ²統計量
にある。Wilksの定理を使うには最尤法の計算が必要になる。
タグ： 統計
posted at 13:32:54
EARLの医学ツイート（一時避難所） @EARL_med_tw2
これ訴えたら勝てるんじゃないかしら twitter.com/nasuyuko/statu...
タグ：
posted at 13:48:37
#Julia言語
「自動微分でHessianを求めるようにして、Hessianが正値かどうかを数値的に確認する」という方法である程度確認できると思います。
ケースバイケースで試行錯誤が必要な面倒な問題だと思います。
自動微分のパッケージは色々ある。
最も素朴なやつが
↓
github.com/JuliaDiff/Forw... twitter.com/physics303/sta...
タグ： Julia言語
posted at 13:52:59
#Julia言語 #自動微分 最近公開された開発中の新世代の自動微分パッケージの
github.com/JuliaDiff/Diff...
のソースコードを見ると、TangentBundle, TaylorBundle, JetBundle のような幾何用語の型名がガンガン出て来ます。
そういう幾何的概念とコンパイラ技術の橋渡しをやっているとみなされる。
posted at 14:02:36
#Julia言語 Julia discourseを見ていても、Pythonからの移行組がPythonへの過学習が原因で変な癖がついてしまっていて、Juliaで困難を覚えている場合が非常に多いという印象があります。
最初に
docs.julialang.org/en/v1/manual/n...
と
docs.julialang.org/en/v1/manual/p...
を見ておくと苦労が大幅に減ると思う。
タグ： Julia言語
posted at 14:16:44
五輪中、メディアのリソースも情報を受ける側の関心のリソースも薄くなりますね。でも今起きているのは今までにない形（40-50代を中心にした）での逼迫の入り口なので、注意を喚起し続けるべしです──新型コロナ: ［社説］政府は感染急拡大への危機意識高めよ: 日本経済新聞 www.nikkei.com/article/DGXZQO...
タグ：
posted at 14:50:14
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#数楽 「常微分方程式を解くこと」は「矢線に沿って矢線の長さの速さで進むこと」だと知っていることはとても大事。
「手計算で解ける易しい微分方程式達」と上の段落の知識の組み合わせがあれば、その後やれることが増える。 twitter.com/74wtungsteno/s...
タグ： 数楽
posted at 15:18:24
ほんと、五輪前ならトップニュース扱いのものが、紙面でもポータルサイトのトップにも出てこないくらいになっていると思うんです。「市民が危機感を持つ」ことは感染制御の大きな手立てなので、現状はやはりまずいです。
タグ：
posted at 15:23:14
現在では #Julia言語 でもBLAS-likeライブラリを書けそうなことも実証済み。
github.com/JuliaLinearAlg...
すでに3年前に時点で、気候モデルのようなスーパーコンピューティングの分野でFortranでモデルを書いていたときよりJuliaで3倍速く計算できていた(既出)。
medium.com/codex/is-julia...
タグ： Julia言語
posted at 15:46:02
#Julia言語 の
docs.julialang.org/en/v1/manual/p...
の最初の項目は「グローバル変数を避けよ」です。より正確には「グローバル変数達を函数の中で使うときにはすべてを引数として渡せ」。それが面倒だと思う人はスタイルに問題があるから。以下のリンク先のスタイルに従えば解決。
github.com/genkuroki/publ...
タグ： Julia言語
posted at 17:31:37
五輪実施に賛成した人達は特に「無事に済むこと」を心から願っていると思います。
リスクを下げるためには、感染制御のための危機感を喚起する努力を率先して行った方がよいです。
五輪実施に賛成した人達が実際にそれをやると非難されるでしょうが、その手の非難を気にしてはいけないと思います。 twitter.com/rsider/status/...
タグ：
posted at 17:40:14
#JuliaLang
Solving differential equations in parallel on GPUs | Workshop | 2021
youtu.be/DvlM0w6lYEY
↑
github.com/luraess/parall...
↑
discourse.julialang.org/t/scaling-for-...
タグ： JuliaLang
posted at 18:08:17
Togetter（トゥギャッター） @togetter_jp
登山の前では日焼け止めは大して役に立たないと一目でわかる写真がこちらです「ひいっ…」「ゲンスルー戦のゴンだ」
togetter.com/li/1749873.
タグ：
posted at 19:30:07
算数用語の「式の意味」は常識的な感覚ではその意味を正しくとらえることができない特殊な言葉。非常識な世界で常識は通用しない。要注目。 #超算数 twitter.com/ookubotact/sta...
タグ： 超算数
posted at 19:31:04
#統計 再度強調しておきますが、95%信頼区間の95%について「確率か、割合か？」という問いを発すること自体ひどくミスリーディングです。「頻度論云々」とか言い出すのもおそろしく低レベル。
その95%はモデル内標本分布で測った確率=割合にすぎません。
モデル内標本分布が見えていることが大事。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 19:50:07
私の経験では、「頻度論」という用語を使っているだけで、伝統的に言われているおかしな考え方を受け入れることによって、シンプルに確実に正しい考え方をできなくなっている人になっていることが多い。当該の議論で必須の概念かどうかよく考えるべき。
ちょー算数における「内包量」に似ている。
タグ：
posted at 19:50:08
自分から見て低レベルではないと評価されているように見える立派な人達が使っている用語を使っているから、自分も低レベルではない、というような判断の仕方は論理的ではありません。
もっとシンプルに「何なら確実に正しいと言えるか」を、「95%」の出所を定義に戻って確認することが大事。
タグ：
posted at 19:53:53
ポイント：シンプルに「何なら確実に正しいと言えるか」を定義に戻って確認することが大事。
これ、できるようになるまで、めっちゃ努力しても数年以上かかる。すぐにできることではないので、その基本的重要さを正しく認識すること自体、相当に難しい。
タグ：
posted at 20:00:20
多くの人は、有名教科書や有名査読論文や有名同業者達を参考にし、シンプルに定義に戻って確認することをしない。
だから、非常に優れた人であっても、数学的にややこしい事柄については、論理と証拠ではなく、伝統と権威に従うという明瞭に非科学的な誤りを犯しがちになる。
タグ：
posted at 20:03:34
定義に戻って確認すれば、95%信頼区間の95%は信頼区間を定義するために使ったモデル内部での標本分布で測った確率になっていることがわかる。
確率の値は確率測度を決めないと決まらない。確率測度の意味での確率は割合と同じ意味。
モデル内での確率=割合の現実での解釈は別の問題。
タグ：
posted at 20:11:04
95%信頼区間の95%が信頼区間を定義するために使ったモデル内部での標本分布で測った確率(もしくはその近似値)になっていることは、信頼区間の定義の仕方から純粋に論理的に出て来る絶対的な真理です。
こういう絶対的基盤があるのに、それを無視して、伝統的な用語法でふわふわした議論をしちゃダメ。
タグ：
posted at 20:15:25
#数楽 確率測度の意味での確率の定義(=割合概念の抽象化)は、ランダム性がなくても使えるので、応用先はものすごく広くなります。
決定論的な手続きで生成される擬似乱数列にも適用できる(笑)
実際、私は #Julia言語 を使ってそういう数値実験のデモを数え切れないくらい見せているし。 twitter.com/toyo9/status/1...
posted at 20:27:59
私のパソコンの師匠がやねさん(第29回世界コンピュータ将棋選手権優勝ソフトやねうら王の開発者)でとても感謝している。
人間性とか人格とか色々認めたくない部分はあれどプログラミングに関してあれほど純粋かつ真面目で情報に精通している人がよく教えてくれるな🤔
今度焼肉(安いの)おごったるよ！
タグ：
posted at 20:41:59
因みに教え方
私「やねさん、この動画ソフトなんだけど」
やね「ふむふむ・・・チュートリアル見たらええやん」
私「このサイトでこれ買うにはどうしたら」
やね「ググれば買い方書いてあるやん」
です。
懇切丁寧に教えてくれる、マジでいい師匠です🙇✨
タグ：
posted at 20:44:14
非常に困ったことに、フィクションに過ぎない正規分布モデル内での確率=割合に過ぎないものについて、正規母集団の仮定を正当化抜きに行って議論を先に進める入門的教科書が多い。
実際、多くの人が「学生時代に、正規分布を仮定してよい理由がわからなかった」と言う。
この点を看過しちゃダメ。
タグ：
posted at 20:52:50
現実には、母集団分布が正規分布になっている保証がないどころか確実に正規分布からずれていることがわかっている場合にも、正規分布モデルに基く信頼区間が使われたりする。
そういう実態を隠して「正規分布を仮定する」と言うのは不誠実だと思います。
タグ：
posted at 20:52:50
単なるフィクションでしかない数学的モデル(例えば正規分布モデル)内での確率に過ぎないものについて、「正規分布を仮定する」のような言い方でモデルと現実の区別が曖昧になるような前提を置いて教わると、シンプルでクリアに考えることができなくなるのは当然だと思います。
タグ：
posted at 20:52:51
かけ算順序指導に強く影響を受けてしまった人達が大人になってからもシンプルな事柄について訳の分からないことを繰り返し述べているのと同じように、高等教育で不適切な教育を受けてしまった人達が95%信頼区間の95%についてもシンプルにクリアに考えることができなくなっているのだと思います。
タグ：
posted at 20:55:41
高等教育の内容が不適切であっても、自分でモデルをいじって、自分で信頼区間を定義してコンピュータに計算させて、同じデータに基く信頼区間がモデルの変更によってどのように変わるかを調べたことがあれば、95%信頼区間の95%について誤解しようがないです。
しかし、そこまでやっている人は少数派。
タグ：
posted at 21:00:31
95%信頼区間の95%はモデル内確率なので、モデルを変えれば、信頼区間も当然変わります。
この点が曖昧なままだと、重要な政策に関わる区間推定の結果の評価に失敗する人達が高等教育を受けた人達の中からも続出することになり、社会的に大変まずいことになる危険が増えます。
タグ：
posted at 21:04:00
シンプルに定義に戻ってクリアに考えれば、95%信頼区間の95%がモデル内確率に過ぎないことにすぐに気付き、モデルを変えれば信頼区間も変わることも単なる常識的な考え方に過ぎなくなります。
「確率か、割合か」とか「頻度論云々」はくだらない議論に誘導する社会的に有害な言説だと私は思います。
タグ：
posted at 21:07:44
@benatawaomi 再びのリプです😅兄弟子筋のたややんさん(世界コンピュータ将棋オンライン大会優勝ソフト水匠の開発者)とても優しくて
私「○○より△△の方が使いやすいんですか？使い方難しくて・・・」
たややんさん「○○でいいです。使い方はチュートリアル観たり誰かの解説動画観ればわかると思います😄」
タグ：
posted at 21:26:11
本日、日本将棋連盟様との間で、関西将棋会館の高槻移転についての合意書を締結いたしました。
明日より、会館建設資金について、ふるさと納税制度を活用した寄付受付を開始したします。ご支援よろしくお願い申し上げます！camp-fire.jp/projects/view/...
タグ：
posted at 22:26:52
OokuboTact 大久保中二病中年 @OokuboTact
小学校の算数では「式」へのこだわりが異常だったりする。
これは最近の傾向ではなく、戦前からあったりする。
#超算数 のタグで天むす名古屋さんが資料を公開しているので、興味のある人はツイートを探してください。
しかし現代の「式」のごだわり教育が戦前のままかと言えば、そうではなく
タグ： 超算数
posted at 22:39:47
Norimitsu Nishida @NorimitsuNishi1
@genkuroki @genkuroki さん
情報ありがとうございます！Jupyter Notebookの内容をみて、JuliaでRを簡単に呼び出せるのすごいと思いました。ただ、私の環境でエラーがおこってRCallが実行できず、ちょっとJuliaで初グラフアニメーションを試みてみようと思います。
タグ：
posted at 23:22:49
OokuboTact 大久保中二病中年 @OokuboTact
#超算数
しかし現代の算数教育の「式へのこだわり」は、実は小学校だけではなく、中学校の数学教育とも大きく関係している。
最近の「式の意味」については「新版 算数科教育研究」＜式の表現と意味＞ （文：清野辰彦）が参考になった。
タグ： 超算数
posted at 23:29:15
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
この本を教科書として使っている大学の講義もある
pic.twitter.com/xmHSsYBITK
タグ： 超算数
posted at 23:34:04
