黒木玄 Gen Kuroki(@genkuroki)/2021年01月/Page 26

あとは、データプロット系の機能も充実してるようだった。これは全く触ってないのであまりわからないけど。
ゴリゴリの数値計算をしてから結果をDBに放り込んで加工してプロット、という流れがシームレスにできるのは気持ちよさそう。
#Julia言語

タグ： Julia言語

posted at 00:03:04

ʇɥƃıluooɯ ǝıʇɐs @tsatie

@genkuroki ちょっと見難いかもしれませんがダメでした... pic.twitter.com/dEHgMuTtPY

タグ：

posted at 00:14:10

ʇɥƃıluooɯ ǝıʇɐs @tsatie

@genkuroki カーネルに接続できないという事です...

タグ：

posted at 00:15:25

p*@stastastastastast @p_star__0per

Juliaのスッキリしたコード見てからどうもPythonのが汚く感じる

タグ：

posted at 04:36:55

イベルメクチンのmeta-analysisで死亡率75％低下が報告されました。forest plotに介入群とコントロール群の情報を加えたのが以下。期待が持てる結果ですが、NIHはまだ推奨に至らず。理由の一つは介入群や対照群が受けた治療が様々なため。1/5
www.researchsquare.com/article/rs-148... pic.twitter.com/NPqIDBHUZU

タグ：

posted at 05:32:02

各研究の参加者数が少ないのも理由の一つ。ランダム化は介入群、対象群の背景因子の偏りをなくすためだが、参加者数が少ないと偶然で偏りが出る場合がある。劇的な効果が見られたElgazzarの背景因子をみるとHydoroxychloloquineに割り付けられた重症者で虚血性疾患などの割合が多いなど偏りが。2/5 pic.twitter.com/hf2lP4xEcy

タグ：

posted at 05:32:03

この偏りだけで劇的な死亡率の差は説明はできないと思いますが、やや懸念が。meta-analysisに加えられた研究全てを合わせても死亡者数は71人。GRADEグループが指摘するように、小規模RCTで示された劇的な結果が後の大規模RCTで否定された例が過去にも散見されている。linkinghub.elsevier.com/retrieve/pii/S... 3/5

タグ：

posted at 05:32:03

そこでイベント数が少ないと、信頼区間の推定値が vulnerableなので統計的に有意であっても、より現実的な治療の効果を仮定して有意差を検出するのに必要なサンプルサイズ（Optimal Information Size）を計算し、これを満たしていたかを考慮することが必要。これは死亡率が50％減と仮定しても　4/5

タグ：

posted at 05:32:03

河北新報オンライン @kahoku_shimpo

検出力90％でもほぼ達成している模様（50%の妥当性は置いておく）。あとは各研究にバイアスがあったり論文の記載が不十分でバイアスの有無が判断しがたい例が多い。ということで期待が持てるが、もう一つ確信が持てない理由はこんなところだと思う。5/5

タグ：

posted at 05:32:03

とんこつラーメン「一蘭」東北初出店　仙台・名掛丁に２８日、感染防止徹底 kahoku.news/articles/20210...

タグ：

posted at 06:02:10

Cookie Monster @MeCookieMonster

Me no geologist, but me think dat rock look a lot like me... twitter.com/ImoveCar/statu...

タグ：

posted at 06:42:57

Y Tambe @y_tambe

mRNAのキモは、生ワクチンと同様に、細胞性免疫も誘導できること。ウイルス感染防御には、これがかなり効く。

タグ：

posted at 07:40:03

@tsatie 結局、役に立てなくてごめんなさい。

どうして1.6.0-beta1だとJupyterでカーネルが死ぬんですかね？

タグ：

posted at 07:50:55

#数楽統計学で出て来やすい確率分布達について「曼荼羅」的な「整理」をするのってもしかして流行っている？

個人的にはあんまりお勧めしたくない考え方。

2次元平面上の「曼荼羅」は確率分布達の整理において不適切な世界観だと思う。

2次元じゃ足りないし、雑な見方になり易いと思う。

タグ：数楽

posted at 09:09:43

#数楽本気で整理したいなら以下のように考えたい。

いきなり「曼荼羅」を描いたりしない。

確率分布(達)から他の確率分布を作る沢山の方法を整理したい。直積、制限(条件付確率分布)、射影(一般には像、積分)、極限、…。

そして、タネになる基本的な確率分布達として何が取れるかも整理したい。

タグ：数楽

posted at 09:09:43

#数楽分布の制限の例

r個のポアソン分布の直積をx_1+…+x_r=nに制限すると、(条件付き確率分布として)多項分布が出て来る。

4個のポアソン分布の直積(座標としてa,b,c,dを取る)をa+b,c+d,a+cが一定の部分に制限すると超幾何分布が出て来る。

この意味でポアソン分布は離散分布の「基本的な部品」。

タグ：数楽

posted at 09:09:44

#数楽球対称な分布(例えば標準正規分布の直積)を、球面上に制限したり、半径方向に積分したりすれば、球面上の一様分布を作れる。

次のツイートで説明するMaxwell-Boltzmann則をコンピュータによるモンテカルロ法で数値的に確認したい人にとって、これは役に立つ知識。

randn(n)をノルムで割る。

タグ：数楽

posted at 09:09:45

#数楽 nが大きいとき、n次元空間中の球面上の一様分布を部分空間に射影してできる分布はその部分空間上の正規分布で近似される(Maxwell-Boltzmann則)。

直積、射影、極限(近似)と色々出て来ている。

正規分布はこういう形でも出て来る。

タグ：数楽

posted at 09:09:45

#数楽共役事前分布関連の知識も有用だよね。

そういう諸々のことを「曼荼羅」に描いても他人には分かりやすくはならないと思う。

しかし「曼荼羅」を描いた人だけは多分非常に勉強になっている(笑)

だから、必ずしも適切な整理の仕方ではないということを承知で描いてみるのはありかもしれない。

タグ：数楽

posted at 09:09:46

これ、大事な問題。

例えば私のような人間が頭をフル回転させる楽しみを失ったときにどうすればよいのか？

これについて私は「かわいいを探すことを楽しみにすればよい」と考えています。😊😊😊

私のTLの動物動画の割合は結構多めだと思う。

この世には楽しいことは沢山ある。 twitter.com/__dalong/statu...

タグ：

posted at 09:16:56

trichoides @benzyl_acetate

Julia set
f(z) = z^2 + C
C = - 1.29

グラデーションがうまくいったと思う
#切り絵 pic.twitter.com/WQypEU3l5L

タグ：切り絵

posted at 10:00:51

非公開

タグ：

posted at xx:xx:xx

#Julia言語

高次元球面上の一様分布の低次元部分空間への射影は正規分布でよく近似される(Maxwell-Boltzmann則)。

たった数行のコードで確認できる。

球面への制限は物理的にはエネルギー保存則で、高次元化は熱浴を考えることに対応している。

gist.github.com/genkuroki/a917... twitter.com/genkuroki/stat... pic.twitter.com/SmaGzIKSOy

タグ： Julia言語

posted at 10:19:04

Hiroo Yamagata @hiyori13

日本のシュンペーター通俗解説書って、「コンドラチェフとか波動マンセー」のトンデモ景気循環とか、「イノベーションなんて余裕ある大企業しかできないよねー中小企業ゴミクズベンチャー死ね」論とか、都合の悪いところはほっかむりしていいことになってるのか。

タグ：

posted at 11:37:15

#数楽一般に○○の(逆温度の函数としての分配函数やカノニカル分布の密度函数の)Mellin変換(=○○のゼータ函数)を考える御利益はMellin逆変換の式を書けば分かる。

Mellin逆変換は積分路が虚軸に垂直な複素積分になるので、もとの○○側の漸近挙動の解析ではゼータ函数の極の情報が有用になる。 twitter.com/genkuroki/stat...

タグ：数楽

posted at 11:53:13

チョー算数問題についてここ10年語って来た人達は多くの小学生とその保護者を救って来たと思う。

教科書の指導書の問題点に関する情報は、特に小中学生の保護者達や中高生にとって緊急かつ重要な情報になり得ます。

内輪でしかその話をしないようになると、「嫌がらせ」が成功したことになります。 twitter.com/f_sei/status/1...

タグ：

posted at 12:15:56

一部の人しか閲覧困難な教科書の指導書に有害なことがずっと書かれてきたことがばれると本当に困る相手は、何も知らない一般人でしょう。

特に小中高生に直接知られてしまうのは痛いでしょう。

嫌がらせに屈すると、情報を必要としている人達に緊急かつ重要な情報が届かなくなります。 twitter.com/f_sei/status/1...

タグ：

posted at 12:15:57

現在教育を受けている小中高生にとっては今すぐに問題が解決してくれなければ、被害を受けっぱなしになってしまうんですね。

私もチョー算数問題については長期的な視点を述べたときにそういうことを言われて考えさせられました。

問題を解決できなくても、問題を緩和するための情報は拡散できます。

タグ：

posted at 12:15:58

#超算数|タグの常連たちが、教師用指導書について積極的に議論した成果の一つがこの秀逸なまとめです。
「かけ算の順序」なんてもう古い⁈　今や時代は「足し算の順序」‼︎ togetter.com/li/901635
教科書通りに算数を習うと、2+6=8, 6+2=8に区別があるとか、増加と合併は違うとか信じてしまうかも。

タグ：超算数

posted at 12:30:44

#超算数更に重要なこと。教師用指導書は現場で物凄く重宝されているだけでなく、教育委員会レベルで使用が推奨されています。【教科書が宝石箱なら、教師用指導書は、その宝石を磨く上質のクロスなのである】www.edu-ctr.pref.okayama.jp/chousa/study/k...と2007年に書いたのは岡山県教育センターtwitter.com/temmusu_n/stat...。

タグ：超算数

posted at 12:38:15

全卓樹 @Quantum_Zen

離散化反応拡散方程式２ビット版でチューリング的キリン模様出したやつ、パラメータをローカルな発色量に依存して変化させたら（一種の場の理論になるけど）、分裂生成する細胞のコロニーの成長と破壊みたいな「人工生命（ALife）」できたぞ。こんなシンプルなモデル誰かがもおやってるとは思うけど... pic.twitter.com/9lRHgzxDxQ

タグ：

posted at 13:13:23

#超算数教師用指導書の研究は、人の恥を暴くというスキャンダルとかゴシップ報道のためにやっているのではありません。ひどい内容の文書があり、かつそれが全国津々浦々の教室で使用されているのです。その内容を理解してどう対処するか議論するべきです。

タグ：超算数

posted at 13:40:58

全卓樹 @Quantum_Zen

@nabeso 実際これ2次元オートマトンです。ライフとの違いは、恣意的ルールではなく由来が自然な結合反応拡散方程式系な事と、メッシュ細かいとゆうか相互作用が長距離(この例だと8セルまで及ぶ)な事です。

Mami Tanaka先生とNabeso先生のツッコミのせいで、ツイタに投げた学生演習の題材からここまで発展...

タグ：

posted at 14:09:33

yudai.jl @physics303

fMRIの多次元配列データがいくつか欲しんだけど，ぐぐってでてくるデータセットが割と巨大で，どれをダウンロードしようか迷っている．

しばらくディープラーニングやらないと，数GBのファイルをダウンロードするのに抵抗感がでるようになるな

タグ：

posted at 15:33:05

yudai.jl @physics303

これが軽くて良かった．using NIfTI で開ける．というが，juliaさん，ライブラリ充実し過ぎでは？

exhibits.stanford.edu/data/catalog/x...

タグ：

posted at 15:56:11

yudai.jl @physics303

あー，これ.matファイルかー．Juliaじゃ開けんなー
⇒ using MAT で解決

あー，これNIfTI formatかー．Juliaじゃ開けんなー
⇒ using NIfTI で解決

あー，これGeoTIFF formatかー．Juliaじゃ開けんなー
⇒ using GeoArrays で解決

Juliaさん，ライブラリ充実具合すげぇな

タグ：

posted at 15:58:43

ちょうど研究で使う機会があったので#julia言語でshooting methodを実装してみました。
二次の非線形光学媒質中でのソリトン解(quadratic optical soliton)を数値的に求めています。
#julialang 1/n pic.twitter.com/jSLo3booWF

タグ： julialang

posted at 17:14:34

shooting methodでは、u(z)とv(z)がz→∞で0になるような孤立解を探すために、
適当なz=0での初期値に対して微分方程式を解く→十分大きなz=z0に対してu,vの値を得る→二分法でu(z0)=v(z0)=0になる初期値を探す
と言った実装になります。(微分方程式はDifferentialEquations.jlで解いています)
2/n

タグ：

posted at 17:14:34

ここではρ=1の場合を解いていますが、ρの値によってはいくつかの山のある高次の解(multi-hump soliton)に捕まってしまったりするので、傾きや曲率に追加の条件を入れて強引に山一つの解に収束させてます。（もっと賢いやり方ありそう、、）
3/n

タグ：

posted at 17:14:34

ちなみにここでのuとvはそれぞれ周波数fの基本波と周波数2fの二次高調波に対応しています。最もよく見るKerr solitonと比べると知名度が低めですが、個人的には実用上quadratic solitonはとても面白いと思っているところです。
4/n

タグ：

posted at 17:14:35

Reference
www.sciencedirect.com/science/articl...
5/n

タグ：

posted at 17:14:35

journals.aps.org/pre/abstract/1...
6/n

タグ：

posted at 17:14:35

Re: RTs #Julia言語

どうしても微分方程式のソルバを自分で書きたい人以外には、DifferentialEquations.jl github.com/SciML/Differen... は非常にお勧め。

RやPythonでも使えます。

ハミルトニアンを与えて正準方程式を解かせることもできるし、多彩な数値解法が実装済みで、衝突なども記述可能。

タグ： Julia言語

posted at 17:29:58

DifferentialEquations.jlはほんとそれだけでも #julialang を使う価値がありそうなライブラリだと思います！（多機能すぎてprecompilationが永遠に終わらないこと以外は、、） twitter.com/genkuroki/stat...

タグ： julialang

posted at 17:39:08

#Julia言語 v1.5.xでは

pkg> add DifferentialEquations
pkg> precompile

でインストールするとよいです。こうしておいた方が、using ～のときにプレコンパイルを待つより精神的によいです。

v1.6.0-beta1では

pkg> add DifferentialEquations

だけで、全力で(笑)プレコンパイルしてくれます。

タグ： Julia言語

posted at 17:46:15

yng @yana_phys

Julia 1.6 (beta) 触ってみたけど plot() がO(1) secでできるようになっていて感動した

タグ：

posted at 17:48:24

universato @universato

PythonとRuby、mod powが標準であるの便利すぎる。

タグ：

posted at 17:53:09

universato @universato

「mod powが標準である言語って何があるかな、juliaかな」って探したら、juliaのpowめちゃくちゃバリエーションに富んでそう
docs.julialang.org/en/v1/base/mat...

タグ：

posted at 17:57:27

#超算数教師用指導書の研究は例えば、足順こだわりの理解に不可欠。学校図書は2015年に【合併と増加の場面性の違いを強調したい。そのことが、加法の問題作りや絵本作りにおいて、より深まりのある問題作りに発展するものと考える。】twitter.com/temmusu_n/stat...と指導書に書いた。教科書にない記述。

タグ：超算数

posted at 18:07:56

ひのみ @HlNOMl

Julia Python R

合わせてJupyter

タグ：

posted at 18:09:42

#超算数

「算数教育の専門家」がろくでもないから心配。

小学校5・6年生に「教科担任制」本格導入を中教審が答申 | NHKニュース www3.nhk.or.jp/news/html/2021...

タグ：超算数

posted at 19:21:30

#超算数　「専門は数学」の現役小学校教師 pic.twitter.com/NSmHkEbjzj

タグ：超算数

posted at 19:23:30

@kuri_kurita

一見朗報のようでいて、チョー算数ウォッチャーとしては…

何しろ日本の場合は「算数教育」の専門家たちが徹底的に狂っているので、その狂った教義を無批判にクソ真面目に専門的に習得してしまった精神的奴隷のようなゾンビたちによる子供への被害が拡大するのではないかという恐怖も…😱 twitter.com/nhk_news/statu...

タグ：

posted at 19:24:00

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

#Julia言語私の環境ではプレコンパイル済みのDifferentialEquatons.jlについて

julia> @ time using DifferentialEquatons.jl

にかかる時間は

35.1秒 ← v1.5.3
10.5秒 ← v1.6.0-beta1

でv1.6.0-beta1が圧倒的に速いです！

タグ： Julia言語

posted at 19:54:47

非公開

タグ：

posted at xx:xx:xx

#Julia言語ハミルトニアンだけを与えて、正準方程式を DifferentialEquations.jl で解くことの例。

オープン戸田格子を解かせています。ついでに、離散オープン戸田格子も解いて、解がぴったり一致することも確認。

nbviewer.jupyter.org/gist/genkuroki...
連続オープン戸田格子と離散オープン戸田格子の解の一致 pic.twitter.com/Ss8W26SkL2

タグ： Julia言語

posted at 19:57:30

#Julia言語添付画像のプロットの仕方は、色々参考になると思います。

①DifferentialEquations.jlで得た解のプロットを行っている。正準方程式の解なのでpとqを別々にプロット。

②離散版の解と連続版の解のプロットの色を合わせている。

③ラベルをまとめて作っている。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/AMilE6EWpO

タグ： Julia言語

posted at 19:59:49

非公開

タグ：

posted at xx:xx:xx

#Julia言語 DifferentialEquations.jlのAPIはあるときに大幅に変更(breaking change!)したので、2年くらい前のコードはそのままでは動かないと思います。

Juliaでは検索して見つけた古いコードは大抵そのままでは動きません。過去1年以内のような制限を付けて検索した方が安全。

タグ： Julia言語

posted at 20:01:52

amaryllis @bylloop

なるへそ参考になる。 twitter.com/genkuroki/stat...

タグ：

posted at 20:03:04

#Julia言語オープン戸田格子は連続時間発展での時刻整数での値と、離散版オープン戸田格子の解がぴったり誤差無しで一致しています。

こういう例を知っていると、「正しく数値計算できているか」を確認するときに便利。初期値のパラメータの順番を逆にしてしまったりするミスはすぐに見付かる。

タグ： Julia言語

posted at 20:04:14

#Julia言語ハミルトニアンを与えるだけでDifferentialEquations.jlは正準方程式を数値的に解いてくれます。しかし、ドキュメントに詳しい説明がないような気がします。だから、この例は貴重かも。

ハミルトニアン函数Hの引数は

H(p, q, param, t)

の順序で4つになります。これ大事。 pic.twitter.com/IHjMEiJvwp

タグ： Julia言語

posted at 20:08:38

@musorami え、まさか、教科書ですか？（おそるおそる）

タグ：

posted at 20:09:40

非公開

タグ：

posted at xx:xx:xx

@musorami やはり、そこから残念なことになっているか…
非常にまずいですねえ…

タグ：

posted at 20:13:12

#Julia言語このハミルトニアンを与えるだけでよいというのは非常に便利です。

ODEのソルバも沢山実装されています。ユーザーは選び放題↓

diffeq.sciml.ai/stable/solvers...

タグ： Julia言語

posted at 20:13:17

早稲田治慶（本名）無制限メタバース専門家 @waseda_fablab

#Julia言語このDifferentialEquations.jlは超巨大パッケージで機能がものすごくたくさんあって、全体像をまったく把握できていません。

みんなでちょっとずつ紹介するといいかも。

タグ： Julia言語

posted at 20:14:40

非公開

タグ：

posted at xx:xx:xx

これ日本の理数系教育の敵、掛け算順序強制問題と同根では
twitter.com/musorami/statu...

タグ：

posted at 20:33:49

@musorami なるほどです。
自分は、現状認識を誤っていたようです。
このレベルだと、確かに教科書をどう使っても駄目そうですねえ…

タグ：

posted at 20:33:57

#超算数正しい推測。おおくぼさんtwitter.com/OokuboTact/sta...に文科省の委託事業でみずほ総合研究所が2007年に行った質問調査を教えてもらったのですがwww.mext.go.jp/a_menu/shotou/...、教員の94.3%が算数の指導は好きで、同じく83.8%が算数は指導しやすい教科だと感じている(ともに感じる、まあ感じるを合算)。 twitter.com/C_space630/sta... pic.twitter.com/keoHURDEXQ

タグ：超算数

posted at 20:53:17

#超算数この傾向は、同年ベネッセによって行われた質問調査でも明らかですberd.benesse.jp/berd/center/op...。何度も言及していますがtwitter.com/temmusu_n/stat...、算数の指導が得意と感じる教員は、他教科に比べてかなり多く、経験年数の浅い教員でも自信があるのが特徴です。 pic.twitter.com/yYT9pWIPYk

タグ：超算数

posted at 20:53:18

非公開

タグ：

posted at xx:xx:xx

@musorami 「たすのかな？　ひくのかな？」みたいなことやってたら、できるものもできなくなりますからね。
あと、子供がわかっている概念をこねくり回して、むしろわからなくさせているものも多数ありますし。
困るのが、こちらがちゃんと教えておいたのに、学校で上書きインストール…

タグ：

posted at 21:00:02

非公開

タグ：

posted at xx:xx:xx

#超算数　

ここに挙がっている事例の多くは、「アホ教師のアホ授業」じゃなくて、算数教育の専門家によるもの。
togetter.com/li/901635

タグ：超算数

posted at 21:10:58

最新のJupyter Notebookとnbextensionsの組み合わせでうまく動かない場合の解決法
【Jupyter Notebookのバージョンを6.1.5に下げる】
qiita.com/SolKul/items/4... #Qiita @ICR3_5より

タグ： Qiita

posted at 21:20:12

ʇɥƃıluooɯ ǝıʇɐs @tsatie

そんなこともあるのですね twitter.com/genkuroki/stat...

タグ：

posted at 21:24:34

非公開

タグ：

posted at xx:xx:xx

Jiahao Chen 陈家豪 @acidflask

Tagging some the actual people doing the work @AlanEdelmanMIT @SimonPJByrne @jakebolewski @JuliaLanguage #julialang @vchuravy twitter.com/MIT_CSAIL/stat...

タグ： julialang

posted at 22:46:30

thread: Sorting algorithms visualized
some famous algo visualized using python, numpy and matplotlib😋

first is Odd-even sort!

#Algorithm #python #visualization pic.twitter.com/Ay6hz7wl5h

タグ： Algorithm python visualization

posted at 23:01:23

Heap sort up next! pic.twitter.com/83yJYaWD6k

タグ：

posted at 23:01:26

Bubble sort 🎈🐡 pic.twitter.com/v5Xstmsww6

タグ：

posted at 23:01:30

Stooge sort! not so useful.. but pretty! pic.twitter.com/o3JaTkXN1j

タグ：

posted at 23:01:34

quicksort ⚡️ a classic pic.twitter.com/itEJKWbzCq

タグ：

posted at 23:01:37

Cocktail sort 🍸 pic.twitter.com/L2M9FRlKl6

タグ：

posted at 23:01:41

Comb sort pic.twitter.com/2TDwoveULF

タグ：

posted at 23:01:45

Shell sort 🐚 pic.twitter.com/hRRNiUM2Rr

タグ：

posted at 23:01:48

insertion! pic.twitter.com/jyULzrLjvs

タグ：

posted at 23:01:51

all the (Bad!) code used can be found here:
github.com/mportesi/sorti...

タグ：

posted at 23:14:22

#超算数非常識・非論理的・不合理なチョー算数の体系を作ったり、チョー算数を小学生に教え込むために教科書や教科書のマニュアル本(教師用指導書、一般人購入不可)を整備したりしているのは、その筋の専門家達。

教科担任制を導入するとチョー算数教育が強化され、算数教育は悪化するだろう。 twitter.com/nhk_news/statu...

タグ：超算数

posted at 23:20:23

#超算数算数の教科書やその指導書がひどいことを当たり前の話だと感じるレベルの人達だけで算数を教えることができるように、小学校に算数についてまともな考え方をできる人達を大量に送り込む政策を先行させないと、教科担任制によって算数教育は悲惨なことになる可能性があると思います。

タグ：超算数

posted at 23:23:39

#超算数掛算順序が逆なら理解していないとみなす教え方には100年以上の歴史があります。

伝統的におかしな分野になってしまっている。

100年以上にわたる人脈と十分に整備された社会的仕組みが子供達を傷付けるチョー算数を支えている。 twitter.com/genkuroki/stat...

タグ：超算数

posted at 23:29:16