黒木玄 Gen Kuroki(@genkuroki)/2020年03月27日

ごまふあざらし(GomahuAzaras @MathSorcerer

Qiitaの記事、全部githubに移行しようかな

タグ：

posted at 00:02:59

Franklin.jl twitter.com/hika_ruriruri/...

タグ：

posted at 00:04:04

Arturo Erdély @ArturoErdely

ごまふあざらし(GomahuAzaras @MathSorcerer

Improved thread management in #julialang for background and nonuniform tasks: github.com/tro3/ThreadPoo...

タグ： julialang

posted at 00:17:00

@hika_ruriruri 今のところC,Python,Juliaの実行結果を表示する機能は確認できています．自分のHPのコードは

github.com/terasakisatosh...

です．Franklin.jl のドキュメントも充実しています．

タグ：

posted at 00:48:47

Yahoo!ニュース @YahooNewsTopics

【米失業率世界恐慌匹敵の恐れ】
yahoo.jp/ob9MdV

新型コロナ拡大に伴う休業などの影響で、米国では新たに失業保険を申請した件数が21日までの1週間で328万3千件に上った。失業率は4月以降、30年代の世界恐慌時の水準まで急速に悪化する可能性が出てきた。

タグ：

posted at 01:08:15

ceptree @ceptree

ヴォルテラ積分方程式とヴォルテラ積分微分方程式のコロケーション法による数値計算のコードをPythonからJuliaに移植してた

タグ：

posted at 01:09:17

Stanford University @Stanford

A motorized exoskeleton that attaches around the ankle and foot can drastically reduce the energy cost of running, researchers in Stanford's biomechatronics lab found. stanford.io/3bEBm1l pic.twitter.com/IIsZMaYsAy

タグ：

posted at 03:01:27

#数楽こんな深夜に何をやっているか。ググれば大量に情報が出て来ると思っていたが、そうでもなかった不完全ガンマ函数

P(a,z) = γ(a,z)/Γ(a)

のaに関する偏導函数の数値計算法。

ググったのは

derivative of incomplete gamma function

見つけた情報をメモっておきます。続く

タグ：数楽

posted at 03:06:45

#数楽 P(a,z)はガンマ分布Gamma(a,1)のcdf(累積分布函数)です。ガンマ分布のcdfにautomatic differentiationを作用させるために、その偏導函数の情報が欲しいと思いました。

これがモチベーション。続く

タグ：数楽

posted at 03:09:19

#数楽 γ(a,z)=∫_0^z e^{-t} t^{a-1} dt = Γ(a) - Γ(a,z) のΓ(a)の導函数はdigamma函数ψ(a)とΓ(a)の積になるので、Γ(a,z)のaに関する偏導函数を知りたい。

答えは

functions.wolfram.com/GammaBetaErf/G...
↓
functions.wolfram.com/GammaBetaErf/G...

に書いてあります。

タグ：数楽

posted at 03:15:45

#数楽使えそうな公式を添付画像として引用。

2Ftilde2(a,a;a+1,a+1;-z)の部分は一般超幾何函数です。tildeの意味はΓ(a+1)²でtildeなしのやつを割った結果です。その辺の記号の意味は

functions.wolfram.com/GammaBetaErf/G...

を見ながら、Wolfram言語のドキュメントを検索して確認しました。 pic.twitter.com/uzI2WGwDQa

タグ：数楽

posted at 03:20:41

#数楽 Γ(a,0,z)はγ(a,z)のことです。

問題1: この公式は信用できるか？

問題2: 一般超幾何函数2F2を #Julia で数値計算するにはどうすればよいか？

問題1の答えは詳しい説明は省きますがYesです。

問題2の答えは

github.com/JuliaApproxima...

にありました。

タグ： Julia 数楽

posted at 03:24:11

#数楽以上を読めば、#Julia言語でガンマ分布のcdf (累積分布函数)のパラメータに関する偏導函数を計算する方法が一応分かるということになります。

一般超幾何函数mFnの部分は

github.com/JuliaApproxima...

で計算できる。

タグ： Julia言語数楽

posted at 03:27:02

#数楽 #Julia言語現在のDistributions.jlでのガンマ分布のcdfはStatsFuns.jlで定義されているgammacdfを使っており、それはRmath.jlから来ています。

gammacdf(a, 1, x) = P(a, x) = γ(a, x)/Γ(a)

です。P(a, x)はSpecialFunctions.jlを使って

P(a, x) = gamna_inc(a, x, 0)[1]

でも計算可能。

タグ： Julia言語数楽

posted at 03:31:16

#数楽 #Julia言語ガンマ分布のcdfのパラメータに関する偏導函数の数値計算ぐらい楽勝だろうと思ってやろうとしたら、超幾何の2F2が必要になって「油断し過ぎていたこと」を思い知らされました。

一般超幾何の高速数値計算は実用面でも結構重要ということのようです。

タグ： Julia言語数楽

posted at 03:33:49

#数楽 #Julia言語 Distributions.jlにおける関連issue

github.com/JuliaStats/Dis...
Limited Support for ForwardDiff on CDFs #638 [Open]

タグ： Julia言語数楽

posted at 03:45:18

buvery @buvery

米国の新規失業者が爆増加傾向しているとの図。米国は、雇用が特に不安定なサービス業が新型肺炎の自粛や封鎖で直撃され、莫大な失業者を産み出している。戦争はもう始まっているので、日本も早く対応しないとヤバイよ。 twitter.com/hirosetakao/st...

タグ：

posted at 04:31:11

なないお @Nanaio627

金持ちには配らなくていいじゃんって意見も見かけるし気持ちがわからないでもないけど、世の中金持ちの方が圧倒的に人数少ないし、配る人を選別するための手続きや作業にかかる経費や人件費の方がよっぽど高くつくと思うんだよね。

タグ：

posted at 05:19:34

なないお @Nanaio627

なので対象者を絞っても1人が貰える分が増えるわけでも国の支出の節約になるわけでもないと思うの。かえって負担増えるよ。

タグ：

posted at 05:21:43

田中秀臣 @hidetomitanaka

こんな経済危機でも財界は増税主義者を日銀に送り込もうとしている。それを後押しする財務省と問題意識が微塵もない国会。委員の大勢に従うだけの高給取りを生み出すだけ恥ずべき選択だ

コロナ不況でも消費増税？お粗末すぎる日銀の「族委員」 (田中秀臣) ironna.jp/article/14608 #iRONNA

タグ： iRONNA

posted at 05:26:31

積分定数 @sekibunnteisuu

@amiami114114 @akita_kia togetter.com/li/1241319　「算数教育の専門家」がろくでもないことになっている訳で、「算数の専門」という教師は、掛け算の順序も足し算の順序も拘ることになるでしょう。

タグ：

posted at 05:45:36

Ivan Danyliuk @idanyliuk

pic.twitter.com/qMAFITkF2j

タグ：

posted at 05:47:36

積分定数 @sekibunnteisuu

@amiami114114 @akita_kia 「細菌の若い教師は等分除も包含除も知らない」と嘆く、筑波大学附属小学校算数部の田中博史氏
tsukuba-sansubu.cocolog-nifty.com/blog/2010/08/p...

マナー講師や宗教と同じで、非合理な教義が事細かく決まっていた方が「専門家」の権威が保たれる。

タグ：

posted at 05:51:08

#jupyter #termux #ubuntu-in-termux Termuxにストレージ権限許可してあるからJupyter Notebookから/storage/emulated/0/にアクセスできる。 pic.twitter.com/aqV7mpa63x

タグ： jupyter termux ubuntu

posted at 08:50:49

他のアプリからもアクセスできる場所だからなんか便利になりそうな気がする。

タグ：

posted at 08:54:35

#jupyter #termux #android ちなみに/storage/と/storage/emulated/0/にはアクセスできますが/storage/emulated/にはアクセスできません（1枚目が/storage/、2枚目が/storage/emulated/、3枚目が/storage/emulated/0/）。（9EBE-170Bは多分マイクロSDカード） pic.twitter.com/Tt2Xe5fwIZ

タグ： android jupyter termux

posted at 09:10:12

#jupyter #julialang（この辺がJupyter Notebookでいじれるのは怖いので実際には0の下にサブディレクトリを作ってそこで起動してます

julia> notebook(dir="/storage/emulated/0/サブディレクトリ")
）

タグ： julialang jupyter

posted at 09:19:24

#termux #julialang JupyterじゃなくてJuliaのコードを書くときは、termuxのフォントサイズの変え方がわからないので、変え方がわかってるDroidVimを使いたい。でtermuxからスマホのストレージにアクセスできるとDroidVimで編集したファイルを直接実行できて便利。

タグ： julialang termux

posted at 09:32:16

（julia-vim github.com/JuliaEditorSup... がどう考えてもDroidVimには対応してないけどVim scriptで書かれてない部分をVim scriptで書き直せば多分大丈夫でしょ。 #簡単に言うんじゃありません） #julialang #vim #droidvim

タグ： droidvim julialang vim 簡単に言うんじゃありません

posted at 09:38:53

#統計モデル内に現実から得たデータの情報をを持ち込む方法は複数ある。

* モデル内確率分布を現実と同じデータがモデル内で生成された場合に制限した条件付き確率分布を考える。

* モデル内確率分布が現実から得たデータを生成する確率を最大化するようにモデルのパラメータを調節する。 twitter.com/genkuroki/stat...

タグ：統計

posted at 09:58:36

#julialang #Jupyter 1.4.0は普通に入れられました。前はbuildを忘れてたみたい😅
twitter.com/MathSorcerer/s...

タグ： julialang Jupyter

posted at 09:59:52

名取宏（なとろむ） @NATROM

[布マスクはないよりマシなのか？ natrom.hatenablog.com/entry/2020/03/... ]。布マスクは「ないよりマシ」とは言えない。つけたい人がつけるのはいいが、文部科学省が普及を呼びかけるようなものではない。

タグ：

posted at 11:18:13

フォントサイズ普通に変えられるじゃないですかsites.google.com/site/fudist/Ho...
twitter.com/brv00/status/1...

タグ：

posted at 12:00:34

#termux ゲームの実況動画とか見てると「…はできないのか」みたいなこと言ってて実際はできるけど少し位置がずれてるだけというようなことよくある。（訳:ピンチインは最初に試したときはできませんでした）

タグ： termux

posted at 12:06:27

（でもDroidVimでjulia-vimの機能が使えるようにはしたいけどね）

タグ：

posted at 12:10:24

#数楽 #Julia言語コードを見たい人がいると思うのでノートブックを作りました。

nbviewer.jupyter.org/gist/genkuroki...
ForwardDiff.gradient of StatsFuns.gammacdf

Forward.Dual{T}型対応していない函数を自動微分対応させる方法の実例としても貴重かもしれません。

タグ： Julia言語数楽

posted at 12:35:35

Nandesukakoreha😵TeXdakejanaindesuka #julialang pic.twitter.com/UbTBzWEuo9

タグ： julialang

posted at 12:52:08

闇のapj @apj

@sugikota 私的な海外旅行を届け出ろというのはプライベートに干渉しすぎじゃないんですかね。そもそもの規則がおかしい。なんでこんなプライバシー侵害を許しておくのか。

タグ：

posted at 12:52:59

ゆべし @yubessy

柳　時熏【囲碁プロ】【YouTuber @ryu_shikun

暇な時に言語やライブラリのドキュメントを読む趣味があるんだけど、最近このおかげでJuliaのvalue typeを導入できる箇所を仕事で発見した

タグ：

posted at 14:40:58

昨日昼休憩の時、色々な棋士と話をしましたが、ニコちゃんは難しすぎるとの事。

動画を見る人も半分が20代から30代。
どうやら若手棋士に愛用されているらしいですね。
もう少し簡単な問題を作ってみます。

今日の問題です
www.ryu-igo.com/entry/2020/03/...
www.ryu-igo.com/entry/2020/03/...

#囲碁　 #ニコちゃん問題

タグ：ニコちゃん問題囲碁

posted at 16:02:27

A(うすいオレンジ)の国はBCGワクチンを定期接種、B(紫)の国は以前していたが今は接種せず、C(濃いオレンジ)の国は定期接種した事がない。

米国・イタリアはC、欧州もBが多くを占めている。とはいえ、中国はあれだけ新型コロナウイルス流行したけれど、日本と同じAです

www.researchgate.net/publication/50... pic.twitter.com/nvp2B7fPGv

タグ：

posted at 16:14:42

あおじるPPPP @kale_aojiru

他の生徒と比較してもよく学んでいる生徒が「うちで全然勉強しないでスマホやパソコンばっかり見てる」と保護者に言われていました twitter.com/TJO_datasci/st...

タグ：

posted at 16:19:48

非公開

タグ：

posted at xx:xx:xx

#統計「尤度」(ゆうど、likelihood)は「モデルのサンプルへの適合度の指標」として使われます。

モデルがサンプルにぴったり適合していると、尤度は高くなりますが、モデルの通常の意味でのもっともらしさは下がるので、尤度は通常の意味でもっともらしさを意味しません。

詳しくはスレッドを参照。 twitter.com/noritokuta/sta...

タグ：統計

posted at 16:22:14

非公開

タグ：

posted at xx:xx:xx

実はBCGワクチンと言っても、色々株があります(www.ncbi.nlm.nih.gov/pubmed/19540166)。

日本は東京172、中国はデンマーク株の亜型の上海株です(www.ncbi.nlm.nih.gov/pmc/articles/P...)。イランはパスツール1173株で、隣国イラクは東京172株で、イランに比べ流行数は抑えられています(www.worldometers.info/coronavirus/)。 pic.twitter.com/cEGTrd9mD2

タグ：

posted at 16:37:22

当然『BCGワクチン接種＝新型コロナウイルス感染予防』は短絡的過ぎますし、単なる相関関係の可能性はありますが、臨床研究で効果が判明すれば、我々小児科医が行っていた事が少し役に立てている事になり、あの大変な接種をした甲斐があります。

タグ：

posted at 16:42:02

TaKu @takusansu

コミさん@K Squad LLC. @komi_edtr_1230

#超算数算数・数学を「好きになって欲しい指導」と「嫌いにならないでという指導」は別物だと思う。
「好きになって欲しい指導」と「理解して欲しい指導」は一致しているとは限らない思う。

タグ：超算数

posted at 17:09:54

Boxプロットができるようになったよ！！！
github.com/komi1230/kai pic.twitter.com/a0oDNHdhyy

タグ：

posted at 17:29:00

3pF @3pF

Misaki Ohta(太田　三砂貴) @misaki_ohta_

受けなくてもいい甲状腺検査を受けて、しなくてもいい摘出手術をして、一生甲状腺ホルモン薬を飲み続けなければいけなくなる羽目になる人もいるのに、それが「福島の人たちを傷つけないでほしいのだ」は、ちょっとないかな。 twitter.com/arai_fukushima...

タグ：

posted at 18:08:37

sinの無限積展開の収束の様子 #Julialang #plot pic.twitter.com/dANEPxjfWC

タグ： Julialang plot

posted at 18:12:10

『接種率で分けると信憑性』

2016年データ(www.who.int/immunization/s...)ですが、WHOの各国接種率です。

実は中国より日本の方がBCG接種率は低く(8割)、岩田健太郎のように『娘にBCG接種しない』派が多いんですかね？接種率は中国はよいので、株の違いなのかな？という印象です

twitter.com/nodonodo1109/s... pic.twitter.com/8fsJk3yqnB

タグ：

posted at 18:43:14

『結核に関しては』10年程と思いますが、新型コロナウイルス感染はまた別になるので、何か別の免疫機序なのかな？と。ま、あくまで仮説なので、臨床研究結果待ちです。

twitter.com/sunring_rise/s...

タグ：

posted at 18:45:09

Kiichi @Ki_chi

蛇蔵@天地創造デザイン部8巻発売中 @nyorozo

Julia本、ポチった

タグ：

posted at 18:47:14

蛇蔵@天地創造デザイン部8巻発売中 @nyorozo

「手を洗えば感染の予防になる」
これを最初に唱えた医者が、まったく信用されず亡くなった話（1/2）
この2ツイートだけでラストまで読めます。現代科学のねもとには、彼のような苦難の一里塚があったことを知ってほしくて書きました。 pic.twitter.com/I6h6peRMEn

タグ：

posted at 19:02:50

これを最初に唱えた医者が、まったく信用されず亡くなった話（2/2） pic.twitter.com/EflScl41VZ

タグ：

posted at 19:03:18

bao @baobabustroll

@sugikota もしもこの女性教授が「いたたまれないから懲罰を受けたい」と言ったとしても、大学は外部から見て懲罰的に見える対応をこの教授に対して一切行うべきでは無いと私は考えます。
理由は「感染者に対する感染を理由にした懲罰の前例」を作ってはならないからです。
これは、私刑の容認へとつながります。

タグ：

posted at 19:07:52

optical_frog @optical_frog

通常の景気後退でやる景気対策と，今回のパンデミックでの経済的支援は，別物だというお話．ニューヨークタイムズのコラムでも出てきた．www.nytimes.com/2020/03/26/opi...

タグ：

posted at 19:21:34

optical_frog @optical_frog

大意：ふつうの景気後退なら，失業してしまった人たちがふたたび仕事に就けるようにマクロ経済政策をはたらかせる．他方で，パンデミックでは一時的な失業が対策の一環として必須になる．政府は，この失業の痛手が人々の生活を壊さないようにお金で支援しないといけない．twitter.com/paulkrugman/st...

タグ：

posted at 19:24:10

地下猫 @tikani_nemuru_M

マスクはほとんどの場合、呪術的武装か、配慮していますアピール以上の意味はないと考えているが、職場で上司が血走った目で「マスクをつけろおお」と叫ぶのとケンカする気力もないのである。 / “布マスクはないよりマシなのか？ - NATROMのブログ” htn.to/3YCRCPvaEt

タグ：

posted at 19:34:10

l_ppp @ppp3141592ppp

1から始めるJuliaプログラミング
発送されないな～ pic.twitter.com/IqDmbcYmj6

タグ：

posted at 19:53:06

PIXSEL@巴拿马文件 @pixselP

@nyorozo かのナイチンゲールが近代医学を推し進めようとしたときも
旧来の慣例主義、利権、身分主義で幾度となく横槍があり、阻まれてきたといいますし、
それを片っ端から力づくでへし折ってきて今の医学の勝利があったとも。

タグ：

posted at 20:26:21

optical_frog @optical_frog

※ニューヨークタイムズのコラムとあわせて参照してくださいまし．www.nytimes.com/2020/03/26/opi... pic.twitter.com/PWdmRH2czg

タグ：

posted at 20:41:46

ふぁっふぉい @sugikota

Misaki Ohta(太田　三砂貴) @misaki_ohta_

@masahirono しかしながら患者本人の利益にならない検査はヘルシンキ宣言に反するので、公衆衛生上の必要性でそれをするには新たな法的根拠が必要と思います。患者全体の利益と患者個人の利益を厳格に分けるヘルシンキ宣言の理念の限界が改めて明らかになったと感じます。

タグ：

posted at 21:21:43

pic.twitter.com/oaURmvyFSo

タグ：

posted at 21:36:30

gearvoice @gearvoice1

@baobabustroll @sugikota 実に興味深い
「心ない人たちがいる。その行動は
リンチだ」との記事で
ある人は責任に目が行き
関知できないルールを問題視し

ある人は個人・関係者の心理に目が行き
卒業式の中止といった事象での
心的負荷を心配し

ある人は感染に目が行き
感染を理由にした懲罰と受け取られることを危惧している

タグ：

posted at 21:44:24

gearvoice @gearvoice1

@baobabustroll @sugikota そして共通しているのが
それぞれが興味を向けたのが
コントロールが
難しいモノであるというコト
その事実が面白い

タグ：

posted at 21:48:46

#超算数 1951年指導要領よりも後だが、アメリカの算数教育書に掛算の式と状況を結びつきを発見した。
Grossnickle, Foster E., Leo J. Brueckner. _Discovering meanings in arithmetic_. Philadelphia: Winston. 1959. archive.org/details/discov...
4×0は的当てゲームで4回の0点を取ることを表せるが、 pic.twitter.com/aANS9BXrCL

タグ：超算数

posted at 22:09:41

#超算数 0×4は状況を表わすことができないといっている(177頁)。日本式の指導法を採用するなら、2×5のような式を併用して、4点の的には一回も当てられなかったが、5点の的には2回当てたという状況を式で表現するところだ。ブルクナーとグロスニケルは別の共著を紹介済みtwitter.com/temmusu_n/stat...。

タグ：超算数

posted at 22:09:42

#超算数 1951年指導要領の一部が、この本に基づいていることは、関係者の証言によって確実である。1947年出版の本には掛順を見いだすことが出来なかった。しかし、今回みた1959年の共著には上のように式と状況の一対一の結びつきが見られる。とはいえ、それも日本の指導ほど頑健なものではない。例えば

タグ：超算数

posted at 22:09:43

#超算数 187頁の九九表では行と列のどちらが乗数でどちらが被乗数なのか指定がなく、積を作る数はどちらも因数と呼ばれている。また、同じ累加を表わす掛算の式が縦書きと横書きで乗数の位置が変わるし(181頁)、176頁は36×7を日本と同じ乗数の位置で読ませている。 pic.twitter.com/qIeezCmakP

タグ：超算数

posted at 22:09:44

#超算数引算の使われる状況、用法として求差、求残など4つを紹介している(156頁)ことは日本の指導法とよく似ているが、足し算においては増加と合併を区別することはしていなかった。当時の日本ではすでに足順の芽生えが確認できるのだが。 pic.twitter.com/jkqJAfvwd6

タグ：超算数

posted at 22:09:45