黒木玄 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年06月27日(土)
#Julia言語
Juliaの文法自体が多くのプログラミング言語へのオマージュになっている感じ。
変数 a が "" や `` の間で $ a ($ の後の空白は除く)で a の中身が展開されることは、shell へのオマージュ。
最近 shell script の代わりに Julia script を使い始めた。
twitter.com/aki202/status/...
タグ： Julia言語
posted at 00:27:32
サンパウロのとあるビルに毎日朝夕200羽以上のメジロメキシコインコが集まってはレンガをかじっていくそうな。一般にインコの仲間はミネラル摂取や毒素排出のために土を食べるけど、これまた都市適応・意図せざる餌付けの特殊例やね verdesp.com.br/o-predio-que-a...
タグ：
posted at 00:31:10
ごまふあざらし(GomahuAzaras @MathSorcerer
Julia コミュニティ茶目っ気成分多め（真面目にするとすごい人）な人が多い気がする． twitter.com/genkuroki/stat...
タグ：
posted at 00:31:38
#Julia言語 つい最近更新された2013年のブログ記事
julialang.org/blog/2013/04/p...
では Julia v0.1 の書き方で、shell script のようなことを Julia でやる方法が解説されていた。更新は現在の Julia v1 で動くようなコードの変更。
Juilia開発者達は開発初期から shell にも敬意を払っていた！
タグ： Julia言語
posted at 00:31:48
Ligun (りぐん) ☠️ Nenio @choonggeun
Juliaやろうと思って早一年が過ぎた
そろそろ手を出す時か？ twitter.com/aki202/status/...
タグ：
posted at 00:41:04
@mtkharu3 @kennel_org #超算数 「チョー算数」指導は、現場の教師が自分自身の経験と考察で行なっている教え方ではなくて、算数教育界でのカースト上位の人たちが算数教科書の教師用指導書(一般人購入不可)などで広めている教え方なんです。
「一般教師よりもずっと悪い人達がいる」という見方の方が正しいと思います。 twitter.com/kennel_org/sta...
タグ： 超算数
posted at 00:56:17
ごまふあざらし(GomahuAzaras @MathSorcerer
alan-turing-institute 結構Julia で数理的な教育をしていこうという意志を感じていてすごくいい．
タグ：
posted at 01:32:55
"What's Bad About #JuliaLang"
Presented by Jeff Bezanson
Livestream: www.youtube.com/watch?v=TPuJsg...
Abstract: pretalx.com/juliacon2019/t...
#JuliaCon #JuliaCon2019
タグ： JuliaCon JuliaCon2019 JuliaLang
posted at 02:20:00
Hideki Kawahara: WAS @hidekikawahara
これだけニッチな仕事だと、論文を書いてそれなりに使えるツールを公開するところまで、自分で仕上げなければならない。FVNは本当に面白くてお薦めです。MATLABに加えてJuliaでも実装しておくと、面白がる人が増えるだろうか？
#MATLAB #Julialang twitter.com/hidekikawahara...
posted at 03:27:43
@neuralengine @MatthiasHennig6 @NPRougier @neuralreckoning @khinsen @gnode I'm trying to move from R to Julia; Julia is good friends with R and python, so I'm hoping this will be a good result.
タグ：
posted at 04:23:14
I am currently writing a sort of parser generator in #JuliaLang and have developed the term "Language Definition Language Language Definition."
The definition of a language used to define another language. 🙃
タグ： JuliaLang
posted at 05:53:51
財政出動をバカにするリフレ派は居ない。居たとしたらそれはリフレ派ではない。リフレ派の言い分は「金融緩和をせずに財政出動だけをしてもダメ」と言う主張。実際にこれをやったのが麻生政権。結果は極度の円高になり物価は上がるどころかむしろ下がった。つまりデフレを加速させたんだよ。 twitter.com/habana1000/sta...
タグ：
posted at 06:05:53
非公開
タグ：
posted at xx:xx:xx
@mtkharu3 @kennel_org #超算数 チョー算数の目的を【児童に対する精神的支配】と勝手に想像する自由もありますが、そういうことを公言する場合には証拠に当たる情報があることを必ず述べるべきだと思います。
それができないなら、証拠がないのにひどいことを言っただけになってしまいます。
twitter.com/kennel_org/sta...
タグ： 超算数
posted at 08:59:13
非公開
タグ：
posted at xx:xx:xx
@mtkharu3 @kennel_org #超算数 掛け算順序固定強制指導は100年以上前からあります。
各種の資料を見ると、掛け算順序を逆にすると誤りになると本気で信じているように見えます。
そして「児童に対する精神的支配」が目的であることを示唆する資料を私は見たことがありません。
twitter.com/genkuroki/stat...
タグ： 超算数
posted at 09:05:03
@mtkharu3 @kennel_org #超算数 あと、チョー算数問題について「教師バカ論」には与するつもりがないとはっきりと再度述べておきます。
平均以上の知性の持ち主でも(チョー)算数教育を真面目に勉強するとおかしくなるという問題。
たとえば、オウム真理教の理系高学歴者達も知能が低いという意味での馬鹿ではなかった。
タグ： 超算数
posted at 09:12:19
@mtkharu3 @kennel_org #超算数 実際、掛け算順序固定強制指導を氷山の一角とするチョー算数問題への理解は、小学校の先生からの情報提供(特に現場の様子や教科書の教師用指導書の情報)のおかげで深まって来ました。
掛け算順序でバツをつけたくない先生も周囲に圧力でそうせざるを得なくなっている場合もある。
タグ： 超算数
posted at 09:15:32
@mtkharu3 @kennel_org #超算数 ツイッターなんかでチョー算数擁護で目立っている先生個人を厳しく批判するのは構わないと思うのですが、それを安易に一般論に拡張するのはひどい誤りであり、そういうひどい誤りをおかしたら、誤りを認めて反省するべきだと思います。
タグ： 超算数
posted at 09:18:28
@mtkharu3 @kennel_org #超算数 例えば、漢字のトメハネがちょっと違っているだけで誤りとみなすことは漢字文化的に誤りであることを政府機関(文化庁)はきちんと広報しています。専門家達の助言が効いている。
それに対して、チョー算数は主に専門家達が宣伝普及につとめています。ここにチョー算数特有の問題がある。
タグ： 超算数
posted at 09:25:27
#Julia言語 MITのJuliaを使った講義
github.com/mitmath/18S096...
では最初の方で「メモリ上での配列の実現の仕方を意識しないといけない理由」が説明されています。
github.com/mitmath/18S096...
nbviewer.jupyter.org/github/mitmath... twitter.com/aioooououou/st...
タグ： Julia言語
posted at 09:42:26
新型コロナ専門家会議廃止に関する新情報で驚き。この動画は見たほうがいいです。与野党ともに聞いていないと。公明党は持ち帰ると。さらに、実はまだ廃止はしていないと内閣府。何がなんだか。不安。
専門家会議廃止 西村担当相の表明について 2020.6.26 youtu.be/il4Xddey8ek @YouTubeより
タグ：
posted at 10:10:20
Progress is slow because of the upcoming job market, but I finally got this code working in #julialang, so it can be called from R/Python! Feeling good, will post soon I hope. Hoping to make using @GioCompiani's estimator (and my in-progress model selection work) easy twitter.com/jamesbrandecon...
タグ： julialang
posted at 10:32:23
#Julia言語 版
nbviewer.jupyter.org/gist/genkuroki...
JuliaでSymPyの式exprを函数化するには
juliaFunc = eval(Meta.parse("() -> $expr"))
が速いです。
Juliaでも普通にSymPyを使えることにも注目。Juliaは多くのプログラミング言語へのオマージュにもなっている高速強力な糊言語。糊言語も速い方がよい。 twitter.com/GenmVisionary/... pic.twitter.com/2mH7F5tzuY
タグ： Julia言語
posted at 10:56:44
#Julia言語 SymPyの式 symbolFunc = 2x+y について
”(x, y) -> $symbolFunc"
が
"(x, y) -> 2*x + y"
に変換されていることにも注目。
こういうことをするときには
気分は shell script
なのですが、高速な数値計算も可能になる！
様々な道具を貼り合わせての高速計算にJuliaは向いている。 pic.twitter.com/7SOlbZcJw9
タグ： Julia言語
posted at 11:03:31
#Julia言語 特殊函数の計算で連分数展開を用いるコードを書く場合があるのですが、分かり難くてバグが発生し易いです。
Juliaにおける数値計算用のコードにSymPyのシンボルを引数として与えると数式を表示してくれます！
式を見れば正しい計算をしているかが一目瞭然！
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 11:06:58
#Julia言語 2つ前のツイートの
気分はshell script
については以下のリンク先スレッドも参照。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 11:08:37
#Julia言語 で指数積分函数のコードを書くという宿題がMITの講義で出されていたのですが、その模範解答はものすごくて、Juliaを使ったアルゴリズムの最適化によって、scipyで採用されているFortranライブラリの5～6倍の速さを持つ函数とpure Juliaで書くことに成功！
nbviewer.jupyter.org/github/steveng... pic.twitter.com/3E1mErAY05
タグ： Julia言語
posted at 11:13:10
#Julia言語 1つ前のツイートで紹介した
nbviewer.jupyter.org/github/steveng...
は特殊函数の数値計算のコードの最適化の良い手本になっています。なぜならば、最適化の過程は原理的にJulia以外にも通用するからです。
続く
タグ： Julia言語
posted at 11:15:20
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
「物理の問題を機械学習を使って解いてみました」，にはあまり関心がないんだけど，「物理と機械学習で両方に定義できる量について考察しました」「物理の手法を機械学習に流入しました」，みたいな話にはものすごく興味がある．エントロピーとか，何かの保存則とか，平均場近似とか．
タグ：
posted at 12:26:35
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
ごまふあざらし(GomahuAzaras @MathSorcerer
julia -O0 (オーの後にゼロ）をつけてREPL起動して
```
using Plots
plot(something)
```
するとちょっとだけJIT待ちが減るね．
タグ：
posted at 13:59:44
なるほどです！
バーゼル問題が解けました！
twitter.com/genkuroki/stat...
#toyoの数学メモ pic.twitter.com/DgTBRcaoY0
タグ： toyoの数学メモ
posted at 15:36:49
∫(xsin(x)/(1+cos(x)^2))dx
ものすごく複雑に見えます！！！
www.wolframalpha.com/input/?i=%5Cin...
#toyoの数学メモ twitter.com/genkuroki/stat... pic.twitter.com/CbSnzmARtK
タグ： toyoの数学メモ
posted at 15:58:56
Perfectionatic @perfectionatic
@randal_olson Mine in #julialang worked much better.
twitter.com/perfectionatic...
タグ： julialang
posted at 16:43:45
#数楽 ものすごく複雑で人間が手で扱うのはしんどいですが、その式の本質的部分はダイログLi₂(z)が出て来ているところだと思います。それ以外の部分は初等函数と初等函数の逆函数で書けている。
個人的な意見では単に面白いという理由でLi₂(z)も微積分での基本的な例としたい。😊 twitter.com/toyo9/status/1...
タグ： 数楽
posted at 16:47:26
JuliaのドキュメントのTypesのページがあまりにも長くて自分の英語力だと読んでる最中に忘れてきそうだから書き込みをしようとPDFにしてiPadに読み込んだのだけどいざ読み込むと長すぎて笑っちゃった pic.twitter.com/6GN87ZEzWG
タグ：
posted at 16:49:05
#数楽 Li₂(z)などなどを「教育用の基本例」に含めたいという意図で書かれた微積分のノート群が
github.com/genkuroki/Calc...
にあります。
一般論よりも具体例の方が重要という主張もこめたつもり。
ダイログの五項関係式の証明は
nbviewer.jupyter.org/github/genkuro...
にあります。
タグ： 数楽
posted at 16:56:38
ごまふあざらし(GomahuAzaras @MathSorcerer
Snappy のツイートを見てなんか動いたのはわかった。
（これが何を意味するのかがわからない。）
```console
$ pip install snappy
$ python
>>> import snappy
>>> M = snappy.ManifoldHP("14n12345")
>>> M.plink()
``` pic.twitter.com/FI4nNMR2pI
タグ：
posted at 17:00:36
#数楽 続き～、Lerchの超越函数(Hurwitzゼータとポリログの両方の一般化)まで拡張しておきたい(函数等式はLerchの超越函数だけで閉じる)。
対数ガンマ函数がHurwitzゼータから出てくることも説明しておきたい(Stirlingの公式より深い)。
などなどが
github.com/genkuroki/Calc...
のノートでは実現済。
タグ： 数楽
posted at 17:04:09
#数楽 #Python
結び目の図を見るだけでなぜか楽しい。
これは面白そう。
github.com/3-manifolds/Sn...
【SnapPy is a package for studying the topology and geometry of 3-manifolds, with a focus on hyperbolic structures. It is based on the SnapPea kernel written by Jeff Weeks.】 twitter.com/mathsorcerer/s...
posted at 17:14:07
非公開
タグ：
posted at xx:xx:xx
#統計 「不偏分散の定義でn-1で割る理由」とか、「正規分布のサンプルX_1,…,X_nについて、その標本平均をMと書くとき、
(1/σ²)((X_1 - M)²+…+(X_n - M)²)
が自由度n-1のχ²分布に従うこと」は本質的に内積の線形代数です。雑な解説が添付画像にあります。
線形代数知らないと無駄に面倒になる。 pic.twitter.com/Na1nNGF4ru
タグ： 統計
posted at 17:24:47
#Julia言語 #JuliaLang
Why Julia - Package Development in Under 4 Minutes youtu.be/vVywlAgyedI
posted at 17:36:12
Hiromitsu Takagi @HiromitsuTakagi
いや、ある意味間違っていませんよ。（国がとらえた）ある人に接触したそれぞれの人に対して警告メッセージを（国が）発出できるわけなので。だからこそ目的外利用が厳格に禁止されなければならないし、目的外利用しないことを国が保証しなければならない。
twitter.com/shinout/status...
タグ：
posted at 17:40:10
非公開
タグ：
posted at xx:xx:xx
#Julia言語
Juliaのconstは「型変更が不可能なおかげで型推論を邪魔しないで済むグローバル変数で値を変更すると警告が出るもの」という感じ。 twitter.com/ultimatile/sta...
タグ： Julia言語
posted at 17:45:39
Hiromitsu Takagi @HiromitsuTakagi
このことは、この番組内で述べていたところ（15:38から）。録画放送が明日まで。
twitter.com/suzukimasatomo...
タグ：
posted at 17:45:56
遅くなりましたが、これの講義資料(本編27ページ+付録7ページ)です。(量子鍵配送は未収録)
researchmap.jp/multidatabases...
量子情報といいつつ、メインは量子論部分です。素直に量子情報をやれと言われそうですが、これがないと始まりませんし、量子論特有の話をやるには結局基礎が重要になります。
タグ：
posted at 17:48:18
量子力学の授業でも内部的に量子論の話をすることになるかと思いますが、大抵の場合、量子論の部分が結構おざなりだと感じます。それは結局のところ、現状で実験可能な散乱問題の取扱いをゴールとするならそれで十分だったからだと理解しています。ただ、今後は変わっていくでしょう。
タグ：
posted at 17:48:18
通常の量子力学で採用している公理系は von Neumann の量子力学の数学的基礎がベースになっていると理解していますが、あまりそこら辺を明確にしていないので、直感的な理解云々の前にそもそも量子力学がどういう理論なのか(例えば現実とどう対応するか)すらよくわからなくなるのではないでしょうか。
タグ：
posted at 17:48:18
一方で、なまじ von Neumann の公理系でやろうとすると射影仮説が気持ち悪いので、無駄にそこでハマりそうになります。
そうなので、本講義では一般確率論の特殊例としての公理系を採用して量子論を導入しています。これは、より操作主義的な公理系なので(個人的には)より受け入れやすいものです。
タグ：
posted at 17:48:19
もっとぶっちゃけた理由としては、von Neumann の本をちゃんと読んだことがなくて、よく分かってないというのもあります。(一般確率論もよくわかっていないとかいい始めると悲しくなってきますが...)
タグ：
posted at 17:48:19
観測可能な物理量や実現可能な測定の制限とかをシリアスに考えなければ、これら2つの公理系は間接測定を通していったり来たりできるはずなので、どちらかに優越性があるわけではないですが、一般確率論的に行う(個人的な)メリットとしては、物理量という概念が公理に出てこないということがあります。
タグ：
posted at 17:48:19
個人的な感想でいえば、物理量というのは極めて古典的な概念に思えるので、そのようなものを排除した公理系で(少なくとも)枠組みは記述できることを見せることは、理解のためにもいいのではないかと思います。
タグ：
posted at 17:48:20
あと、この講義の中身は通常の量子情報ネタに芝浦工大の木村先生とKEKの筒井先生から教わった内容を時間コスパ重視で適宜加えたものになっています。
そうなので量子情報科学入門の木村先生が書いている部分と似通っているのは必然です...。筒井先生の方はいつになったら教科書を出すのか...。
タグ：
posted at 17:48:20
話が長いですが、量子論の話を授業レベルに落とし込むというのはあっても数が少なかったと思うので、今回の資料がそのような試みの参考になれば幸いです。
ただ私の場合は有限次元限定です。そうじゃない場合は堀田さんが教科書を書いていると言っていますので、それ待ちですかね。
タグ：
posted at 17:48:20
#Julia言語
Googleの設定を
地域：アメリカ合衆国
にすれば、Juliaの検索をキーワードJuliaで便利に行えるようになります。(Julialangをキーワードに使うよりずっと優れた方法)
プログラミング関連の検索が多いなら「地域：アメリカ合衆国」としておく方が合理的かも。
twitter.com/katsuhitosudoh...
タグ： Julia言語
posted at 17:50:41
Hiromitsu Takagi @HiromitsuTakagi
高木「本当に目的外の利用ができないのか、悪用ができないのかというと、たしかに、手元にあるデータだけ見て何かをしようとしてできるかというと、できそうもない気がするでしょう。でも、システム全体を見れば、悪用の余地はある。例えば香港みたいな国や地域の状況、つまり政府に対する抗議活動が…
タグ：
posted at 18:07:13
うっかりリンクを踏んだら画面からはみ出しちゃいましたよ！😱ヒェｪ twitter.com/genkuroki/stat... pic.twitter.com/5IrGx0jAmB
タグ：
posted at 18:08:29
Hiromitsu Takagi @HiromitsuTakagi
…抗議活動が盛んな状態に国がなっているとき、活動を先導しているリーダーが捕まったり、政府がその人のスマホを手にしてそこにこのアプリが入っている場合、そのリーダーと2週間以内に15分以上接触した人に対してメッセージを送れるわけですよね。そのメッセージの内容がカスタマイズできるなら、…
タグ：
posted at 18:09:42
Hiromitsu Takagi @HiromitsuTakagi
…『お前ら、分かっているぞ。あの抗議活動に参加しただろう』と忠告を出すことができる。誰であるか『特定』できてなくても、そういうメッセージを政府が出すことができる。そういうシステムですからね。メッセージが出るだけでも、行動の自由を委縮させて人々の自由がなくなるという効果は発生する…
タグ：
posted at 18:14:00
Hiromitsu Takagi @HiromitsuTakagi
…発生するので、悪用の余地がないとはいえない。これが公衆衛生のためであれば、「あなたは接触しています」と出ることは本人も望んで入れているわけですから問題とならないわけであり、目的外利用の余地があることはそこは隠さずに言っていかないといけないと思います。」
タグ：
posted at 18:15:28
Hiromitsu Takagi @HiromitsuTakagi
鈴木「面白い使い方ですよね。牽制効果は十分に発揮できそう。捕捉されてると思ってビビる。デモへの参加を委縮させるコントロールには確かに使える。そこまで手間はかけないとの信頼はあるでしょうけどね。ただ、リスクは全部洗いだして開示とかして逆に信頼感を出していくのは政府の仕事ですよね。」
タグ：
posted at 18:21:24
ランダムウォーク失敗中 #fractal #Julia言語 pic.twitter.com/31tMeDjUnH
posted at 19:09:16
#Julia言語 私も遊んでみた。
添付画像1: StackOverflowError
添付画像2: ifの分岐でがんばった改良版
添付画像3: multiple dispatch でシンプルに書ける
他人が書いたif文を読むのは結構難しいと思う。
nbviewer.jupyter.org/gist/genkuroki...
twitter.com/abap34/status/... pic.twitter.com/pqtv1NUtvB
タグ： Julia言語
posted at 19:22:52
山中さんがメンバーだと何がいいのか、全然わからない。感染症・感染症対策の専門家がいなくちゃ話にならないでしょうに。富岳だろうがなんだろうが、専門家抜きには意味がないと思いますけどねえ twitter.com/nishy03/status...
タグ：
posted at 19:47:29
をー、Perlまだ現役なのかー
約25年ぶりのメジャーバージョンアップ「Perl 7」が発表される、Perl 5への下位互換性は維持 - GIGAZINE gigazine.net/news/20200625-...
タグ：
posted at 19:48:50
@mtkharu3 @kennel_org #超算数 添付画像のようなクズは結構見かけますよね。
学校の先生であるか否かと無関係に結構見る。
そういう個人を名指しでクズ扱いして行く社会であってほしいと思います。
この手の人物は教育に関する議論に関わるべきではない。
twitter.com/mtkharu3/statu... pic.twitter.com/GhZr18swWq
タグ： 超算数
posted at 19:54:57
@mtkharu3 @kennel_org スレッドが分岐しすぎると後で面倒なので以下のリンク先でコメントしておきました。
twitter.com/genkuroki/stat... twitter.com/genkuroki/stat...
タグ：
posted at 19:56:43
@genkuroki @mtkharu3 @kennel_org #超算数 算数教育には【[式は]自分の考えを間違なく伝えるために、重要なものであると考える】傾向が70年以上前から散見します(和田義信 1950)twitter.com/temmusu_n/stat...。子供への精神的支配は掛順の目的」ではなくても、そのような「利用」を許してしまう土壌は、公開資料に現れていると思います。
タグ： 超算数
posted at 20:45:35
#統計 なぜか清水裕士さんにブロックされているのでこちらにコメント。
docs.google.com/presentation/d...
の件について。
添付画像中のナマの階層モデルはサンプルサイズ1の場合になるのでWBICはwell-definedではなくなります。
だから、周辺化せずにWBICを計算するのは無理です。続く twitter.com/ohkubo_yusaku/... pic.twitter.com/RSl6XglBu9
タグ： 統計
posted at 21:24:32
非公開
タグ：
posted at xx:xx:xx
#統計 β = 1/log 1 = ∞ だと破綻すると書いたのですが、β→∞ の絶対零度への極限は、最尤法への極限なので、その場合の「ほぼ定義通りの事後分布」は最尤法の解で(λ^, θ^)のデルタ分布とみなすべきでしたね。続く
twitter.com/genkuroki/stat...
タグ： 統計
posted at 22:00:50
#統計 結論を訂正しました。
(サンプルサイズ1のWBIC)
= (サンプルサイズ1のBIC)
= - log p(Y|λ^, θ^)
ここで(λ^, θ^)はナマの階層モデルの最尤法の解。
これでどれだけ正しくモデル選択できるかの数値実験ならばそう難しくないと思います。
twitter.com/ohkubo_yusaku/...
タグ： 統計
posted at 22:10:21
オープンソースで、ライセンスは緩く、Rubyのように動的で、Cのように速く、Pythonのように汎用的で、Matlabのように数式が描け、Perlのように文字列を操作できる言語、がJuliaだって。本当ならみんな使うだろう
タグ：
posted at 23:28:21
#統計 論文
arxiv.org/abs/1908.10572
ではWBICには少し小さめになるバイアスがあることが指摘されているのですが(バイアスの推定量も作っている)、私の #Julia言語 を使ったノート
nbviewer.jupyter.org/gist/genkuroki...
にあるBernoulli分布モデルでの公式でもそうなっています。 pic.twitter.com/s0JylIILXw
posted at 23:34:32
数式描画や３Dプロットを追加して現在２３個です。
#Julia言語 #julialang
koisignal.com/?p=422
posted at 23:57:49
