黒木玄 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年11月29日(月)
例えば、データの数値に台を持つデルタ分布モデルの尤度は無限大になるが、そのようにして尤度を大きくしても何の意味もない。
モデルを非常に複雑にしてパラメータを大量に入れておけば、最尤法によってモデルのデータへの適合度を高くできるが、オーバーフィッティングを引き起こす。
タグ：
posted at 23:28:55
@takitoshi ubuntuだと、/usr/lib/x86_64-linux-gnu/julia/libopenblas64_.so が実体っぽいので、これを上書きしちゃうなり、LD_LIBRARY_PATHで誤魔化すなりすればいけるかも… まあ、こういう泥縄で動いてもそれは動いたとは言わない気がする笑
タグ：
posted at 23:26:12
そういう初歩的な理解があれば
尤度の最大化でモデルのパラメータを決める最尤法はどのような場合に有効なのか？
何か特別な条件が必要なのではないか？
という健全な疑問を持つことができるようになる。
Fisherさんという偉大な統計学者が開発して広めた方法をそのまま信じてしまうのは非科学的。
タグ：
posted at 23:24:30
見込みも可能性も皆無のモデルで尤度が巨大なものを作れるので、それもアウト。
尤度は単に「モデルのデータへのフィッティングの良さ」の指標の1つ。
単にデータにフィットしているだけでは、もっともらしくもなく、見込みも可能性も皆無かもしれないということは統計学の重要な教えの1つだと思う。 twitter.com/itoh_kiminori/...
タグ：
posted at 23:19:17
@komatsu_pretty twitter.com/komatsu_pretty...
このツイートの意味は、「一旦は1つ分×いくつ分じゃないと駄目」という意味でしょうか？
タグ：
posted at 23:09:48
#統計 以下のリンク先では、2つ山の分布から(乱数の分布として)最も出て来やすい1つ山の正規分布を求めています。
2つ山の分布をほんの少し変えただけで、そこから最も出て来やすい正規分布が大幅に変わることがある。その理由はSanovの定理を知っていると直観的に把握できる。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 23:07:16
#統計 続き。この考え方は、事後分布pの変分近似で使われます。事後分布を直接計算するのが大変なときに、何らかの意味でシンプルな分布qで事後分布pから最も出て来やすい分布を見つけるという手法を「変分ベイズ」「変分推論」などと呼びます。
そのときSanovの定理は直観の基本的源泉になる。
↓ twitter.com/genkuroki/stat...
タグ： 統計
posted at 23:03:11
#超算数
2年前にツイートしたこれだけど、著者の意図とは別の意味で楽しめる、という点では、正しく、トンデモ本だね。
落語「一眼国」を思い出す。
著者がアジア諸国の算数教育をどう見ているかを見ることで、日本の算数教育の奇妙さが浮かび上がる。 twitter.com/sekibunnteisuu...
タグ： 超算数
posted at 23:01:06
#数楽 繰り返し書いていることですが，Sanovの定理によれば、大雑把に、分布pのサイズnの標本について
(その標本の経験分布と分布qの違いが小さな確率)
= exp(-n D(q||p) + o(n))
が成立しています。
KL情報量D(q||p)が大きいほど、この確率は速く小さくなります。続く twitter.com/genkuroki/stat...
タグ： 数楽
posted at 22:59:10
@takitoshi BLASにopenblasを使っている以上、これは予測できたことで、CONSISTENT_FPCSR=1でopenblasをmakeし直せば直りそう、という気もしますね。Juliaが使ってるopenblasをどうすればmakeしなおせるのか、はちょっと調べる気がしない笑
タグ：
posted at 22:58:17
知ったかぶって指導要領に書いてあると言い出す人の99%以上は指導要領解説を出してくるの法則（しかも書いてない） twitter.com/rinnetuki/stat...
タグ：
posted at 22:53:51
I started using #JuliaLang in my research. I like the concept of multiple dispatch. An unexpected benefit from using multiple dispatch is that it often removes the need to worry about clashing function names. Functions solve(x) and solve(m::MyModelType) can coexist peacefully. pic.twitter.com/W6xV5kqbTr
タグ： JuliaLang
posted at 22:47:18
こちら、早とちりしました。。all(Cup .> Cdown)がtrueでないと精度保証はされていません。やってみると、n=128からCupとCdownの値が同じになる場合があります。BLASで並列計算するときに、丸めの方向指定に失敗してますね。。 twitter.com/takitoshi/stat... pic.twitter.com/9vrUvELNya
タグ：
posted at 22:45:42
@genkuroki @antimon2 ご指摘通り、強い不等号でやるべきでした。そして精度保証するためには mean(Cup.>Cdown)の値が1でないと精度保証されてはいません…
タグ：
posted at 22:33:32
そりゃ納得いかないですよね。
こういう小学校教育特有の異常な約束事に翻弄される子供たちが、本当に気の毒です。
数学の力ではなく、先生に忖度する能力をひたすら伸ばされている感があります。
#超算数 twitter.com/nao_p_on/statu...
タグ： 超算数
posted at 21:21:33
大阪市の塾代助成のクーポンは、15億円配布するのに5億円かかってますやん。吉村知事は大阪市長もやってましたやん。配るのにお金かけて、学力向上の効果はさっぱりなら、完全な愚策。
維新・吉村氏「完全な愚策」 １０万円給付の巨額経費を批判(産経新聞)
#Yahooニュース
news.yahoo.co.jp/articles/09001...
タグ： Yahooニュース
posted at 21:12:46
完全に超算数の被害者が加害者になるパターンの動画。所々小さな子の声が聞こえてすごく可愛いが、(多分親子?)配信者も子供もとてもかわいそう。 twitter.com/pmanvitaminc1/...
タグ：
posted at 21:12:04
#超算数 わざわざ見に行きましたが【世界で一番成績がよい児童が育つシンガポール算数で、順序が使われていないのは、順序派としては残念】という発言があった。こういうのを贔屓の引き倒し、本末転倒などと称します。 twitter.com/shoyugi/status...
タグ： 超算数
posted at 21:04:40
ごまふあざらし(GomahuAzaras @MathSorcerer
パラメータ設定ファイル群の管理は
TOML.jl (標準ライブラリ)
Preferences.jl
Configurations.jl
ToStruct.jl
Parameters.jl
を組み合わせると良さそうだけれど皆さん何を使ってらっしゃいます？
#Julia言語
タグ： Julia言語
posted at 20:43:08
これを見て「合ってた。よかった」と考える人間が少なからず存在するということが怖ろしい。
#超算数 twitter.com/oomurasakioo/s...
タグ： 超算数
posted at 20:33:05
OokuboTact 大久保中二病中年 @OokuboTact
嘘出鱈目を（論理整合的に）説明できる人はいない。
逆に言えば、論理的整合性を無視すれば、どんな嘘デタラメでも説明可能。 twitter.com/sekibunnteisuu...
タグ：
posted at 20:28:16
@Maiky97588215 #超算数
足し算の順序のような、間違ったことを工夫して教えられても・・・
www.shinko-keirin.co.jp/keirinkan/sans...
タグ： 超算数
posted at 20:23:27
@Maiky97588215 #超算数
かけ算の順序のような、間違ったことを工夫して教えられても・・・
www.asahi.com/edu/student/te...
タグ： 超算数
posted at 20:22:13
超算数を見ていたら超英語に行きついたwww
おれも、小さい頃に習った謎のハワイ人の文字癖がついてグラムだとかメートルとか言われて英語がめんどくさくなって興味を失ったwww
ameblo.jp/hitomi-nakahas...
タグ：
posted at 19:53:00
#超算数
#秋の掛け算順序祭り
今年は大漁！
twitter.com/nao_p_on/statu...
posted at 19:42:42
@antimon2 @takitoshi #Julia言語 Windows 11 環境でも試してみました。
①たし算
②かけ算
③行列の積
強い不等号 < 成立の割合も計算。
それぞれで現れる0.268…や0.5や0.125は統計的にこうなるべき値になっていますか？
Windowsだと色々不安な話題。
github.com/genkuroki/publ... pic.twitter.com/E572d5lvID
タグ： Julia言語
posted at 18:58:54
非公開
タグ：
posted at xx:xx:xx
#julialang
#julia言語
内包表記で書いたこれを、内包表記を使わずfor文中で書くには、どうしたらよいでしょうか？ pic.twitter.com/sfQYWljLbK
posted at 18:12:25
全農広報部【公式】日本の食を味わう @zennoh_food
あの、キャベツって刻むのが大変(高確率でまな板からあふれる)じゃないですか。刻むのがしんど〜って時は、葉を剥がさず丸ごとをくし形に切って煮物にするのも一案。楽だし、たくさん食べられるし、意外と豪華見えします。今回は鶏モモ肉とコンソメ煮に。くたくたのキャベツ、うまいんだよなあ。 pic.twitter.com/9dCb169kSj
タグ：
posted at 18:07:45
#数楽
pronouncekiwi - How To Pronounce
Jensen's Inequality
www.pronouncekiwi.com/Jensen's%20Ine...
「ジェンセン」または「イェンゼン」または「イェンセン」という感じ。
Danish Pronunciationでは「イェンセン」に聞こえる。
タグ： 数楽
posted at 18:02:15
@takitoshi `Cup[:] ≥ Cdown[:]` だと辞書式順序での比較しかしないので、例えば `[1,2,3] ≥ [1,1,999]` とかも true になってしまいますよ。
やりたいことは「全部≥」だと思うので `all(Cup .≥ Cdown)` と書くべきです。
ちなみに手元で↑でやってみたら確かに全部 true でした。
タグ：
posted at 17:38:01
非公開
タグ：
posted at xx:xx:xx
Juliaで丸めモードの向き指定ができない問題。公式のsetrounding_raw関数を使って変えられるのかもしれない。以下のコードの155行目にccallでfesetroundを強引に呼んでいる。
github.com/JuliaLang/juli...
jl_set_fenv_roundingはfesetroundのラッパー。
github.com/JuliaLang/juli... pic.twitter.com/vCCAzehgBS
タグ：
posted at 16:29:05
中村 良幸 (Nakamura Yosh @nakayoshix
1994年から働いていた某建設コンサルで見た津波シミュレーションや港湾の水質解析用のF77のソースコードがまさにそんな感じでしたけど、それから四半世紀経ってもなおまだ大学の研究室内でもそんな感じなのだとしたらそれは悲惨ですね。当時、津波に関してだけはFortran90で全て綺麗に書き直しました。 twitter.com/genkuroki/stat...
タグ：
posted at 15:32:07
中村 良幸 (Nakamura Yosh @nakayoshix
「新しい酒は新しい革袋に盛れ」の現代版ですか…😇 twitter.com/genkuroki/stat...
タグ：
posted at 15:28:33
具体的な実物の例を挙げた方が良いかもしれませんが、ど素人が多方面に喧嘩を売ることになるので、ずるいかもしれないけど避けさせて下さい。
(#Julia言語 のdiscourseでグローバル変数の使用継続にこだわる人のコードをたまに見る。Juliaのような新しい道具を使う所までやって来てもそういう感じ。)
タグ：
posted at 15:23:07
私はプログラミングも数値計算もど素人なのですが、こんなコードを20世紀からの伝統でずっと書き続けているというのはさすがにまずそうだと分かります。
研究室の卒業生の就職にも悪影響を与えそうなレベルでまずいかも。
日本語でツイートしていますが、この話題は日本語圏の話題ではないです。
タグ：
posted at 15:19:24
大学の研究室内での数値計算のコーディングスタイルが
* グローバル変数にモデルのパラメータ達などを全部代入
* グローバル変数を大量に含む函数達
* グローバル変数まみれでしかも長大なmain()函数
になっていることが多いという印象がある。
21世紀のプログラミングでさすがにこれは😭です。
タグ：
posted at 15:19:23
Juliaで作って学ぶベイズ統計学、週末に届いてまだペラペラめくっただけだけど「こういう数式実際どう書くのがいいのかよう知らん」って事柄がたくさん載ってそうなのでちゃんとやろうと思った
www.amazon.co.jp/Julia%E3%81%A7...
タグ：
posted at 14:55:43
これ、大学の研究室での教育が世界的にひどいことになっているのが主な原因じゃないかと思っています。日本語圏だけの問題じゃない。
#Julia言語 のコミュニティでの質問時に示される数値計算のコードが、型を固定したグローバル変数まみれ(constまみれ)になっていることが結構ある。 twitter.com/Hishinuma_t/st...
タグ： Julia言語
posted at 12:59:00
特に自分のダメさを痛感させられるのは、私が積極的に批判している算数教育の問題について、「この問題の解決には権威ある数学者(もしくは権威ある数学の学会)に頼るのが良い」という意見を本気で持っている人を見た場合です。
現実は結構お寒い状態であることを伝えることに失敗している。
タグ：
posted at 12:52:43
こういうのは正直に「自分自身もまた自分の立場でやり難い批判を避ける汚い人間の一人に過ぎないこと」を(ずうずうしく)明言してかつ、やり難い批判も少しはしていることを示すのが良いと思います。
他人のせいにしたりせずに、自分自身がオールクリーンではないことをはっきりさせることが大事。
タグ：
posted at 12:51:44
自分自身の立場に置き換えると、何度か既に述べていることですが、数学の同業者達による算数数学教育分野における悪行の批判は私個人にとってはつらい所で、一応批判はしていますが、結構避けてしまっています。
そこは許して欲しいと思います。 twitter.com/kikumaco/statu...
タグ：
posted at 12:45:07
非公開
タグ：
posted at xx:xx:xx
今井翔太 / Shota Imai@える @ImAI_Eruel
東工大の機械学習の講義で使用された、理論の解説とPythonプログラムの実装・実行が一体になった「機械学習帳」という講義ノートが無料公開されています
chokkan.github.io/mlnote/
解説や図の豊富さやアニメーション、詳細な式展開、実装・実行の簡易さなど、Webの公開資料としては破格の内容です pic.twitter.com/FKIhYXuMz9
タグ：
posted at 12:24:53
#統計 以下のリンク先のヒートマップでプロットされた尤度函数は、理論的にはサンプルサイズ→∞でシンプルな単峰型(2変数Gauss函数で近似される)になるのですが、特異モデルになるパラメータ集合に沿って台が広がってしまっています。
n→∞での結果は安易に有限サイズのデータに適用できない。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 12:22:35
#統計 尤度函数のプロットの例は私のツイログにある。私は自分のツイログ検索のヘビーユーザーで、統計学入門や #Julia言語 関連の情報を得るために非常に便利です。
twilog.org/genkuroki/sear...
私のツイログで「尤度函数 pic」を検索
サンプルツイート
↓
尤度函数の形は単純な単峰型とは限らない。 twitter.com/genkuroki/stat...
posted at 12:17:52
#統計 仮にFisherさんのlikelihoodを日本語に「もっともらしさ」と翻訳していたら、現代においても「もっともらしさ」と言い続けるのは相当に辛かったと思う。
「尤度」と読み方さえすぐに分からない日常で出会わない意味不明の言葉に翻訳してくれたのは不幸中の幸いだったと思います。 twitter.com/stattan/status...
タグ： 統計
posted at 12:02:16
そういう意味では、「一記者の意見」を無理に表明させるよりも、他の第三者の声として、メディア批判の大きさを伝えてもらうためのチャネルを維持することが、建設的じゃないかなと思います。
そう考えると「あなたたちがそれを言うのか」という声をインタビューにそのまま載せたのは英断だったのだな
タグ：
posted at 11:57:30
相互批判不在は非常に大きい問題だとは思います。ただ現実問題、（バズフィードには失礼ながら）「小さなメディア」に所属して、他の全ての大手メディアの向こうを張って批判するのは、とても大変だろうということも想像に難くないわけで。 twitter.com/kikumaco/statu...
タグ：
posted at 11:52:16
岩永さんはHPV問題でメディア批判を避け続けてきたんですよね。これは不自然なほどに避けてきたと思います。
僕は放射能からHPVまで地続きにしか見られないから、「メディアの相互批判不在」は大きな問題と考えています。岩永さんもそこからは踏み出せない twitter.com/y_tambe/status...
タグ：
posted at 11:46:55
OokuboTact 大久保中二病中年 @OokuboTact
落語が教えてくれるダメ教育の見分け方 twitter.com/sekibunnteisuu...
タグ：
posted at 11:43:54
HPVに関しては、メディアの中で唯一頑張ってたのが、バズフィード＆岩永さんであり、逆風の中で頑張り続けてきた医療関係者にとって、心強い存在だったことは（次に似たようなことが起きたときの芽を摘まないためにも）まず認めないと。マスメディア／メディア全体で括ると、その評価も妥当だけど。 twitter.com/kikumaco/statu...
タグ：
posted at 11:43:45
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
小学校の掛け算の教え方のどこかダメかと言えば、
導入で「１つ当たり × 幾つ分という順序で教えること」ではなく、
交換法則を教えながら逆順を間違いにすること。
「１つ当たり × 幾つ分」でも「幾つ分 × １つ当たり」でも同じ意味で、正解なのだから。
タグ： 超算数
posted at 11:37:14
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
＞ 交換法則に入るのは九九表に入ってからだから、先に塾で習って理解しているのか、勘違いしているのか、採点者には分からないよね。
何を言いたいのだろう？
生徒が何を勘違いしているの？
交換法則と何が関係あるの？
塾で予習して何が悪いの？
タグ： 超算数
posted at 11:34:08
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
＞ そもそも、かけ算の順序強制は、文科省が実施する学力調査に関わるものだから、
この人、エビデンスなしで思いつきでデマばっか語っているなあ pic.twitter.com/RMdBsNe3Pv
タグ： 超算数
posted at 11:24:22
亀だってなでられたい。ブラッシングしてくれる女性の前に亀の大行列 karapaia.com/archives/52307...
シンガポールの植物園の池。ブラッシングしてもらうと満足したように帰っていったそう。亀の甲羅表面の下には神経や血管が通っているので、亀はちゃんと甲羅への刺激を感じることができるとのこと。 pic.twitter.com/Tul1Nbh8k3
タグ：
posted at 11:07:20
いやいやいやいや……😌💨(笑)
驚くなかれとか、すくなくないって笑いましたけど、それ以外がいたらむしろびっくりっつー話でしょうが！ twitter.com/sekibunnteisuu... pic.twitter.com/MsFt9OYZt8
タグ：
posted at 10:25:22
@komatsu_pretty 失礼します。
A 行列では一般に交換法則は成り立たない。
B 複素数の掛け算では交換法則が成り立つ。
どちらも正しいことを教えればいいだけでしょう。
Aを教えるためにBを否定して噓出鱈目を敢えて教えるというのは倒錯していると思います。
タグ：
posted at 08:58:20
朝、息子のランドセルからぴらりと出てきたテスト。95点。納得がいかず悔しくて出せなかったらしい。
【問】4人の子供に8枚ずつシールを配るには何枚シールが必要か。
【解】
☓4×8=32枚
○8×4=32枚
大丈夫、母さんも納得いかない。
おはようございます。
最高裁までいく。 pic.twitter.com/NoemLquPtV
タグ：
posted at 08:27:17
須山さんの新しい本、尤度という言葉は他で使ってるからしょうがなく使ってたり、分布の更新は理解のためにはミスリーディング、などなかなか須山節が効いていてよい。
タグ：
posted at 06:58:25
honto電子書籍期間限定セール (12/09)まで 冬☆電書2022テクノロジーとビジネスの最前線 須山 敦志,杉山 将『【期間限定価格】機械学習スタートアップシリーズ ベイズ推論によ...』
30%OFF
amasale.newif.net/sale/hdetail/9...
タグ：
posted at 06:25:14
#数楽 個人的には、リー群論入門書としては
山内恭彦・杉浦光夫、連続群論入門
www.amazon.co.jp/dp/4563003298
が非常に良いと思うのですが、絶賛品切れ中古本の値段がべらぼうに高い！
こういうのは結構困る。
いずれにせよ、SO(3)とSU(2)から入るのが基本定跡だと思います。 twitter.com/toyo9/status/1...
タグ： 数楽
posted at 01:52:02
#統計 KL情報量とは全然違う種類の分布間の違いを測る指標にWasserstein距離がありますが、ℝ上の確率分布の間のp=1のWasserstein距離は特別に簡単で累積確率分布函数の間のL¹距離に一致します。
例の計算については以下のリンク先スレッドを参照。
↓ twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:35:29
#統計 その説明の仕方を発展させれば様々な直観が得られます。
個人的にちょっと面白いと思う例については以下のリンク先を参照。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 00:29:48
#統計 リンク先の説明はいまいち。これだとKL情報量の非対称性を説明できないし、KL情報量の基本声質への直観も得られない。
個人的にはSanovの定理の線に沿った説明がよいと思う。続く twitter.com/of_complex/sta...
タグ： 統計
posted at 00:24:55
#数楽 何度も繰り返し述べていることですが、行列を基底を取ると得られる線形写像の表現行列としてのみ理解しようとするのは誤りだと思います。
双線形形式も基底を取れば行列で表現できるし、Hesse行列として自然によく現れます。
いずれにせよ、行列は何かを表現するために使われる。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 00:09:15
