黒木玄 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年03月16日(月)
ごまふあざらし(GomahuAzaras @MathSorcerer
んー．plot 系のライブラリ，
大量の引数オプション要請されること多いからキーワードで順序変えられたほうが絶対いい．
タグ：
posted at 00:07:00
Windows で Linux バイナリを実行するための技術である WSL 2 (Windows Subsystem for Linux version 2) の
更新の話💡奇態
🎀 / マイクロソフト、WSL 2でファイルシステムI/O性能最大20倍高速に (マイナビニュース) #NewsPicks npx.me/F18w/avP5?from...
タグ： NewsPicks
posted at 00:13:50
ちゃんちゃん | JSとKotlinした @fashioncrazy66
@chomado WSL1に比べてDocker周りが
かなり強化されてそうで嬉しい！
タグ：
posted at 00:19:09
ごまふあざらし(GomahuAzaras @MathSorcerer
それは，通常のLinuxで動かした時と比べてどうなるんですか？って誰かメンションして？
タグ：
posted at 00:25:55
ごまふあざらし(GomahuAzaras @MathSorcerer
ホォ？ HCubature.jl とどう違うんや？
github.com/stevengj/Cubat...
タグ：
posted at 00:34:16
#Julia言語
Cubature.jl はCで書いた数値積分のライブラリをJuliaで使えるようにしたもの。
HCubature.jl は pure Julia で書かれた数値積分パッケージ。
multiple dispatch が効くので適用可能範囲が広い。
で、非常にびっくりなことに、C版よりもpure Julia版の方が計算が速い！ twitter.com/MathSorcerer/s...
タグ： Julia言語
posted at 00:41:02
ごまふあざらし(GomahuAzaras @MathSorcerer
@genkuroki 多次元での計算だとHCubature.jlが有利ということがわかりました．しゅごい．
gist.github.com/terasakisatosh...
タグ：
posted at 01:46:22
ごまふあざらし(GomahuAzaras @MathSorcerer
扱うターゲット次第だけれど，Pythonみたいな設計をJuliaで期待するとすごく苦しくなる．
＞RT
タグ：
posted at 01:52:01
ごまふあざらし(GomahuAzaras @MathSorcerer
たった３行で雛形ができちゃうできゅ〜
julia> using Franklin
julia> newsite("mySite", template="pure-sm")
julia> serve() twitter.com/hs_heddy/statu...
タグ：
posted at 01:54:17
ごまふあざらし(GomahuAzaras @MathSorcerer
... は "splat" operator と呼ばれるものですね．
REPLで ? を入力した後 ... を押すとわかる． twitter.com/misaki70533995... pic.twitter.com/NJA77AgW7J
タグ：
posted at 01:58:58
英首相の「降伏」演説と集団免疫にたよる英国コロナウイルス政策(小野昌弘) - Y!ニュース news.yahoo.co.jp/byline/onomasa...
タグ：
posted at 02:27:28
ごまふあざらし(GomahuAzaras @MathSorcerer
完全にこれにハマってる・・・．
x=3
Atmark m x
で An Integer のプリント出させられないのか． pic.twitter.com/Pw8oRzx2nW
タグ：
posted at 03:11:57
Takashi Watari M.D, @wataritari1
これ勘違いで、間違いです。イタリア北部に滞在と言う『検査前確率が異常に高い帰国者』に適切に行ったからです。EBM＋診断学の基本なんです。検査は必要な時に確率を高めて適切に行う、コレが大事。
私の医師免許は舛添さんが厚生労働大臣でした。それでもこの科学的思考はむずかいしいのでしょうか。 twitter.com/masuzoeyoichi/...
タグ：
posted at 06:56:44
中央集権的な制度はサムシング・グレートにやられ、地方分権的な制度はゲーム脳にやられ、いずれにしても日本会議の功績が大きいですね。気息奄々たる21世紀の初頭です。
タグ：
posted at 07:51:36
JunoチームでJuliaConのproposalを出した :)
通りますように．そしてonsiteで開催できますように ...
pretalx.com/juliacon2020/t...
タグ：
posted at 08:44:58
Takashi Watari M.D, @wataritari1
コメント感謝です。厚生労働省より:こちらに掲載。我々が徹底的に疑って検査を選別しても非流行地域ではこのような結果でした。『検査は必要な時に検査前確率を高めて適切に行うべき』です。TVの意見とは異なるかもしれません。我々は手洗い厳守、人混みを避けて、勉強する、コレしかありません。 pic.twitter.com/DX9u5E7onV
タグ：
posted at 09:44:03
秋田県の知事らの会見ではマスクを着用していない。
・感染の疑いのある人が（感染を）広げるのを防ぐための着用は必要
・明確に感染を予防できるという医学的知見は確立していない
・マスク不足、効果の薄い対応は行うべきではない
・ちょいちょいあごに下げるのは最もやってはいけない
ド正論。
タグ：
posted at 10:02:07
このコメントを出した諸冨伸夫健康福祉部長は、2019年度に部長昇格した元健康医療技監で、医師免許を持つとか。写真見つけました、お若いように見受けられます。
www.sakigake.jp/news/article/2...
タグ：
posted at 10:15:39
#Julia言語 統計がらみの数値計算でlog(1+x), log(exp(x)+exp(y)), log Σ_i exp(X[i])という計算がよく出て来ます。それらを数値計算するときには専用の函数
log1p(x)
StatsFuns.logaddexp(x, y)
StatsFuns.logsumexp(X)
を使うべき。
github.com/JuliaLang/juli...
github.com/JuliaStats/Sta...
タグ： Julia言語
posted at 11:08:40
#Julia言語 log1p, StatsFuns.logaddexp, StatsFuns.logsumexp の使用例。exp(x)は簡単に0.0やInfになってしまうので要注意。
logintexp(f, a, b) = log ∫_a^b f(x) dx
もあると便利そうですね。階層ベイズでのWAIC, LOOCV, などの計算で出て来る。 pic.twitter.com/MQUIwiCvwO
タグ： Julia言語
posted at 11:36:27
#Julia言語 あと
x log y
の数値計算では 0 log 0 = 0 と計算するようにしておかないとはまる。例えばベルヌーイ分布モデルの対数尤度函数
log(pᵏ(1-p)ⁿ⁻ᵏ) = k log p + (n-k)log(1-p)
では、p = k/n, k = 0 のとき 0 log 0 が出て来て、NaNになる。
StatsFuns.xlogy が使える。 pic.twitter.com/DgThg65Pem
タグ： Julia言語
posted at 11:59:57
#Julia言語 あと、階乗 n! の計算も要注意。n! をナマで扱うのではなく
log n! = log Γ(n+1)
を扱わないと簡単にオーバーフローしてしまう。ただし
log(factorial(n))
は n=21 でオーバーフローする。
SpecialFunctions.logfactorial(n)
や
SpecialFunctions.loggamma(n+1)
を使うとよい。 pic.twitter.com/f8HtWeIh4f
タグ： Julia言語
posted at 12:11:46
私は「エビデンス」という言い方は昔から嫌い。
「証拠」でいいんじゃね？
統計的証拠の意味なら「統計的証拠」とか「スタティスティカル・エビデンス」とか言ってほしーでしゅ。
統計的証拠を使う判断は常にギャンブルになる。
そしてギャンブルで確率やデータを無視する奴は負け組確定。
タグ：
posted at 12:37:45
JuliaのFloat64とBigFloat使った計算の比較、久しぶりだったのでそれぞれ別の関数名で実装してたんだけど、途中で多重ディスパッチ案件やんけとなった
タグ：
posted at 13:05:25
#Julia言語 これ、尤度函数
L(w) = p(X_1|w)…p(X_n|w)
のプロットで要注意。nが大きいと容易に0.0やInfになる。
うまく定数倍を繰り込んで計算して(対数をとって最大値を引いてからexpするとかして)プロットする必要がある。
ベイズ分配函数も大変。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 13:47:13
統計力学の講義ノート、図を入れたバージョンを公開します。不適切な書き方をしている部分やケアレスミスも相当直しましたが、まだ間違いがあるかもしれません。完成間近のベータ版です。お暇なかたは間違いを見つけてください。
www.cp.cmc.osaka-u.ac.jp/~kikuchi/kougi...
タグ：
posted at 14:27:07
柳 時熏【囲碁プロ】 【YouTuber @ryu_shikun
磁石大盤を製作しました。
家を建てた時に余った壁紙に線を引いて、マグネットボードに貼りました。
解説用なので、線の番号と星印は必要かな？
白黒石はとりあえず100円均一で購入したマグネットで、代用しました。
現在ブログの準備をしています。
お楽しみに😊✨
#囲碁 #大盤 pic.twitter.com/xECfnZwQHa
posted at 14:50:27
#Julia言語 でPlotsが使えず，エラーメッセージを読むとFixedPointNumbersパッケージの読み込みでEOFエラーが出ているようです．何かご存じの方がおられればご教示いただきたいです… pic.twitter.com/cTWlZTPwBP
タグ： Julia言語
posted at 15:46:34
#Julia言語 ~/.julia/compiled を削除することで解決しました！ twitter.com/dolphin7473/st...
タグ： Julia言語
posted at 16:24:23
#統計 え？褒めても大丈夫な本なのですか？
公開されている『瀕死本』
www.asakura.co.jp/books/isbn/978...
のプロローグや
www.jstage.jst.go.jp/article/sjpr/6...
p 値を使って学術論文を書くのは止めよう
豊田秀樹
を見るだけでもトンデモないことが書いてある。
詳しくはこのツイートが繋がっているスレッドを参照。 twitter.com/mutopsy/status... pic.twitter.com/ln1JWR2Upw
タグ： 統計
posted at 16:30:25
@MathSorcerer #Julia言語 ただし、
using HCubature
hcubature(x -> exp(704 - (x[1]^2 + x[2]^2)), [-10, -10], [10, 10])
とすると、計算が止まらなくなる問題があって、誰か何とかしてほしーです。
maxevals=100 とかにしても解決しない問題なので、たぶんバグだと思う。おそらく、Inf がらみの問題。
タグ： Julia言語
posted at 16:39:02
#Julia言語 logintexpの私による実装例
nbviewer.jupyter.org/gist/genkuroki...
logintexp.jpynb pic.twitter.com/jIPdnU7EsT
タグ： Julia言語
posted at 16:43:17
#Julia言語
log ∫_a^b exp(f(x)) dx
を max_{x∈[a,b]} f(x) = f(x₀) と書くとき
f(x₀) + log ∫_a^b exp(f(x) - f(x₀)) dx
で計算することは、統計学の文脈では最尤法やMAP法の解 (x₀, f(x₀)) を途中で求めていることにもなりますね。
gist.github.com/genkuroki/9024...
logintexp.ipynb
タグ： Julia言語
posted at 16:55:01
やっとPlots使えた(*´-`)
github.com/jheinen/GR.jl/... のGRは gr-framework.org/building.html でcloneしてるgrのことだったんですね(￣ー￣) GRというフォルダも（あちこちで）勝手にできてるからCMakeLists.txtが入ってるのがないか探したり何かしたら入るんじゃないかってググったりしてた #julialang pic.twitter.com/icq42OaktF
タグ： julialang
posted at 17:33:49
（GIFアニメは昔JuliaBoxで作ったやつ（gist.github.com/brv00/118ed8ec...）をダウンロードして（少し書き換えて）実行して作った）
タグ：
posted at 17:40:18
@MathSorcerer #Julia言語 訂正: maxevals を設定すれば止まります。
hcubature(x -> exp(704 - (x[1]^2 + x[2]^2)), [-10, -10], [10, 10]; maxevals=10^7)
で2秒弱。maxevalsの設定はこの辺が安全ということかも。理論値の exp(704)π が正確に計算できています(対数を取って確認している)。 pic.twitter.com/iRGGXkOXHq
タグ： Julia言語
posted at 17:40:45
非公開
タグ：
posted at xx:xx:xx
#Julia言語 関連情報
HCubature.hcubature で数値積分をするときには、計算が全然止まらなくなる場合があるので、maxevals=10^7 くらいにしておけば、数秒程度で計算が止まる(誤差は大きくなるかもしれない)ので、ちょっと安心かも。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 17:45:36
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@mutopsy それはわかります(→添付画像の内容)。
しかし、その本の内容に問題があることに気付いていなかったり(無知無能)、問題に気付いているのにそのことに何も触れずに褒めるだけで済ます(通常悪い行為とみなされる)のはまずいと思いました。 pic.twitter.com/8MQ8O2ddr5
タグ：
posted at 18:10:30
@genkuroki なるほど，意図を理解いたしました。僕自身は実際にこの本から学べることがたくさんあると確信しているので素直に好意的なコメントを付したのですが，誤解を避けるために注を付けたほうが良かったかもしれませんね。
タグ：
posted at 18:20:49
#統計 『瀕死本』を沢山売ることは社会的には負の貢献になると思う。
豊田氏の意味での「仮説が正しい確率」の中には片側検定などのp値と近似的に一致することが分かっているものが多い。
豊田さんは「p値は仮説が正しい確率の近似値である」と言っているとみなされます。
twitter.com/rikoushonotana...
タグ： 統計
posted at 18:51:21
#統計 『瀕死本』のp.86には「仮説が正しい確率」については
【何回計算しても多重比較のような補正が必要ありません】
と書いてあって、非常にまずいと思いました。
「仮説が正しい確率」が何らかの意味でp値と近似的に一致する場合があるのでそれは本当にまずい。
twitter.com/genkuroki/stat...
タグ： 統計
posted at 18:54:50
新型コロナを理由とする学生への内定取消について、「学生が悲鳴」「これをチャンスとみる企業も」「国は支援を」など書く記事などもあるけど、そもそも内定取消は通常の労働者の解雇と同等の要件がない限り認められない労働基準法違反で無効なものということを明確にすべき。
www.kigyou-houmu.com/labor-4/ pic.twitter.com/WZ3ANQswQC
タグ：
posted at 19:13:42
電子版！他ストアにもあるので検索してみて下さい🙇♂️
BookLive! →
booklive.jp/product/index/...
bookwalker→
bookwalker.jp/de0adc8608-c24...
LINEマンガ→
manga.line.me/book/detail?id...
楽天ブックス→
books.rakuten.co.jp/rk/0aeb1e84d25...
コミックシーモア→ www.cmoa.jp/title/195153/
amazon→www.amazon.co.jp/dp/B085HDHPRQ/... pic.twitter.com/e2MqjUkS9f
タグ：
posted at 19:33:23
コミック版、Amazonは現在待ち時間があるみたいなのですが、こちらも送料無料で！！すぐ届くので…！！
楽天ブックス→books.rakuten.co.jp/rb/16230927/?l...
ヨドバシ.com→ www.yodobashi.com/product/100000... pic.twitter.com/ORAgnEKTwO
タグ：
posted at 19:35:00
@temmusu_n @LimgTW @hgn_no_otaku @miyuki_MathT マッチポンプの可能性は大きいと思ってます。
水道方式も駄目な教え方をしている。
togetter.com/li/1411620
タグ：
posted at 19:36:06
@temmusu_n @LimgTW @hgn_no_otaku @miyuki_MathT twitter.com/sekibunnteisuu... にも書いたけど、この問題で一桁違った値を出してしまうと言うことは、
「小数の計算が分かっていない」じゃなくて、「ナントカ算をすればいい」となってしまっているから。
それまでの教え方を反省すべきだと思う。
タグ：
posted at 19:40:17
非公開
タグ：
posted at xx:xx:xx
@sekibunnteisuu @LimgTW @hgn_no_otaku @miyuki_MathT 理解を謳うのになんでマナー商法になっちゃうんでしょうね。いつも不思議です。
タグ：
posted at 19:56:46
すみません、あのー、立命館の教員のものですが…… webの授業ってどうやったらいいか誰か教えてくれませんか……………… 何したらええんや……… twitter.com/ritsumeikan_pr...
タグ：
posted at 20:06:13
ごまふあざらし(GomahuAzaras @MathSorcerer
あれは思い返すと伝説級の珍事件でしたね． twitter.com/kado_judo0312/...
タグ：
posted at 20:53:18
ごまふあざらし(GomahuAzaras @MathSorcerer
電子版， 購入しました！ twitter.com/DEATHcotori/st...
タグ：
posted at 20:55:12
非公開
タグ：
posted at xx:xx:xx
#Julia言語
そういう用途には最も初心者向きで易しいJuliaをどうして使わないのか？
JuliaとC++やRustの計算速度は同じ程度。
JuliaならmatplotlibでのプロットもPyPlot.jl経由で易しい。
csvファイルに計算結果を保存とかも当然いらない。
Juliaの公式ドキュメントは読み易いです。 twitter.com/phasetrbot/sta...
タグ： Julia言語
posted at 21:08:05
「瀕死本」の序盤を読んでみたが、そもそも正答率50%を帰無仮説にした二項検定と「超能力があれば70%以上は正答するという仮説の正しい確率」なるものを比較して二項検定をこき下ろすのはアンフェアでしょうとか思っていたら、@genkuroki さんがわかりやすいデモをあげていた
twitter.com/genkuroki/stat...
タグ：
posted at 21:21:16
ごまふあざらし(GomahuAzaras @MathSorcerer
「コーディング規約：
startup.jl にローレンツアトラクターを書いてはいけない」
って前も呟いた気がするw
タグ：
posted at 21:41:28
科学としての実質的な効果を議論しないとダメだよね、というのはその通りだと思うのであるが、社会心理学でよく使われる「質問紙」得点の「実質的な効果」とはなんだろうか。
タグ：
posted at 21:48:11
裏を返すと、我々が必要としているのは、実験室環境で得られる指標や質問紙で測られるものと実世界とのつながりなのでしょうな(うつ病得点みたいなものはともかくとして)
タグ：
posted at 21:55:15
非公開
タグ：
posted at xx:xx:xx
ざっと読んだところ「瀕死本」の第一部で言いたいのは
・有意性検定でのみ科学的主張を行うのは色々まずい
・その結果、実質的な効果がないような論文が量産される
・n設計や事前登録にはしばしば困難さがつきまとう
ということで、扇情的な表現が多いことに目を瞑れば言いたいこと自体はわかる
タグ：
posted at 22:12:14
「検査しろ派」の物理学者がいるようですが、僕はブロックしてると思うので知らん。物理学者がみんな「検査しろ派」なわけではないので、そこんとこ誤解ないようお願いします
タグ：
posted at 22:12:56
二部からが本番で、タイトルは「ベイズの定理・仮説の正しい確率」である。最初のページの例が迷惑メールフィルタなので、その時点でだいぶ暗雲が立ち込めているといった趣だ。
タグ：
posted at 22:16:18
ごまふあざらし(GomahuAzaras @MathSorcerer
LLVMのバージョン事情把握し切れてないんだけれど、Juliaと対応するLLVMのバージョンゲキアゲしてない？
タグ：
posted at 22:19:53
@metameta007 2点確認なのですが、10未満の数を数えるのでも1あたりは1として、7なら「1が7つ」ということで、被乗数と乗数の概念が必要なのでしょうか？
2点目ですが、今この段階でこの2点目は、1点が2つということで、被乗数と乗数の概念を必要とした2点のうちの2点目でしょうか？
タグ：
posted at 22:25:06
