黒木玄 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月07日(木)
長いキャッシュフロー計算をする時、for文が何重かになっている上に中身が複雑になっていて、やたらと計算遅いことがそこそこある気がします。僕は、最近まで自堕落に碌に最適化されていないpythonコードを使っていましたが、juliaで書き直したら、お茶飲んでる時間がなくなりました。
タグ：
posted at 00:11:46
A pool is a good context to test some fluid dynamics concepts, like this water wall, interestingly visualized in slow motion [source of the clip, TikTok's @don.hesel: buff.ly/3nhxcBB] pic.twitter.com/xSyAFI72Sq
タグ：
posted at 00:40:06
非公開
タグ：
posted at xx:xx:xx
メタプログラミングすごい！、とは思うんだけど、自分の乏しいプログラミング経験では、数値計算においてどう役に立つかはよく分からない…。
マクロを自作したり、シンボルを駆使する状況が思い浮かばないのである…。
タグ：
posted at 01:47:47
@mayataka_ #Julia言語 行列Aの対角成分のみを置き換えたい場合には
A[diagind(A)] .= 0 # 対角成分を0に
のように、diagind(A)を使えます。続く
gist.github.com/genkuroki/f65f... pic.twitter.com/C8D9iVo6Q3
タグ： Julia言語
posted at 04:10:56
@mayataka_ #Julia言語 公式ドキュメントで探すのは大変なのですが、添付画像のように適当にヘルプを調べているうちに「もしかしてこれでは？」のようなキーワードが見つかる場合が結構多いと思います。 pic.twitter.com/UMLsRwkDuA
タグ： Julia言語
posted at 04:12:53
#Julia言語
メタプログラミング = コードの自動生成
は数値計算高速化の定跡の1つになっているようです。
Julia界の有名人のSteven G. Johnsonさんはまさにその専門家で、高速フーリエ変換のライブラリFFTWもコードの自動生成を使っています。
この動画がおすすめ！↓
youtu.be/mSgXWpvQEHE?t=85 twitter.com/physics303/sta...
タグ： Julia言語
posted at 04:45:46
#Julia言語 繰り返し述べていることですが、私が「Juliaはよい」と確信したのはSGJさんのMITでの講義の宿題の模範解答
nbviewer.jupyter.org/github/steveng...
を見たときです。
そこでは指数積分函数E₁(z)をJuliaのメタプログラミングで書いて、あるFortranライブラリの5～6倍(❗️)の速さを実現しています。
タグ： Julia言語
posted at 04:45:46
#Julia言語 上で引用したMITの講義の宿題の答えは現在のJulia v1.x では動きません。私がそのコア部分を現在のJuliaで動くようにしたものが
nbviewer.jupyter.org/gist/genkuroki...
にあります。 pic.twitter.com/GPuwVMUbXi
タグ： Julia言語
posted at 04:45:48
各節のタイトルが【10,000回 電源をぶちぎってみる】などになりまくっていて笑える。娯楽の1つとして楽しんで読める。
こういう話は楽しい。 twitter.com/ksksue/status/...
タグ：
posted at 05:00:31
【産技研なんかで1時間700円くらいで静電気ガンを貸してくれる】←こういう話も楽しい
Raspberry Piは本当に壊れやすいのか｜su zenn.dev/su_/articles/9... #zenn
タグ： zenn
posted at 05:06:32
#Julia言語 初期化されていない配列を
A = Array{Float64, 3}(undef, 4, 3, 2)
のように作れますが、よりシンプルに
B = similar([0.0], 4, 3, 2)
で作れます。
similar(X, n₁, n₂, ...) は「Xと似たものをサイズ(n₁, n₂, ...)で作る」という意味です。
gist.github.com/genkuroki/0e29... pic.twitter.com/oOzwN4Yz7w
タグ： Julia言語
posted at 05:54:57
#Julia言語
個人的には、multiple dispatchを使う場合以外には、できるだけ、Float64とか書きたくない。
similar(a)
zero(x), one(x), oneunit(x)
eltype(a), typeof(x)
などで済ませたい。
docs.julialang.org/en/v1/base/num...
タグ： Julia言語
posted at 06:12:10
@kuchibirusensei @tarookung @May_Roma 失礼します。
twitter.com/genkuroki/stat...
黒木玄さんのツイートの画像にあるレシートは会計学に反しているのでしょうか？ pic.twitter.com/61DT2pMilU
タグ：
posted at 06:47:32
"紙書籍は84万部、電子書籍は13万部超、Audibleは3万超"
＊ ＊ ＊
『ファクトフルネス』が100万部突破！中学・高校の先生向けガイドを配布 prtimes.jp/main/html/rd/p... @PRTIMES_JPより
タグ：
posted at 07:04:58
@sekibunnteisuu @kuchibirusensei @tarookung @May_Roma 失礼します。
国税庁の説明では「単価×数量」と「数量×単価」の順序が混在しています。国税庁は会計学に反しているのでしょうか？
🤣🤣🤣🤣🤣🤣
www.nta.go.jp/law/joho-zeika... twitter.com/temmusu_n/stat... pic.twitter.com/G6Uov0EZc4
タグ：
posted at 07:10:32
イギリスでの感染爆発は変異ウイルス「だから」ではないでしょうね👶なぜイギリスの科学者はそこを強調する人が少ないのか。なにか圧力があるの？それとも本気で変異の「せい」だと信じてるの？…解せぬ
タグ：
posted at 08:05:37
非公開
タグ：
posted at xx:xx:xx
@to_dried_bonito 「感染力」が「1.7倍」という科学的根拠はどれなんでしょう👶科学者の「一部」のグループが言いつのっているのを無批判に信じる以上に検討していない、ということはないですか？
タグ：
posted at 08:50:47
@to_dried_bonito ばぶは結構引いて見ていますし聞いています👶
イギリスの複数の学者が「感染力」が上がっていると示唆していますが、いずれも統計・モデリングによるもので、メカニズムや実証的検証はほぼないと見ています。報道は「科学者」を信用した権威主義に引っ張られている、とこのケースは見ています。
タグ：
posted at 08:55:45
@to_dried_bonito これについてはここにリンクをおいています👶どう読み解くかが大事ですが、読み解けないから科学者を信じる、も危ういというよいケーススタディになると思って注目しています。
covid-19.babubabulog.com/%e3%82%a6%e3%8...
タグ：
posted at 09:00:10
まあ、そういう意味で言うと10年と言ってること変わりませんね。
「皆さんの言いたいことはよくわかるし、自分も同じ気持ちだ。だが、皆さんのいう《銀の弾丸》は解決策として現実的に機能しない。皆さんの望む結果を得るにはこちらの現実的な策を使ったほうがいい」という構造は今も昔も変わらない。
タグ：
posted at 11:26:00
今日、ファイザー社の新型コロナワクチンを接種してきました。気になる症状はなく、一日中いつも通り働いてました。十分な数の臨床研究データも公表されてたので、安心して接種できました
驚異的な速さの開発の影には、研究者/医療者/製薬会社/政府の方々の凄まじい努力があります。心から感謝します
タグ：
posted at 11:41:58
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Common Lispの=、eq、eql、equalの使い分けみたいなものが（ポリモーフィズムのおかげで）暗黙のうちに行われている感じ？
#Julia言語 にも==、===、isequal()がある。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 13:00:16
モグラやケラなど土を潜る生き物は頑丈な手で土を掘り進めますし、ミミズのようにまっすぐに伸びた紐状の生き物は体を収縮して頭から少しずつ潜ります。
でも、カブトムシの幼虫には頑丈な手もなければ、真っすぐなることもできません。
どうやって潜るのだろうか、、とずっと不思議でした。
続く pic.twitter.com/mpC7apBAFI
タグ：
posted at 13:03:33
不思議でしたので、観察しました。
すると、頭から土に潜り、ショベルカーのように土を掘り、潜り、ある程度土に潜ったら、体をねじって反転して、回転してから頭を下に向けて、また同じように掘って潜るようでした。
体をイモムシみたいに伸ばせないカブトムシの幼虫の潜り方がやっと理解できました pic.twitter.com/i1YWkLKyr7
タグ：
posted at 13:06:36
非公開
タグ：
posted at xx:xx:xx
末石さん @yaonacs が年末に呟いていた件．Bruce E. Hansen, A Modern Gauss-Markov Theorem. ガウスマルコフの定理で「線形」を外せる，つまりOLSは非線形な不偏推定量を含めたクラスでもベスト．２次のモーメント条件のみ必要．かっこいい結果です．www.ssc.wisc.edu/~bhansen/paper... pic.twitter.com/EI3SBVwkgx
タグ：
posted at 14:12:35
証明は自分には難しそう．based on a trick used by Newey (1990) in his development of asymptotic semi-parametric efficiency bounds. とのこと．
タグ：
posted at 14:14:00
Hiroyuki Chihara @HiroyukiChihara
I computed the multiresolution approximation of a relatively uncomfortable audio signal given by
x(t)=sin(2t(t-3)(t-6)(t-9)(t-12)).
These pictures might be helpful to understand the multiresolution decomposition.
github.com/JuliaDSP/Wavel...
#julialang twitter.com/HiroyukiChihar... pic.twitter.com/Mv7K9ZAgaI
タグ： julialang
posted at 14:22:56
非公開
タグ：
posted at xx:xx:xx
岩田先生は決して抗菌「薬」だけの専門家ではなく、疑いようなく日本を代表とする感染症専門家の一人なので山本氏がそこを矮小化するのはちょっとお門違いだと思います。先生の問題点はそこではなくてｺﾞﾆｮｺﾞﾆｮ…。 twitter.com/ichiro_leadoff...
タグ：
posted at 15:15:06
外科志望の医学生が白い巨塔を読んで傲慢でない財前五郎になりたいのと思うのと同じように、感染症科志望の医学生が名著「抗菌薬の考え方、使い方」を読んで(ここで文は途切れている)
タグ：
posted at 15:23:57
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
実際にはグルーミング（体表のお掃除）ですね。口脚で腹部の遊泳肢をグルーミングするにはこういう体制になります。 twitter.com/uxotto/status/...
タグ：
posted at 16:50:26
Juliaでテンソル．これ，もっとかっこよい方法ないん？（2×2行列3つから2×2×2テンソルを作ってる）
ttmとかうまく使おうと思ったけどよくわからなくなってしまった． pic.twitter.com/CXN70WyFIM
タグ：
posted at 16:57:48
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Juliaで書かれた第一原理計算のコードを発見した。まだ開発段階みたいなので、開発者の方にはぜひ頑張って欲しい。
github.com/JuliaMolSim/DF...
タグ：
posted at 17:34:55
非公開
タグ：
posted at xx:xx:xx
コロナワクチン 米で約190万人中21人に激しいアレルギー反応 | NHKニュース www3.nhk.or.jp/news/html/2021...
アナフィラキシーを重篤な副反応考えても、21人/190万人×100＝0.001 ％、効果が確実ならめちゃくちゃ安全性の高いワクチンじゃん。
タグ：
posted at 18:13:38
#統計 興味深い話。
しかし、いつも理解できないのは unbiased という強烈に強い仮定の重要性。
実際に数値実験してみると、残差が正規分布でない有限の分散を持つ分布q(e|θ₀)に従う場合の中には、最小二乗法の誤差が残差がq(e|θ)に従うとするモデルよりも非常に大きくなる場合があります。 twitter.com/umaruyama/stat...
タグ： 統計
posted at 18:16:34
#統計 分散共分散行列が単位行列で期待値μ=(μ_1,…,μ_d)が未知のd≥3次元の正規分布のサイズn=1のサンプルX=(X_1,…,X_d)からμを推定せよという苦しい状況で、μの推定量としてX(=標本平均!)を採用することが平均二乗誤差の意味で最良では__ない__、というのがStein推定の話。
nbviewer.jupyter.org/gist/genkuroki...
タグ： 統計
posted at 18:30:14
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#Julia言語 Juliaの配列はbegin始まりインデックスなので、0始まりインデックスのa[i]をJuliaではa[begin+i]に置き換えて行けばミスを減らせると思います。続く twitter.com/kivantium/stat...
タグ： Julia言語
posted at 19:02:37
#Julia言語 もしくはOffsetArrays.jlの利用も検討してみる。Juliaは配列のインデックスの始まりを任意に動かせます。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 19:02:38
#Julia言語 関連
インデックスの始まりを0や1に固定されること自体が不自由で不便だと思う。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 19:02:39
Zenn, QiitaにJuliaの記事が投稿されると自動でツイートするbotです！
気が向いたらその他のサイトにも対応する予定です。
不具合報告などあれば: @abap34 pic.twitter.com/8xhjS531Dp
タグ：
posted at 19:02:48
非公開
タグ：
posted at xx:xx:xx
Raspberry Pi 4とYamaha THR-10を接続してみたら、何事もなくつながってしまった。REAPERからも普通にALSAで使えるし、これは便利だな。
タグ：
posted at 21:47:46
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Ninja DAO | CryptoNi @CryptoNlnjaNFT
日本語での「トランプ」＋ANTIFA で詮索すると、予めお付き合いをお断りしたい方々のリストがほぼ自動的にできる事がわかった。
タグ：
posted at 21:54:31
Raspberry PiでREAPERもJuliaもRも動くし、TeXとEmacsもあるし、PureDataもGIMPもInkscapeもあるので、これで良いのでは。ChromeからTwitterできるし。
タグ：
posted at 22:01:38
岩田健太郎がとくダネ出演とかw
逃げろ！伊藤隼也！！
医療ジャーナリストの化けの皮が剥がされるからっ！！！ twitter.com/TOTB1984/statu...
タグ：
posted at 23:00:38
しかし岩田健太郎も伊藤隼也も反自民党だから、イデオロギーで共闘する可能性もあるねw
伊藤隼也
『政府が勧めるインフルエンザワクチン効果なし！』
岩田健太郎
『おおっ、そうだな！政府はインフルエンザ製造機だっ！！動画Upしよ！』
タグ：
posted at 23:05:38
>西浦さんが検査抑制論
いつでも誰でも何度でも、根拠なくレッテルを貼ることだけは達者
>少なくとも検査拡充を《積極的に》主張していない
「検査拡充を主張していても『検査抑制論』」ということですね。はい、終いにはこう来ると思っていました。 twitter.com/masakioshikawa...
タグ：
posted at 23:07:28
伊藤隼也さんという反医療の人(及びそういう人を使っているマスコミ)を名指しでボコって下さるという噂を聞きました。
期待しています！ twitter.com/georgebest1969...
タグ：
posted at 23:11:50
非公開
タグ：
posted at xx:xx:xx
@aquablau @togetter_jp 専門分野では業績を残している方ですが、なぜわざわざよく知らない分野で他者を断罪したり不当に貶めたりするのでしょうか。
タグ：
posted at 23:16:45
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
私は確信していましたよ。
「検査拡充を主張していたとしても『抑制論』」
と彼らが絶対に言い出すと。 twitter.com/tsukuru_ouu/st...
タグ：
posted at 23:24:37
@sabamisolover あと、#超算数 とも絡んできますが、2.305 これを、「にいてんさんぜろご」と読むのは誤り、「にいてんさんれいご」が正しい、というのもあります。
しかしそれ言ったら、「24時間」を「にじゅうよじかん」というのもおかしいという話にもなります。
タグ： 超算数
posted at 23:33:10
@nekomath271828 @taketo1024 わざと改良の余地のある状態でプログラムを公表すると、「1週間待っていろ。本物の〇〇を見せてやる」な人が勝手に改良してくれます。
タグ：
posted at 23:48:06
