黒木玄 Gen Kuroki(@genkuroki)/2022年03月23日

#統計正規分布モデルに代表される統計モデルによって、現実がぴったり記述できると考えるのは馬鹿げています。

だから、統計モデルが的を外していたときに何が起こるかについて(教科書の内容とは別に)気にする癖をつけておく必要がある。

そのための題材としてこのスレッドの内容が使えます。

タグ：統計

posted at 23:58:02

#統計無作為抽出で得られた標本の平均と不偏分散がどのように確率的に揺らぐか、という基礎的な話題の理解への近道は #Julia言語などを使って様子をグラフに描いてみること。

それを実行した人は知らないうちに、真の分布が統計モデルで実現不可能な場合の最尤法の様子を見てしまっているわけです！

タグ： Julia言語統計

posted at 23:53:48

#統計以上では、大量の標本の平均と不偏分散の組の散布図を描いているだけなので、平均と分散の記述統計という最も基礎的な話題に過ぎないのですが、実質的に

　正規分布モデルの最尤法の漸近挙動

を見ていることになります！

しかも、真の分布が統計モデルで実現不可能な場合になっている！

タグ：統計

posted at 23:48:27

せめてゼレンスキー大統領の原語の音声（同時通訳を被せていないやつ）が聞ければ、どこが和訳から抜けているかぐらいは僕にも判断できるんだけど…
（僕は琉球語ぐらいにはウクライナ語が理解できる、はず。）

タグ：

posted at 23:30:22

あれ、アエラの全文翻訳、文面が変わってない？日本語は自然になったけど、内容が薄くなったような…

タグ：

posted at 23:25:35

ToruStark トールスターク @ToruStark

どうなんだろう。明日辺りには大統領の演説全文の公式の翻訳が出るんだろうか。出せるんだろうか。

タグ：

posted at 23:19:12

マイクラのステーキが美味しそうに見えるｗ pic.twitter.com/GCg1IG9J9l

タグ：

posted at 23:11:49

ナザレンコ氏の翻訳がなければわからないことばかりだった。今日の朝に、日本におけるウクライナ語の専門家の不足を心配してみたけど、実際に日本は語学の専門家の深刻な人材不足を抱えていると思う。ロシア語の専門家だって多くはないんだぜ。

タグ：

posted at 23:04:14

すごいレトリックだ…

タグ：

posted at 22:46:37

要するに、日本もウクライナと同じようにロシアに領土を奪われている。その領土（北方領土）は返還されるべきだ。ウクライナが海上封鎖されつつあるように、日本も他の侵略国家に海上封鎖されるかもしれない。それを起こさせないように、侵略の企みは失敗し高い代償を払わされることを共に示せと。

タグ：

posted at 22:45:58

この戦争は資源を持たず貿易に頼る国々（当然日本も入る）への被害は甚大だと言ったり、ロシアは、海上封鎖によって圧力かけるという方法を、侵略を企んでいるかもしれない他の国に示したと言ったりして、かなり過激に踏み込んだ発言をしているじゃないか。

タグ：

posted at 22:38:51

F. Takagi @cactus_f

とある旅人（とあるヒモ） @toaruhimo

ゼレンスキー大統領のスピーチライター超優秀よなぁ

タグ：

posted at 22:37:17

温かい皆様のおかげでただの無職一般人の僕も8,000フォロワーを達成することができました
ありがとうございます

何もしないと決めてた1年のギャップイヤーも今月で終了します
4月からは日本、世界問わず暴れたいと思いますのでどうぞよろしくお願いいたします pic.twitter.com/wzQ4WtnU0l

タグ：

posted at 22:36:13

adhara_mathphys @adhara_mathphys

天むす名古屋 Temmus @temmusu_n

Heine-Stieltjes多項式の根とvan Vleck多項式の係数の間には次のような関係式が成立しています． pic.twitter.com/AYwkIhTQif

タグ：

posted at 22:12:30

改正教育基本法の毒が回りまくっている。47年の原始条文に戻すべき。教育勅語を現職の国会議員が称揚するなんて、あるべきではない事態ですよ。

タグ：

posted at 22:08:14

#統計添付画像は、左半分がサイズnの標本の平均と不偏分散の分布で、右半分は中心極限定理による2変量正規分布による近似です(少し上の方の「答え」の結果を使っている)。

n=320くらいになるとその近似がかなりうまく行っていそうに見えます。

github.com/genkuroki/publ... pic.twitter.com/3UeGpwiYo7

タグ：統計

posted at 21:30:56

Julia日本語記事まとめ @julia_kizi

#統計そういう地道な計算を実際にやり始めてみたり、やってみた人であれば、結構面倒だったと感じると思います。

しかし、コンピュータで数値的な確認を気軽にできれば、面倒な計算を安全にとばせる可能性も出て来る。

#Julia言語のコード↓
github.com/genkuroki/publ... pic.twitter.com/8lCjZFj3Ss

タグ： Julia言語統計

posted at 21:27:09

新しい記事がQiitaにアップされました！#Julia言語 #Julia日本語記事
 qiita.com/1s22s1/items/6...👈

タグ： Julia日本語記事 Julia言語

posted at 21:05:20

れたす⚪️ @circle_retas

相手の動物がレベルで表されてて可愛らしいのもいいね
ぷよ碁の次はこれが最適解っぽい

ぷよ碁→これ→COSUMI のステップアップカリキュラムが良さそう
このサイトキャッチーな名前がないので呼び名に困る
#囲碁 pic.twitter.com/cnjPRX3z8n

タグ：囲碁

posted at 19:59:08

れたす⚪️ @circle_retas

ここやってみたけど最後の二眼だけの状態みたいになっても一眼にしたりしないし
最弱でも比較的まともなAIでよい twitter.com/mhlpanda_limin...

タグ：

posted at 19:56:17

#統計証明には特別な工夫は必要なくて、定義に戻って地道に計算すればよい。μ=0、σ=1の場合に帰着することを使えば書く分量は少し減らせる。

タグ：統計

posted at 18:31:14

もるてん @iitenki_moruten

#julialang

function makerand(n)
rng = MersenneTwister(200)
for i=1:n
println(rand(rng, Uniform(0, 2.0)))
end
end

こういう関数を作って実行すると，Uniformの範囲が[0,1.0]のときのちょうど２倍の値の乱数が出力されるのは仕様ですか？ pic.twitter.com/2fN5hcLF0a

タグ： julialang

posted at 17:42:33

#統計答え続き

var(X̅) = σ²/n

cov(X̅, S²) = (κ₃/n)σ³

var(S²) = (κ₄/n + 2/(n-1))σ⁴

例：Dが正規分布ならばκ₃=κ₄=0になる。そのとき、標本平均X̅と不偏分散S²の共分散が0になるだけではなく、X̅とS²は独立な確率変数になる。

タグ：統計

posted at 17:34:32

#統計答え続き

分布Dのiid X_1,…,X_nの標本平均X̅と不偏分散S²の定義

X̅=(X_1+…+X_n)/n
S²=((X_1-X̅)²+…+(X_n-X̅)²)/(n-1)

これらはそれぞれ母平均μと母分散σ²の不偏推定量である。すなわち、

E[X̅]=μ

E[S²]=σ²

続く

タグ：統計

posted at 17:34:31

#統計答え

Xが分布Dに従う確率変数のとき、分布Dの平均μ、分散σ²、歪度κ₃、尖度κ₄の定義は

μ=E[X]
σ²=var(X)=E[(X-μ)²]
κ₃=E[((X-μ)/σ)³]
κ₄=E[((X-μ)/σ)⁴]-3

期待値0の確率変数A,Bの共分散の定義は

cov(A,B)=E[AB]

続く

タグ：統計

posted at 17:34:30

#統計 ℝ上の分布よりも、ℝ²上の分布の方が複雑さが増して面白くなります。

基本的な演習問題：母集団分布の平均μ、分散σ²、歪度κ₃、尖度κ₄を使って、サイズnの標本の標本平均と不偏分散の分散と共分散を表せ。

演習問題の続き：上のプロットとこの演習問題の結果を比較せよ。

タグ：統計

posted at 15:37:13

#統計続き

観察

1. 大数の法則によってnを大きくすると、(標本平均, 不偏分散)のばらつきは小さくなり、真の値に集中して行く。

2. 中心極限定理によってnを大きくすると、(標本平均, 不偏分散)の分布は2変量正規分布で近似されるようになる。この場合には標本平均と不偏分散は相関している。 pic.twitter.com/raEZj06fVq

タグ：統計

posted at 15:19:37

#統計 #Julia言語

分布 MixtureModel([Normal(), Normal(20,1)], [0.95, 0.05]) のサイズn=40,80,160,320の標本を1万個生成して、標本平均と不偏分散を計算してドットでプロット。橙のドットは真の値。 twitter.com/genkuroki/stat... pic.twitter.com/ew3xAsl2U0

タグ： Julia言語統計

posted at 15:19:34

Yuichi Tanaka @ychtanaka

一番問題なのは旧スーパードライとの比較がないことなので、この画像みたいなグラフだったらまだわかる。なぜパッケージにするときに旧版との比較を消した…。

special.nikkeibp.co.jp/atclh/NXR/22/a... pic.twitter.com/ACw8k1FpFI

タグ：

posted at 12:45:39

朝日新聞デジタル @asahicom

やる気そぐ｢測りすぎ｣　成果主義の落とし穴、研究の世界にも
www.asahi.com/articles/ASQ3Q...

春は人事異動の季節。米国流の成果主義ですが、その米国では成果主義の評価制度を廃止するＩＴ企業が増えているそうです。記者コラム「多事争論」です。

タグ：

posted at 10:30:04

Daisuke Iizuka @diizuka

@kikumaco 車椅子の東大の先生が、311の時にEVが止まり建物外に避難できなかったので、ああ、障害者支援は健常者の余力がある時にだけなんだと感じだ、みたいな事を話していたのを思い出しました。私はスマホでオンライン音声認識しないと聞き取りできないので、災害時でもネットと電気が必須です

タグ：

posted at 09:52:43

kayoちゅう。 @KAYO_chu_des

１日１ぷよ碁

16対19で敗北しました。 puyogo.app/rp?kf=M1RTJTVE... #ぷよ碁

タグ：ぷよ碁

posted at 09:48:52

ふにふにかめゆき® @funikame

うさぎ林檎＠ししょーPPMP @usg_ringo

@J_kaliy ちなみに前ツイートのグラフは、news.yahoo.co.jp/articles/e5f28... より。

タグ：

posted at 09:35:17

うさぎ林檎＠ししょーPPMP @usg_ringo

@sayakatake 日本助産評価機構が適格認定している助産院をのぞいてみたけど、まぁ自然がどうのとかがウリなんだよな。でもそれって本来の安全なお産とは直結しないんだよね。
そんでカイロプラクティックとか代替医療の名前もちらついてた。
www.josan-hyoka.org/midwiferypract...

タグ：

posted at 09:31:31

@sayakatake 助産師の開業権が歪みを生むんだよねぇ。分娩では侵襲性の高い医行為が必要な場面があるのに医行為が出来ない助産師に独立した場所で行わせればどうしたって無理が出るよ。自律のために医行為紛いの施術を取り込みやすいのも不味い。

タグ：

posted at 09:27:38

ふにふにかめゆき® @funikame

@J_kaliy 去年の第5波で感染して(8/17陽性判明)重症化ギリギリ手前まで行って帰還した私が通ります。政府は昨年7/12から東京都に4回目となる緊急事態宣言を出し、東京五輪開幕まで2週間を切っていました。それでも五輪は開催。その後の感染者数はこの通りですね。 pic.twitter.com/8VXakCRNwt

タグ：

posted at 09:11:01

#ぷよ碁実戦詰碁(正解は1通りではない)

白先如何？ (正解は上の方に書いた)

すべての変化を読み切るのは私の棋力ではかなり大変。 pic.twitter.com/GtVL5DyAaj

タグ：ぷよ碁

posted at 08:52:36

5路盤だと最後の勝負所の局面がガチの詰碁になってしまい、ルールを覚えたての人にとっては難解過ぎる局面になるので、7路盤の方が「相手よりも支配領域の面積を広くする」という戦略で勝ち易いということが十分にあり得ると思います。

タグ：

posted at 08:41:58

人間だと悪手を打たずに済むことはなくて、最後の勝負所で正しく打てた方が勝つ。

勝負強い人はそういう場面で間違えない。

私は性格的に無理(笑)

タグ：

posted at 08:30:46

まるで作られたドラマのように逆転が繰り返される大激戦。

この棋譜、囲碁漫画とかでの囲碁を覚えるシーンで使えるんじゃね？

自分も激戦してみたい人は
↓
ぷよ碁
puyogo.app

タグ：

posted at 08:27:38

daikatsu yasuyuki @yasudaidai

#プログラミング勉強中
 #julialang #julia #julia言語
 #random
#vscode #jupyternotebook pic.twitter.com/apH2ujkygn

タグ： julia julialang julia言語 jupyternotebook random vscode プログラミング勉強中

posted at 08:25:06

#ぷよ碁続き

添付画像

①黒の再再悪手で再再再再逆転！！！惜しい！

②白がついに正解手を打ってしまった！

③終わりの局面。白が劫を取った直後。同形反復禁止のルールによって黒はすぐに取り返せない。黒の負けになっている。

最後は劫で決着する大激戦！

puyogo.app/rp?kf=MzJCIyRD... pic.twitter.com/srx8x1WpH4

タグ：ぷよ碁

posted at 08:23:24

#ぷよ碁続き

添付画像

①黒の悪手。これで逆転している。△に打てば白が全滅していた。

②白の悪手で再逆転！

③黒の再悪手で再再逆転！

④白の再悪手で再再再逆転！

△の場所が急所で、先に△に打てた方が勝つ。続く

puyogo.app/rp?kf=MzJCIyRD... pic.twitter.com/wQHWOLn3Zv

タグ：ぷよ碁

posted at 08:23:21

#ぷよ碁この棋譜も非常に面白い！

添付画像

①黒の2段バネ！危険な手ですが、実はこれでも白は全滅している。この後、激戦になる！

囲碁は危険だが成立する手を打てて激戦になるときが一番楽しいと思います。

②次が黒番なので黒の勝ちになっている。

しかし～続く

puyogo.app/rp?kf=MzJCIyRD... twitter.com/mao_prv/status... pic.twitter.com/WHtXOqoNXm

タグ：ぷよ碁

posted at 08:23:15

ナカイサヤカ @sayakatake

うさぎ林檎＠ししょーPPMP @usg_ringo

@usg_ringo 助産師会が桶屋式とか使えばいいと考えちゃったからなんだけどねえ。助産師って全員看護師のはずなんだけど、困るわねえ。

タグ：

posted at 08:20:05

@sayakatake 茨城県助産師会が紹介してる。相変わらず地雷だらけの助産院。
www.sanba-ibaraki.com/map

タグ：

posted at 08:12:40

ナカイサヤカ @sayakatake

怖くてこんな助産院行けないわ。 twitter.com/rubyeyes9/stat...

タグ：

posted at 08:06:50

コノマルスー @konomaru22

まさみつ@囲碁NFT開発 @igokyoto

この白さ、丸さ…　　伝われ！

#猫 #cats pic.twitter.com/OYLuT85Fnt

タグ： cats 猫

posted at 07:51:07

つぶや棋譜で、着手ではない黒石白石（SGFのAB、AWプロパティ）がAndroidスマホの一部でで表示されないという大問題があるようですね…

#Javascriptの何かの仕様の問題なのか🤔 twitter.com/igokyoto/statu...

タグ： Javascriptの何かの仕様の問題なのか

posted at 06:53:44

モト@PPMMPP @29silicon

東北新幹線、福島〜仙台間の復旧は４月20日前後というニュース
あらためて被害の大きさがわかる…

タグ：

posted at 06:49:28

@kuri_kurita

pic.twitter.com/x4SlyfXBuk

タグ：

posted at 06:20:18

まさみつ@囲碁NFT開発 @igokyoto

前野［いろもの物理学者］昌弘 @irobutsu

つぶや棋譜の調査お願い🙏
gokifu.net/t2.php?s=96316...

９路盤３子局で白が左上に打った局面なのですが、置石が表示されない人(androidスマホ)がいました。
白石しかみえない！という方はブラウザについて情報のコメントいただけると嬉しいです。

タグ：

posted at 05:37:38