黒木玄 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年11月05日(木)
#Julia言語 を #R言語 側から使ってみることを試したい人はリンク先動画が参考になりそう。 twitter.com/juliaconorg/st...
posted at 00:32:58
Hideki Kawahara: WAS @hidekikawahara
55年以上前のビデオを見ると,今と(ほとんど)同じマイクを使っている.音の(アナログ部分の)ハードウェアは恐ろしいほど進歩していないように見えるが本当か?これでいいのか?
Dave Brubeck - Take Five youtu.be/tT9Eh8wNMkw @YouTubeより
タグ:
posted at 01:14:25
Ninja DAO | CryptoNi @CryptoNlnjaNFT
数年前にアメリカの大学で講演をした時に、同時に行われたイベントで「模擬国連のジャッジ」をさせてもらったのですが、そのシナリオが「国連で主要国が朝鮮半島問題について話している間に、次から次へとトランプの不規則ツイート発言が飛んできて各国が振り回される」というシュールなものでした。 twitter.com/order1914/stat...
タグ:
posted at 01:15:48
Ninja DAO | CryptoNi @CryptoNlnjaNFT
確か基本ストーリーは韓国で原発事故が起こって、その混乱に乗じて北朝鮮が攪乱戦略を取ろうとしている状況の中、日米中ロの4カ国が対策を協議するというもので代表団が北朝鮮封じ込めで合意寸前に至った時に、”Kim Jong Un is a smart guy, I trust him!"みたいなツイートが飛んできて混乱していた。
タグ:
posted at 01:24:10
統計学にょ哲学するの人、ポストトゥルースだからとか言っておちょくるあたりめちゃくちゃ気に入らない
アカデミアという高級な論理に正しさを(ある程度)保障されたコミュニティーでポストトゥルースとか言っちゃうの、自分が正しいって取り巻きから安心感を得たいいつもの哲学者の構図と同じやん
タグ:
posted at 01:26:58
テレビをつけパなしにしていたら、ケーブルTVのCMで、スマートスピーカー使い方的なものがあり、我が家の電気が消えたり、その他もろもろ部屋のスマート家電が反応し、深夜作業中に焦る。
タグ:
posted at 01:32:46
#統計 こんな感じにグラフを描けば
* AICによるモデル選択は、データが運悪く偏っている場合に脆弱で、そういう場合には大外ししている予測分布を与えるモデルの側を選択してしまう
* AICの差と(実践的場合に知ることができない)真の予測誤差は逆相関している
ということがすぐにわかる。 pic.twitter.com/6hzlPaLiwq
タグ: 統計
posted at 01:35:16
The Julia Language @JuliaLanguage
Check this out: A MIT-led team will be developing software to help forecast space storms... The coolest part, they will be using #JuliaLang to do so! Congrats to @AlanEdelmanMIT and the other @MIT_CSAIL + @MIT collaborators.
See the full article here: news.mit.edu/2020/mit-led-t... pic.twitter.com/56DmcgRf4h
タグ: JuliaLang
posted at 02:00:49
非公開
タグ:
posted at xx:xx:xx
財務省だけの責任にするのは良くないです。これは財政制度審議会の委員の責任です
お名前を貼っておきます twitter.com/y_kaneko/statu... pic.twitter.com/a8HCkvd0Mo
タグ:
posted at 02:33:55
私は,実在あるいはその反映としての認識の一番深いところに,普遍的に存在する記号的な何かがあり,これが数学的な本地だと感じています.この記号的な何かに数学的形式を与えたものが「代数的構造」といえます.(佐藤幹夫) pic.twitter.com/jqlBPaDVnV
タグ:
posted at 02:36:45
この記事に対して「どの面下げて…」「お前が言うな」というような反響が凄くて、それは全くその通り、特に犠牲者を間近に見ておられるお医者さんたちの怒りの程は察するに余りあるのですが、朝日新聞の責任を追求するのとは別に、この記事を書いた記者さん(編集委員)は応援してあげたいと思います。 twitter.com/asahi/status/1...
タグ:
posted at 09:52:15
特にこの方、去年もこういう記事を書いています。
この件に関して朝日新聞の責任はもちろん追求されるべきだけど、「記者としてこの件には触れずにいた方いい」なんてことになって誰もが逃げるようになったら、結果として更に犠牲者を増やす結果にもなりかねませんし。
www.asahi.com/sp/articles/AS...
タグ:
posted at 09:52:15
OKUMURA, Akira(奥村 曉) @AkiraOkumura
選挙のたびに統計誤差の概念は少しずつ広まってきている印象がある。次は系統誤差だ、こっちは難しいぞ。実験であれば装置の特性を熟知している必要があるし、選挙であれば選挙区ごとの特性や、事前投票の層がどうなっているかを理解していないといけない。
タグ:
posted at 09:58:02
OKUMURA, Akira(奥村 曉) @AkiraOkumura
難しすぎて、物の本には長さの狂った二つの物差しの例え話くらいしか出てこない。あれはどうにかしないといけない。
タグ:
posted at 09:59:37
#統計 KL情報量の
D(q||p) = ∫q(y) log(q(y)/p(y)) dy
とD(p||q)のp=qの近くでの比較は
p=q+h, ∫h(y)dy=0
とおいてhについて展開すれば分かります。
q log(q/(q+h)) = -h + h²/q - h³/(3q²) + …
(q+h)log((q+h)/q) = h + h²/q - h³/(6q²) + …
hの1次の項は積分すると消え、2次の項は一致。 twitter.com/mshero_y/statu...
タグ: 統計
posted at 10:03:46
これも去年の記事。
後半(有料部分)、「副反応」に対する評価がちょっと腰が引けているけど、
「効果示す研究発表 相次ぐ」とか良く書いてくれた。
www.asahi.com/sp/articles/AS...
タグ:
posted at 10:07:10
大隅典子@東北大学(『小説みたいに楽しく @sendaitribune
以前にこちらの『日本史サイエンス』の紹介記事をnoteに書いたら、担当編集者の方が新聞広告用に使ってくださって、今日の日経朝刊にも掲載♫ note.com/sendaitribune/... pic.twitter.com/G4a2T3wy7M
タグ:
posted at 10:17:41
Masa Yamamoto予測誤差が大き @mshero_y
先生、ありがとうございます。
等号が成り立つ時は?
まずは確かめるべき条件ではあります!
pとqが近い時の入れ替えでの挙動、自分の分野だとそれだけでも使い道が見えてきました。
ありがとうございました。 twitter.com/genkuroki/stat...
タグ:
posted at 10:17:43
大隅典子@東北大学(『小説みたいに楽しく @sendaitribune
実は、拙著『脳からみた自閉症 「障害」と「個性」のあいだ』のご担当だったという偶然♫ www.amazon.co.jp/dp/B01ESWLF1I/...
タグ:
posted at 10:18:57
ごまふあざらし(GomahuAzaras @MathSorcerer
Juliaのノートブックを作る Pluto.jl を Binder で動かそうというプロジェクト.リポジトリのオーナーは Pluto.jl のメンテナーなので期待が膨らみます(手元の💻で動作も確認ずみ)
github.com/fonsp/pluto-on...
#Julia言語 pic.twitter.com/Xh9QoMMpF1
タグ: Julia言語
posted at 10:19:58
ごまふあざらし(GomahuAzaras @MathSorcerer
Pluto.jl は教育用途でも使われています.
www.youtube.com/watch?v=rRCGNv...
タグ:
posted at 10:19:58
ごまふあざらし(GomahuAzaras @MathSorcerer
ENV["PYTHON"] = Sys.which("python3")
using Pkg; Pkg.build("PyCall")
これで行けますかね? twitter.com/zillione/statu...
タグ:
posted at 10:50:14
#Julia言語
ENV["PYTHON"] を設定するだけでは意味がなくて、build PyCall することが大事。
github.com/JuliaPy/PyCall... にアクセスして Pkg. build に「下線を引いて」読み直すと理解が進む可能性あり。 twitter.com/mathsorcerer/s...
タグ: Julia言語
posted at 10:55:39
ごまふあざらし(GomahuAzaras @MathSorcerer
おっしゃる通りで PyCall に限らず おかしいなぁ〜設定したんだけどなぁ〜と思ったときは Pkg.build すれば解決しますね(笑)
時間がある人は deps/build.jl を見ると理解が深まります.
github.com/JuliaPy/PyCall... twitter.com/genkuroki/stat...
タグ:
posted at 10:59:03
ごまふあざらし(GomahuAzaras @MathSorcerer
計算機 A でうまくいくのに同様のOSの計算機BでどうしてもインストールやDLLの読み込みがうまくいかず行き詰まったときは
~/.julia
を消して計算機Aで行なったインストール作業するとBでもよくなったりする(白目)
#Julia言語
タグ: Julia言語
posted at 11:19:48
#統計 赤池弘次さんの1980年の2つの論説はこれ↓
ismrepo.ism.ac.jp/index.php?acti...
統計的推論のパラダイムの変遷について
www.jstage.jst.go.jp/article/butsur...
エントロピーとモデルの尤度
数え切れないくらい繰り返し、みんなに勧めて来た。
以下ではこれらを読んでいることを前提とします。続く
タグ: 統計
posted at 13:02:48
#統計 私による「通常のP値」と「ベイズ統計での事後分布を用いたその類似物」の同時プロットの例は以下のリンク先などで紹介しました。
以下のリンク先の例では片側検定のP値とベイズ統計での類似物が数値的によく一致していることを示しています。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 13:17:42
#統計 1980年の赤池弘次さんによる尤度概念の解説を読めば、実際に尤度の概念について一生のあいだ誤解し続けた人達が20世紀に沢山いたであろうことを想像できると思います。
ismrepo.ism.ac.jp/index.php?acti...
統計的推論のパラダイムの変遷について
www.jstage.jst.go.jp/article/butsur...
エントロピーとモデルの尤度
タグ: 統計
posted at 13:38:35
ごまふあざらし(GomahuAzaras @MathSorcerer
素晴らしすぎてワロタ
ronisbr.github.io/PrettyTables.j... pic.twitter.com/Nt040dVE8V
タグ:
posted at 13:44:27
ごまふあざらし(GomahuAzaras @MathSorcerer
LaTeXバックエンドもサポートされているっぽい.
ronisbr.github.io/PrettyTables.j... pic.twitter.com/Jz9KXtfd9u
タグ:
posted at 13:47:01
ごまふあざらし(GomahuAzaras @MathSorcerer
もうそろそろでメジャーリリースされるようなのでフィードバック待ってるよ〜とのことです.
タグ:
posted at 13:47:02
#統計
#Julia言語 で計算して来ました。
nbviewer.jupyter.org/gist/genkuroki...
添付画像
1. n=200での二項分布の正規分布での近似。ほぼぴったり一致。
2. n=20での二項分布の正規分布近似も悪くない。
3. nが200周辺での負の二項分布の正規分布近似。
4. nが20付近だと負の二項分布の正規分布近似は悪い。 pic.twitter.com/9K1IiM3NcY
posted at 14:34:27
#統計 以上によってやっと
oku.edu.mie-u.ac.jp/~okumura/stat/...
に書いてある「止め方で仮説検定の結果が変わる」という話を扱う準備ができました。そこでは、
* コインを投げる回数を12に固定した場合(表は3回出た)
と
* 表がちょうど3回出るまでコインを投げ続けた場合(その回数は12になった)
~続く
タグ: 統計
posted at 14:44:12
#統計 ~の場合の片側検定のP値がそれぞれ 7.3% と 3.3% になり、大きく違うことが紹介されています。
7.3%は二項分布モデルの、3.3%は負の二項分布モデルのP値です。
添付画像はその計算を
nbviewer.jupyter.org/gist/genkuroki...
で確認した場面です。念のために2通りの方法で計算して一致を確認しています。 pic.twitter.com/S9cbi0px4k
タグ: 統計
posted at 14:44:16
ごまふあざらし(GomahuAzaras @MathSorcerer
github.com/terasakisatosh...
ので参考にして自分でも弄ったプライベートにしたのをパブリックにします。
(Juliaのsys imageに自前の関数を焼いて追加し、それをCの関数から呼ぶというもの)
タグ:
posted at 14:45:53
#統計 詳しくは
oku.edu.mie-u.ac.jp/~okumura/stat/...
の解説を見て欲しいのですが、
* コインを投げる回数nを固定して表の出た回数kを観測する場合
と
* 表がちょうどk回出るまでにコインを投げた回数nを観測する場合
では、使うべき適切なモデルの確率分布族は違います。このことは高校数学レベル。
タグ: 統計
posted at 14:49:45
午後眠いけど寝ずにNewton法の続き。normがなかなか小さくならないのは何故?ということで解から定数項を逆算することに。なんか計算式が間違っているような気がしてきた。あとループ回数はかなり大きくしないとダメだね。 #julialang
タグ: julialang
posted at 14:52:09
#統計 添付画像は上のn=12の場合に対応する計算をn=120,1200,12000の場合にもやってみた結果です。
2つのモデルの片側検定のP値の違いが、nを大きくするにつれて小さくなって行く様子が見えています。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/zsQw67T5pI
タグ: 統計
posted at 15:04:39
#統計 せっかくなので、各nごとにkを動かして2つのモデルの片側検定のP値を重ねてプロットしてみましょう。
nbviewer.jupyter.org/gist/genkuroki...
nが小さなときの2つのP値の食い違いは大きいが、nを大きくするとぴったり一致するようになることが分かります。 pic.twitter.com/uakmmShTPV
タグ: 統計
posted at 15:12:12
#統計 さらに、以下のリンク先を見れば分かる様に、二項分布モデルの片側検定のP値とそのベイズ統計版での類似物もnが十分大きければ、近似的に等しくなります。
以上で示したnを大きくしたときの近似的は正規分布近似(中心極限定理)を経由すれば一般的に示せます。続く twitter.com/genkuroki/stat...
タグ: 統計
posted at 15:22:28
#統計 訂正
❌nを大きくしたときの近似的は
⭕️nを大きくしたときの近似的一致は
正規分布近似はほとんど空気を吸うごとく使えるようになっておくと便利です。コンピュータでグラフをプロットして確認した方が経験値が楽に増えます。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 15:44:04
ごまふあざらし(GomahuAzaras @MathSorcerer
sysimage を動的にロードしてるってことでいいのかな?
github.com/GunnarFarnebac...
手元ではテストコードは動いています.
タグ:
posted at 15:46:32
#統計 重要な補足
モデルを固定すればn→∞で尤度は「もっともらしさ」の正しい指標になります。
しかし、サンプルサイズnを大きくしたお陰で詳細な構造が見えて来そうだと考えて、モデルの側をパラメータの多いより複雑なものに変える、というようなことをすると、そうはなりません。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 15:50:25
今個人的に関心を持ってる言語、、、NimとJuliaかなあ、、ちゃんとやったことないので、、あとSwiftは色々入ったタイミングで触ってみたい。
goとかRustはここで使うとバッチリ行けるなみたいなのが見れたので、あの利点をいろんなところにそれぞれのドメインに対して特化させた形で移植してみたい。
タグ:
posted at 16:06:20
Nice complement to #Jupyter - #Pluto. Written in #JuliaLang. Main feature, I think, is tracking of dependencies and updating proper cells on change.
github.com/fonsp/Pluto.jl
posted at 16:44:06
ごまふあざらし(GomahuAzaras @MathSorcerer
DynamicallyLoadedEmbedding.jl works on my RPIZero
discourse.julialang.org/t/ann-dynamica... pic.twitter.com/afr7m20nIG
タグ:
posted at 17:09:17
非公開
タグ:
posted at xx:xx:xx
中村 良幸 (Nakamura Yosh @nakayoshix
今日はハロワに寄った後で「1から始めるJuliaプログラミング」を買いに紀伊国屋書店に来たんですが、いくら探してもみつからないので検索したところ、Pythonコーナーに置いてありました。
周囲には今話題の「Pythonによる」で始まる本が16種類も置いてあり、噂は本当だったんだな…と実感しました。😇 pic.twitter.com/kUiT5Ybfl8
タグ:
posted at 18:05:36
#統計 さすがに、「確率分布を特徴づけるパラメータ」という言い方が頻出する分野について語ろうとしているのに、【確率変数が持つ分布を特徴付ける値を、その期待値~という】(p.31)と説明しちゃうのはまずい。
これ一ヶ所だけが雑なだけなら見逃しても良いが、ほとんど稠密に説明が雑でひどい。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 18:32:44
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
#統計 p.83には【分布族が対象を十全にモデル化】していなくても【ある弱い前提さえおけば、ベイズ流の更新プロセスは最終的に真理へと到達しうる】と書いてある!(おそらくEarman, 1992を不適切に引用している)
もう「ベイズ主義」とか言ってもいいから、こういうのは勘弁して欲しいです。 twitter.com/genkuroki/stat... pic.twitter.com/OSv3YS2TCO
タグ: 統計
posted at 18:38:26
#統計 ベイズ統計は数学的に難しいので色々雑になるのは仕方がないかもしれないが、
p.17の図1.1【「Major axis」と表示されているのが回帰直線】‼️
p.139に【最尤法の他】に【最小二乗法】があるかのように書いてある‼️
pp.142-143では、回帰の誤差項εの平均がμ‼️ twitter.com/genkuroki/stat... pic.twitter.com/GTfUUYwDhg
タグ: 統計
posted at 19:00:19
クルーグマンはQuoraではテクニカルな話はほとんどしないんだけど、これに読者個々人がどう反応するかがリトマス紙みたいに機能しておもしろいかなって思ってやったんだよね
タグ:
posted at 19:15:22
非公開
タグ:
posted at xx:xx:xx
一般に本を読んだ結果を自分が理解しているかを確認するための方法として、「その本に書かれていない例を自分で作って説明できるか?」と自分に問うことは鉄板の定跡だと思います。
単なる印象を受け入れただけでは理解したことにならない。
タグ:
posted at 19:28:06
Noda takayukiさんによる「鬼滅の刃の作者は、どんな人ですか?」への回答 jp.quora.com/%E9%AC%BC%E6%B... すごい!他人とは思えない
タグ:
posted at 19:32:04
willseeds.txt-nifty.com/omoi/2020/11/p...
>「その550gって問題に出てきてないよね?どうやって考えて110gってわかったの?」と尋ねたところ、なぜか完全にフリーズ。何度声をかけても答えてくれることはなく、
タグ:
posted at 20:05:28
>どう考えたのかは分からぬまま。それでも、ノーヒントで解けたのだから、恐らくきちんと考えたのは確かだろうと、ひとまず次の問題に進んでもらったものの、そこでおかしなスイッチが入ってしまったのか、その後引き続き不調のまま浮上せず。
タグ:
posted at 20:06:14
「その550gって問題に出てきてないよね?どうやって考えて110gってわかったの?」を「そのやり方は駄目」と受け取ったんじゃないかな?
こんなこと言っている人だからね。
willseeds.txt-nifty.com/omoi/2009/06/p...
タグ:
posted at 20:11:51
私なら、正しい答えが出ているなら、どうやって求めたのかは問わないで、類題を出す。
それで悉く正解を出したなら、「確実に理解している」と判断する。
それから、どうやって解いたかを聞く。
しかし明確に答えられるとは限らない。その場合は、深く詮索しないで次に行く。
タグ:
posted at 20:15:39
twitter.com/sekibunnteisuu...
>類題を出す。
それで悉く正解を出したなら、「確実に理解している」と判断する。
一部訂正。
類題+αを出さないと「確実に理解している」までは言えないな。
タグ:
posted at 20:34:44
( a )mで( b )kgのロープ( c )mの重さは?
が解ける子には、
一辺の長さが( a )mの正方形の板の重さが( b )kgだとすると、一片の長さが( c )m正方形の板の重さは?
も出すようにしている。
タグ:
posted at 20:35:39
非公開
タグ:
posted at xx:xx:xx
willseeds.txt-nifty.com/omoi/2006/05/p...
>私自身、つい数年前まではそうでした。
子どもに関わる仕事をして、実際、その単元についても何度も何度も指導してきたにも関わらず、意味をあまり考えることがなかったのです。
タグ:
posted at 21:05:36
深く考えることなく「そうなるもの」と覚えていた人が、教える立場になって、おかしな考えに嵌って「ちゃんと理解するようになった」と誤認するというパターンがあるようだ。
タグ:
posted at 21:07:28
以前やり取りした数教協の人も、
「足し算の意味とか考えなかったので算数・数学が苦手だった。今は、増加と合併をしっかり理解している」
と言っていた。
頭が痛くなる。
タグ:
posted at 21:08:41
#DancingLinks な #AlgorithmX の #ruby での実装は #LuaLaTeX と #MetaPost での経過の可視化で遊び続けていてなかなか進まないけど漸く少し目処が付いてきたので前に進まないとな。 www.instagram.com/p/CHNcF1Bg9MV/...
タグ: AlgorithmX DancingLinks LuaLaTeX MetaPost ruby
posted at 22:00:53
須山敦志 Suyama Atsushi @sammy_suyama
Turing.jlで欠損値付きデータに対して線形回帰を行う - my blog tetuppei.hatenablog.com/entry/2020/10/...
タグ:
posted at 22:11:48
須山敦志 Suyama Atsushi @sammy_suyama
lightGBなんちゃらが欠損値を扱えるとかわりとどうでもいい話をよく聞きますが、基本的にあらゆる回帰モデルは容易に欠損値補間ができるように拡張できます。もちろん、欠損の捉え方により無数にバリエーションはありますが。
タグ:
posted at 22:31:16
この大統領選挙(およびトランプ大統領の訴訟戦略)に関連して、「ゲーデルが、民主主義国アメリカが独裁化することが可能だと言った」と書かれる人がtwitter上に徐々に増えてきているようですが、コレは全くそう言う話ではありませんからね。
ytb-logic.blogspot.com/2012/08/blog-p...
タグ:
posted at 22:31:19
#Julia言語 #Turing
Turing.jlで欠損値付きデータに対して線形回帰を行う - my blog tetuppei.hatenablog.com/entry/2020/10/...
posted at 22:36:27
ネバダ州の選挙管理委員会のスタッフが記者会見する中、突然乱入し「バイデンは選挙結果を奪おうとしている。バイデン家は不正だらけだ」と絶叫する男性。このシャツのデザインのメッセージ性にやられた。 twitter.com/evanmcmurry/st...
タグ:
posted at 23:25:42