黒木玄 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年01月04日(月)
It's happening :) Makie just got a new Figure API, which brings us much closer to completely merge it with MakieLayout:
discourse.julialang.org/t/ann-makie-up... pic.twitter.com/rI11zNW1Kf
タグ：
posted at 23:55:26
脳が完全な確率計算をしているとすると、判断・意思決定に関する古典的な知見の多くを説明できないが、このサンプリングモデルならそれが可能。サンプリングモデルで、可用性バイアス、代表性バイアス、係留と調整ヒューリスティックを統合的に説明できるかもしれない。
journals.sagepub.com/doi/full/10.11...
タグ：
posted at 23:01:39
脳は確率を表象してもいないし、確率で計算もしていない。代わりに、記憶か心的シミュレーションからのサンプルを引き合いに確率計算に近いことをやっている。
journals.sagepub.com/doi/full/10.11...
タグ：
posted at 23:01:38
「数学の世界がコンピュータの中に来たぞー！」が好き。
【AbstractAlgebra.jl】Juliaで代数学をやってみたいんじゃ② ～剰余環を使って√2やiをつくろう～｜みぽ zenn.dev/nekomath271828... #zenn
タグ： zenn
posted at 22:59:31
だってアレだけブチブチ電源切っても平気の平左で動くし、インバータとか高周波機器があってもナンノソノな代物がほぼラダー一択とか悲しすぎない？
C/C++でもRustでもElixirでもJuliaでも書けた方が絶対楽しいじゃん！
(「誰がメンテすんのー」とか「言語のしんらいせーは？」とかはおいといてさ)
タグ：
posted at 22:55:34
ついに完成！！！！
今度は
HTMLとCSSだけでマリオつくりました💡
jsはもちろん、画像も一切使ってません！
相当な数のanimationが入ってるので
重いのと僅かな不具合はご了承ください💦
CSSだけで1万行ありますww😂
ash-creator.com/hentai_css/mar...
#Bボタンは飾りです
#2021年も変態をどうぞよろしく pic.twitter.com/XtxKqEJfxi
タグ： 2021年も変態をどうぞよろしく Bボタンは飾りです
posted at 22:20:02
@s_nkd1024 Wolfram Language(Mathematica): 「macで数式とかグラフとか綺麗に出力したいなら呼ばれたきが。」(？w)
(去年、Wolfram Language + Jupyter でMathematicaの基礎を学習した cs.mkamimura.com/tags/%E3%83%8F... からそう思っただけですが…😅)
タグ：
posted at 22:09:48
面白かったのでアニメーションにしてみた. デッドヒート twitter.com/physics303/sta... pic.twitter.com/PYKa4RsuFz
タグ：
posted at 21:58:19
Julia言語のLibSndFileというパッケージを導入して、音ファイルを扱えるようにしてみた。滝の音とか加湿器とか紙を破く音とか、身の回りの色々な雑音スペクトルをちょっくら調べてみたい。
タグ：
posted at 21:58:15
@stephensenn In Fig.1 in this excellent paper, you place Fisher, E. Pearson, Neyman, Jeffreys and de Finetti on the two dimensional space.
Where do you place G.E.P Box on the diagram if possible?
Is he placed eveywhere as an eclecticismist?
Or, do we need another dimension?
-Yusuke Ono
タグ：
posted at 21:38:26
Julia v1.6 なら↓とか？
```
Threads.foreach(Channel(nthreads()) do chnl
for val in 1:N
put!(chnl, val)
println('T', threadid(), " => ", "put ", val)
end
end) do val
println('T', threadid(), " => ", "took ", val)
end
``` twitter.com/LirimyDh/statu... pic.twitter.com/dqXcD6aF41
タグ：
posted at 21:18:55
#超算数 シェーンフェルドさんの問題提起は大量にRTしたSoraさんの撮要が分かりやすいです。なお原著
Schoenfeld, A. H. (1988). When good teaching leads to bad results: The disasters of "well-taught" mathematics courses. Educational Psychologist, 23(2), 145–166. jwilson.coe.uga.edu/TiMER/Schoenfe...
タグ： 超算数
posted at 20:34:15
#超算数 アラン・H・シェーンフェルドさんは日本の算数・数学教育界でも知られた研究者www.google.com/search?q=%E6%9...。ごまかし勉強の藤沢伸介さんも算数・数学教育界からの言及が多い。問題提起には事欠かないはずだが、認識は普及せず、解決はなお遠いと思う。
タグ： 超算数
posted at 20:27:07
#超算数 2018年に積分定数さんが取り上げていたときは私あまり関心がなかった話題だ。
解法暗記の唱道者、杉山奈津子氏が出版社の連載でアメリカの実例にサラッと言及。
エスカランテ先生gentosha-go.com/articles/-/5010
アラン・シェーンフェルド教授gentosha-go.com/articles/-/5011
どちらも解法暗記と全く無関係。
タグ： 超算数
posted at 20:22:30
#数楽 添付画像は大数の法則と中心極限定理が見えているとも思える。X_iがランダムに同じ確率で±1なら、n→∞で、それぞれ適切な意味で
① (X_1 + … + X_n)/n → 0
② (X_1 + … + X_n)/√n → Z
となる。Zは標準正規分布に従う確率変数。
スケーリングの違いを表す分母のnと√nの違いに注目。 pic.twitter.com/Q8TaRAZYHx
タグ： 数楽
posted at 20:07:05
沼Utsubo@Kummerspeck @epsilonminder
RでIntel MKLやOpenBLASごにょったりするくらいならpythonやjulia覚えて書いた方が圧倒的に速い
タグ：
posted at 20:03:07
年始から重回帰分析（FE+clustered SE付き）を全てJuliaでやってみてる。N=100万を13本回して5秒前後。快適。PythonとかRでもこれぐらいの速さだっけ？
タグ：
posted at 19:36:12
非公開
タグ：
posted at xx:xx:xx
ふと検索してみると、北大水泳部の背泳ぎの記録、2019年に抜かれてたことに気づいた。shirokumakai.org/?page_id=3796 pic.twitter.com/3UfYsfMrmO
タグ：
posted at 19:05:18
はてなブログに投稿しました #はてなブログ, #Julia言語
Juliaで婚活アプリのシステムモデリング - システムとモデリング
otepipi.hatenablog.com/entry/2021/01/...
posted at 18:52:06
めんどくさくて放ってたメモをまとめたやよ: 「○○を軽視している」という言い分と引用の話：ケイト・マンの事例 flipoutcircuits.blogspot.com/2021/01/blog-p...
タグ：
posted at 18:50:56
私もマルチスレッドに詳しいわけではないのですが、参考程度にどうぞ。#Julia言語
gist.github.com/Lirimy/650171e... twitter.com/bicycle1885/st...
タグ： Julia言語
posted at 18:38:28
#Julia言語 だとこんな感じ。この程度なら `@ inline` つけなくてもインライン展開される♪ twitter.com/korinVR/status... pic.twitter.com/EI3a413ylL
タグ： Julia言語
posted at 18:01:32
#数楽 ランダムウォークの逆正弦法則に出て来るBeta(1/2, 1/2)分布は中心(1/2, 0)半径1/2の円周上の一様分布のx軸への射影に等しい。
Beta(1/2, 1/2)分布の「分母」のB(1/2, 1/2)=Γ(1/2)²は円周率に等しい。
Beta(1/2, 1/2)はBernoulli分布モデルのJeffreys事前分布としても出て来る。 pic.twitter.com/PLzV3MjMiU
タグ： 数楽
posted at 16:56:53
JapaneseFoxBoy @orematasaburou
c++でいうデータがない場合にnullptrを返すようなメンバ関数をjuliaで実装するとき、愚直にはnothingを返せばいいと思うのだが、そうすると戻り値がUnion{Nothing, T}になってしまって性能が落ちそうな気がするのが、どうしたものか
タグ：
posted at 16:45:22
#Julia言語 #数楽
シンプルランダムウォークの中心極限定理 pic.twitter.com/1f44TUS8TR
posted at 16:28:59
非公開
タグ：
posted at xx:xx:xx
#Julia言語 ランダム行列も簡単
添付画像
①n×nの実ランダム行列の複素固有値の分布
②n×nの実対称ランダム行列の実固有値の分布
ソースコード
gist.github.com/genkuroki/5fa6... pic.twitter.com/CQFOcAfFXc
タグ： Julia言語
posted at 15:59:44
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#Julia言語 3次元のランダムウォークもこんなに簡単。疑似乱数のシード値の指定も入れて、たったの4行。
Juliaが「超高級電卓」と呼ばれるゆえん。
10行程度で相当に色々な結果をプロットして確認できる。 pic.twitter.com/QdwV3v0Wh9
タグ： Julia言語
posted at 15:32:12
#Julia言語 細くて半透明(lw=0.3, alpha=0.3)の曲線を重ねてプロットするのは現代におけるプロットの基本定跡の1つだと思う。 pic.twitter.com/KAjksWHQG5
タグ： Julia言語
posted at 15:19:24
#数楽 この結果をランダムウォークの逆正弦法則と呼びます。
公平なギャンブルを長時間続けると、勝っている時間と負けている時間の割合が50%対50%に近くなることは少なく、非常に偏ってしまう確率が高い。
これはその手の遊びが好きな人にとっての経験則とも一致していると思います。 pic.twitter.com/E8RRZSkga3
タグ： 数楽
posted at 14:54:38
#数楽 #Julia言語 解説
0から出発してランダムに±1の移動を10⁴回繰り返した結果を見て、0以上になっていた時間の割合を計算(ちょうど0なら0.5と数える)。
これを10⁵回繰り返し、0以上になっていた時間の割合の分布をプロット。
その分布はBeta(1/2, 1/2)分布でよく近似されていることが分かる。 pic.twitter.com/tm5PJl4z7h
posted at 14:49:44
ランダムウォーク懐かしい．大学院の授業でmt.hとGnuplotのシステムコール使ってC言語でゴリゴリ書いたわ twitter.com/genkuroki/stat...
タグ：
posted at 14:49:18
#数楽 逆正弦法則に至る道のまとめが以下のリンク先スレッドにある。
ワンステップで±1の移動をするシンプルランダムウォークの逆正弦法則なら、お絵描きをしながら直観的に証明可能です。 twitter.com/genkuroki/stat... pic.twitter.com/2rxZECx7R7
タグ： 数楽
posted at 14:32:24
#数楽 #Julia言語
そう！その逆正弦則、結構面白いと思います。
添付画像
①ランダムウォークを200本プロット。赤線は±2√x
②逆正弦法則の確認をメモリの節約をしながら8行で実行！
cumsum! とか rand! のような ! 付きの函数も使えると楽しいです。
ソースコード↓
gist.github.com/genkuroki/5fa6... twitter.com/physics303/sta... pic.twitter.com/yA6PHftxk6
posted at 14:17:43
「正の領域に居る時は善行，負の領域に居る時は悪行を行っている時間とすると，善悪半々というマジメ人間は最も少く，悪い奴ほどはびこり，あるいは善人ほど多いということなのである．これは真実かもしれない」p148
タグ：
posted at 13:15:20
この話は，小針先生の「確率・統計入門」の第6章「酔歩」に証明を含め詳しく議論されてる．この章は独立してるので，確率の基本的なことを知っていれば読めると思う．
古い本だけど，とても読みやすい．
独特の言い回しで著者の息遣いが聞こえてくるすごい良い本だと思う．
amzn.to/3ofQeK3
タグ：
posted at 13:15:19
また，どんなに長い時間観測しても，「１度も原点に戻らない」確率と「１度しか原点に戻らない」確率が最も大きく，m回復帰する確率は，mが大きいほど小さいことが示せる．
原点まわりをうろうろするわけじゃないんですよね．
タグ：
posted at 13:15:19
実は，単位時間あたりに原点より東に居る時間をx，西に居る時間を1-xとすると，この分布は
p(x) = (π√(x(1-x)))^(-1)
になって，東か西かのどちらか片方ばかりに多くいる確率の方が大きい． pic.twitter.com/xLIxJbDkbr
タグ：
posted at 13:15:19
酔っ払いをJuliaで数人召喚して試してみた．私はもっと小刻みに原点付近でうろちょろすると思ってたから，この結果は意外． twitter.com/physics303/sta... pic.twitter.com/wywNWEGbnH
タグ：
posted at 13:15:18
この本はいろんな疑念がすっきり解消するし、問題点が理解できるし、なによりシン・ゴジラみたいでめちゃくちゃ面白いのでおすすめ> 理論疫学者・西浦博の挑戦 新型コロナからいのちを守れ！ | 西浦博, 川端裕人 | ノンフィクション | Kindleストア | Amazon www.amazon.co.jp/gp/product/B08...
タグ：
posted at 12:18:57
Our Milky Way Galaxy_ How Big is Space
(Credit: NASA JPL/Caltech) pic.twitter.com/sgJ4IOen7a
タグ：
posted at 11:19:18
@kyouikukadai @modifyso1 @warabekko 受験票を忘れた場合は係のものに早めに申し出て下さい。受験本部のほうで対応させて頂きます。忘れても受験できなくなるということはございませんので、ご安心下さい。
タグ：
posted at 10:39:42
#統計 Mayoさん曰く
* Royallの3つの問いの「戒律」によれば
①何を信じるか→ベイジアンの事後分布
②どう行動するか→Neyman-Pearsonの方法による長期的なパフォーマンスの良さ
③証拠の比較→尤度比較法
となるが、あなたはこれら全部を拒否したいかもしれません。私が拒否しているように。 pic.twitter.com/qDZ3nMVSLW
タグ： 統計
posted at 10:37:33
#統計 統計学の哲学者のMayoさんによれば
* Royallの3つの問いに、びっくりするほど多くの(しかも異なる信念を持つ)統計家達・哲学者達が賛同している。(大変なことになっている、とてもひどいことになっているというニュアンス)
続く
errorstatistics.com/2014/10/10/bre... twitter.com/genkuroki/stat... pic.twitter.com/7LdK5C0OA5
タグ： 統計
posted at 10:37:31
ごまふあざらし(GomahuAzaras @MathSorcerer
ML ops に関する求人出てますね．
#Julia言語
discourse.julialang.org/t/devops-and-m...
タグ： Julia言語
posted at 10:25:12
#Julia言語
A = []
for i in 1N
push!(A, randn())
end
とするとAnyの配列Aができてしまう。
Anyの配列に格納されたサンプルによるモンテカルロ積分はFloat64の配列の場合と比較すると180倍くらい遅くなった！
A = [] → push!(A, ～) の繰り返しは要注意。
gist.github.com/genkuroki/7cd0... pic.twitter.com/BGQ7coCN36
タグ： Julia言語
posted at 09:41:33
息子が漢字テストでほぼ0点だったことが判明して、休み明けのテストでどこを直せばよいかわからないのでおしえてください pic.twitter.com/dMBV5PnCCx
タグ：
posted at 09:38:34
Masahiro Ono 小野 昌弘 @masahirono
人類にとって幸運なことは、現在開発されているワクチンは微調節が比較的容易にできることである。ワクチン耐性の変異株が出現すれば、それにあわせて、ワクチンの中にいれる遺伝子配列を変異株にあわせてやることで、その変異株にも効果があるワクチンに変えることができる
news.yahoo.co.jp/byline/onomasa... twitter.com/masahirono/sta...
タグ：
posted at 08:43:18
Masahiro Ono 小野 昌弘 @masahirono
ウイルスは免疫のはたらきをすりぬける変異を獲得しうる（免疫回避）。英・南アの変異株は、免疫不全患者が回復者血清等で治療される中で免疫回避の変異を獲得した可能性があり、これを試験管内で再現した実験で意義深い。３箇所の変異で複数の抗体に耐性は、変異株の実験的性質検証の喫緊の必要を示す twitter.com/masahirono/sta...
タグ：
posted at 08:23:24
Julia言語で入門するプログラミング(その5) - SE教育パパむううみんのブログ
muuuminsan.hatenablog.com/entry/2021/01/...
タグ：
posted at 08:14:04
Masahiro Ono 小野 昌弘 @masahirono
査読前論文であるが、回復者血清の存在下でSARS-CoV-2を培養すると、ウイルスは変異を徐々に獲得、80日目までに３箇所の変異(１箇所は11aa挿入)を含むようになり、その回復者血清に耐性になった（中和できなくなった）。同時に複数の回復者血清・モノクロ抗体にも耐性に
www.biorxiv.org/content/10.110...
タグ：
posted at 08:01:25
非公開
タグ：
posted at xx:xx:xx
ハジキのテントウムシを書かないとバツ、はいくら何でもないだろうな、いや、長方形を横×縦でバツの例もあるから、あっても不思議はない。いやいくら教師が馬鹿でもハジキがないとバツなんてことはないだろう。でも、予想を超えた馬鹿教師は存在するから・・・
で、やっぱりあった #超算数 twitter.com/takasoroban/st...
タグ： 超算数
posted at 07:56:46
検定いらずのABテスト：ポアソン分布とベータ分布によるサンプルサイズ設計 - ジョンとヨーコのイマジン日記 selfboast.hatenablog.jp/entry/2021/01/...
タグ：
posted at 06:27:29
受験生の方へ。入試当日受験票を忘れてもそんなに大変ではありませんので、万一、そうなっても慌てずに係の人に申し出てください。具体的な対応については「受験上の注意」に書かれています。例えば、令和3年度共通テストなら、§6(1)（9ページ）に書かれている指示に従ってください。 twitter.com/kyouikukadai/s...
タグ：
posted at 04:16:44
富谷(助教);監修 シン仮面ライダー @TomiyaAkio
結局、Julia が速いのって
極限までの型推論→型に応じた最適化→最適化されたバイトコード→最適化されたネイティブコード
ってことなんかな
タグ：
posted at 03:29:10
World and Science @WorldAndScience
Awesome kitchen gadget! Cut Resistant Gloves: bit.ly/2SknuQr -AmazonAfflink
タグ：
posted at 02:40:02
非公開
タグ：
posted at xx:xx:xx
#統計 問題解答例: 2×2の分割表
a b
c d
での仮説「オッズ比=ω」のχ²検定(ω=1の場合が独立性検定):
A=a-δ, B=b+δ, C=c+δ, D=d-δとおき、
AD/(BC)=ω, -min(b,c)<δ<min
でδを決めて、近似的に自由度1のχ²分布に従う統計量
X² = δ²(1/A+1/B+1/C+1/D)
を使う。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat... pic.twitter.com/mJkKUbowrt
タグ： 統計
posted at 02:09:40
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@sumikkophysik plot(x,y,z)で出るみたいですよwww.cas.cmc.osaka-u.ac.jp/~paoon/misc/ju...
タグ：
posted at 01:07:40
水波ユキノ@skeb依頼受付中だいたい絵 @minami_suiha
最近数学だとR言語ばっかり触ってるからちょい、Juliaもおさわり。
Python（常用）
R言語（数学のときだけ常用）
Julia（勉強必須）
どーすんだ・・・ pic.twitter.com/63n5D2beQf
タグ：
posted at 00:59:10
@uKi2wQXyG7rx3gL @taifu21 @machidasensei 「うちの料理がまずかったら、私じゃなく食材納品業者に文句言ってください」と言ったら
「その業者から仕入れて料理を提供しているのはあなたでしょう。」と言われるでしょうね
twitter.com/machidasensei/...
タグ：
posted at 00:35:26
