黒木玄 Gen Kuroki(@genkuroki)/2020年02月18日

Great package I recently discovered ChemometricsTools.jl github.com/caseykneale/Ch... . Has a wealth of useful functions with quite a bit of flexibility and intuitive parameterization #JuliaLang #julia

タグ： julia JuliaLang

posted at 00:03:04

みやを @A_Utsu_tw

こういう解説は本当にためになります。
工学の徒が測度論から確率を学ぼうとすると集合論の沼にはまる(まって沈んだ)ので、、、 twitter.com/genkuroki/stat...

タグ：

posted at 00:03:11

みやを @A_Utsu_tw

まずは考え方のイントロダクションみたいなところを理解して、
こまかいことが気になったら、怖いもの見たさで沼をあさるって感じでやらないと、ついていけないです。。

タグ：

posted at 00:05:41

田嶋陽子で強く記憶に残っている場面。TVタックルだったと思うが、橋下徹が田嶋に「なんで（公共施設の）女性センターしかないんですか？男女平等だったら男性センターもないとおかしいでしょう！」と食ってかかったことがあった。
↓

タグ：

posted at 00:10:26

田嶋「そんなことも分からないの!? この社会そのものが男性センターだからだよ！」
橋下は100%のキョトン顔で、わざとはぐらかしてるんじゃなく心の底から「何を言っているのか分からない」といった表情で、
橋下徹「は!? どこにあるんですか男性センター、どこにもないじゃないですか！」

タグ：

posted at 00:15:07

#統計訂正

❌1個抜き出す交差検証
⭕️1個抜き出し交差検証

WAICは素晴らしいアイデアだが、交差検証一般のシンプルな考え方も魅力的。

❌ダィガンマ函数
⭕️ディガンマ函数 (digamma(x))

トリガンマはtrigamma(x), polygamma(1,x)

twitter.com/genkuroki/stat...

タグ：統計

posted at 00:15:23

田嶋「だから世の中全部が男性センターなんだよ！」
橋下「いやだからどこにあるんですか！場所を教えてくださいよ！」
みたいなやり取りが続いた後、田嶋は「話になんないよ！」とスタジオを出ていった。橋下馬鹿説はあると思うが、もしかして世の中こんな奴らばかりなのかも…と思ってゾッとした。

タグ：

posted at 00:21:34

残念ながらその状況は保守というか、当時の自民党政権も同じなんですよ。でなければ「金丸訪朝団」なんてあり得なかったですからね。 twitter.com/eikouarucroati...

タグ：

posted at 00:46:41

さらに言えば、自民党政権も拉致問題を真剣に考えていなかったからこそ、2002年の小泉訪朝があった。日本政府もメディアもまさかあの段階で北朝鮮から拉致問題に関わる「自白」があると思っていませんでしたからね。保守的な人々が拉致問題を真剣に考えいてた、というのは「都市伝説」です。

タグ：

posted at 00:49:22

Viral B. Shah @Viral_B_Shah

Franklin.jl is a beautiful static site generator written in #julialang with many nice features - math rendering, Julia code evaluation, Literate.jl integration, and much more, developed by @t_lienart

franklinjl.org pic.twitter.com/o98YCyQW8F

タグ： julialang

posted at 00:50:20

この問題で「左右の別」があまり意味がないことを典型的に示すのは、1988年に初めて日本政府が公式に拉致問題に言及する発言を引き出したのが日本共産党の質問だったこと。社会党と違って共産党は北朝鮮の政権政党、朝鮮労働党と対立関係だったからですね。大韓航空機撃墜事件に絡んでの話です。

タグ：

posted at 00:55:43

2年後の金丸訪朝に典型的に表れた様に、当時の自民党は社会党と並んで、北朝鮮との国交樹立に積極的だった。だからこそ、金丸自民党副総裁と田辺社会党副委員長が共に平壌に行くことになる訳です。

タグ：

posted at 00:58:02

なお、1980年代から1990年代の初頭の自民党は北朝鮮のみならず、中国に対しても極めて融和的でした。だからこそ天安門事件に際しても他の先進国に先駆けて、制裁緩和に動く訳です。30年以上も前の政治状況を現在の状況から類推しても何の意味もありません。

タグ：

posted at 01:03:34

例えば、1980年代半ばの段階では、自民党と公明党は中国と韓国双方の関係を重視している一方、北朝鮮との国交正常化を視野に入れていて、社会党は韓国と対立（所属議員は入国すらできなかった）する一方、中国と北朝鮮とは友好関係にあり、共産党は中国とも北朝鮮とも韓国とも対立していた。

タグ：

posted at 01:08:14

基本的に当時の日本共産党はベトナムとルーマニアしか友達がいない状態です。覚えておきましょう。

タグ：

posted at 01:15:32

#Julialang

Cxx.jl
github.com/JuliaInterop/C... pic.twitter.com/WiIXnV3ohz

タグ： Julialang

posted at 02:21:23

久しぶりに動かしたけれどちゃんと動いていてビビる。 twitter.com/MathSorcerer/s...

タグ：

posted at 02:21:59

#統計見逃していた新解説(必読)

watanabe-www.math.dis.titech.ac.jp/users/swatanab...
最適なパラメータと特異点が一致していないとき
渡辺澄夫
q_and_a2020217.pdf

私の言葉遣いでは「測度零集合であっても、モデルの複雑さとサンプルサイズで決まる『解像度』によっては現実的に意味を持つことが大いにあり得る」という話。

タグ：統計

posted at 02:31:08

#統計 WAICを使いたい人にとって次のスライドは必読だと思う。

watanabe-www.math.dis.titech.ac.jp/users/swatanab...
ベイズ統計基礎　講義（修士課程)
↓
watanabe-www.math.dis.titech.ac.jp/users/swatanab...
目標条件つき独立の問題

中身が何も見えないと読む人が減るとおもうので、スクショを一枚だけ引用。

このスライドの方向は応用上とても重要。 pic.twitter.com/HEqawUBixY

タグ：統計

posted at 02:43:11

l_ppp @ppp3141592ppp

@genkuroki 安定のごまさんで見た活用例です。
賢いですよね。

twitter.com/MathSorcerer/s...

タグ：

posted at 06:42:02

centerou @mtanaka160

@takusansu @hgn_no_otaku @genkuroki 考える力をつけたい場合、都合の良くないものになっているんですね。

タグ：

posted at 06:44:41

Tadashi Nagasawa@Arm @Sonopapa

英The Economist紙による、日本の経済失速に関する論評。
安倍政権の失策の結果と、これ以上ないくらい明確に指摘。

Fiscal, natural, viral – Japan’s GDP shrinks dramatically after a tax rise and a typhoon www.economist.com/finance-and-ec...

タグ：

posted at 07:03:59

質問者2 @shinchanchi

田中秀臣さんの記事。安倍晋三総理の国会答弁は緊縮太郎飴のように財務官僚の仲間たちと同じレベル…「君子豹変す」大胆な経済政策を期待する / 1件のコメント b.hatena.ne.jp/entry/s/ironna... “「内需総崩れ」安倍首相の楽観シナリオを壊すのは「桜」ではない” (2 users) htn.to/3rAujDKUah

タグ：

posted at 07:41:14

非公開

タグ：

posted at xx:xx:xx

吉田弘幸 @y__hiroyuki

ことちとパーティーチャーの Game s @kotochi_pa

早稲田の准教授は，実は共通テストに反対でワザとやったのかな。と思うくらい，やり方がユルユル。

タグ：

posted at 08:16:24

なんでこれ式が間違いなの？( ･᷄д･᷅ ）
7✕9、9✕7どってでも良くない？
今の小学校ダメなの？

式を逆にする事で、✕にされて答案が返ってくる事が多々ある！

注目しなければいけないポイントってそこじゃないでしょ？

#ここが変だよ小学校
 #掛け算 pic.twitter.com/P2XHODa2Bd

タグ：ここが変だよ小学校掛け算

posted at 08:17:52

@kuri_kurita

基本的には同意なのですが、あの言葉はもう少し慎重に選んで欲しかったです。

こういう連中にエサを与えちゃったからなぁ… （まだまだ他にも酷いツイートが沢山あるぞ！😣） twitter.com/ebi_kosuke/sta... pic.twitter.com/3HT0VyoDYA

タグ：

posted at 08:44:11

Hiroo Yamagata @hiyori13

買って読んだ。明解です。 twitter.com/rikomrnk/statu...

タグ：

posted at 08:45:48

Hiroo Yamagata @hiyori13

@hidekih 一本だけなら480円。中身は、目新しくはない。そんなに死んでない、子供はなんだか不思議と大丈夫、まだすべての感染源は中国にたどれるからパンデミックとか言って煽るな、クルーズ船で日本はよくやっている、不明点まだ多し、落ち着け、という内容。ソース確認したければお金出す。納得すれば不要。

タグ：

posted at 08:53:40

optical_frog @optical_frog

アオイ模型土-東エ38b @aoi_mokei

フィナンシャルタイムズ「専門的な意味での景気後退に入りつつある日本」英文）www.ft.com/content/561ac7...

タグ：

posted at 08:53:42

ラブライブを応援している地方議員を恫喝してきた国会議員、表現規制推進＋反ワクチンという左翼の悪いところを煮詰めたような人だった

タグ：

posted at 09:13:17

政府側から働きかけがあったという話でしょ。言い換えれば自民党自らは、自党の立場があって提起できなかった。だから北朝鮮と関係が悪かった共産党を使った。そういう話でしょ。まさに「左右の別」がこの問題で意味がなかったことを示すエピソードです。 twitter.com/oyakochoko/sta...

タグ：

posted at 09:48:19

当時、リアルタイムで新聞よんでたら、共産党がソ連や中国、さらには北朝鮮はもちろん、ユーゴスラビアやアルバニア等とも円滑な関係になかったのは常識です。中学生だった僕でも知ってましたので。

タグ：

posted at 09:54:32

羽藤由美 @KITspeakee

つまり、対策問題集を作成した研究者はいまだに共通テストの作問グループのリーダーをしている。今後も、作問に携わる「例年約600人の大学教員」は「類推できるような情報」を漏らさなければ対策問題集の作成に関わってよいということですね。（こんな風に逃げるからこんな目に合ってるのに...絶望的） twitter.com/ut_med_lecture...

タグ：

posted at 10:05:18

kcrt @kcrt

@bicycle1885 買います
kindle版もでます？

タグ：

posted at 10:17:03

#統計続き。添付画像は

watanabe-www.math.dis.titech.ac.jp/users/swatanab...

より。例えば「状態空間モデル」の類では【そもそも汎化誤差は何で測るのでしょうか】(これは「そもそも予測精度を何で測るのでしょうか？」と同じ意味だと思ってよい！）という状況。

「状態空間モデル」ユーザー達はどうしているのでしょうか？ pic.twitter.com/veTkaqkiqV

タグ：統計

posted at 10:19:08

#統計「観測データを生成するという制限による条件付き確率分布を解とする方法」(所謂ベイズ統計の方法)では、条件付き確率分布を作るときの分母がないと困るので、常に「周辺尤度」=「分配函数」やその対数の-1倍の「自由エネルギー」はあると思ってよい。

相対的に「汎化誤差」の概念は難しい。

タグ：統計

posted at 10:23:33

#統計統計学入門の設定である独立同分布確率変数列の統計学でも、その辺の相対的な難易度の違いがあると思う。

タグ：統計

posted at 10:34:37

タグ：統計

posted at 10:34:40

あ〜る菊池誠(反緊縮)公式 @kikumaco

#統計モデル内の仮想世界での確率法則による現実のデータを生成している確率法則のシミュレーションの誤差の大きさは、自由エネルギーの標本分布(未知)に関する期待値で比較可能。標本分布(未知)に関する期待値は知りようがないので、自由エネルギーの実現値の比較で代用するという処方箋が使われる。

タグ：統計

posted at 10:42:05

失政です。民主党政権は愚かな経済政策によって政権を追われました。その後を受けた安倍政権は経済状況を回復させかけましたが、道半ばにして路線を転換し、自分たちがやってきたことを元の木阿弥にしました。日本の経済政策は失政続きです twitter.com/sonopapa/statu...

タグ：

posted at 10:45:40

#統計期待値の推定値として実現値を採用することについては、「それでうまく行くのか？」と疑問を持つべき。その疑問は自由エネルギーの漸近挙動を理解するまで解消されない。

自由エネルギーの漸近挙動が特別な形をしているおかげで、その期待値の推定値として実現値を採用することが実用的になる。

タグ：統計

posted at 10:47:57

Kiyoshi Mizumaru @kmizumar

この本を宣伝しろという電波を受信した気がする：
www.oreilly.co.jp/books/97848731...

Juliaプログラミングクックブック
――言語仕様からデータ分析、機械学習、数値計算まで
Bogumił Kamiński、Przemysław Szufel　著、中田秀基　訳 twitter.com/hidemotoNakada...

タグ：

posted at 10:50:45

Kiyoshi Mizumaru @kmizumar

Juliaはだいぶ昔にちょいと遊んで好きな部類の言語だわって印象は持ってる。面接んときにJuliaの話題振ってくいついてきたのぼくだけでしたよって @chezou に言われたような気がするけど業務で使う使わないは別にして気になったら一通り試して遊ぶよね普通。

タグ：

posted at 10:53:42

maru @Maru11992524D

過不足なく、菊池先生が指摘される通りだと思います。

『安倍政権は経済状況を回復させかけましたが』←この表現を認められない方が多いように感じます。功績を認めてこそ活路が開けるんじゃないのかなぁと twitter.com/kikumaco/statu...

タグ：

posted at 11:04:54

共立出版アリがと蟻 @1738310

本書は、新しいプログラミング言語でアリます「Julia」を学ぶための格好の書です！コロナ社さんのイチ押し『1から始める Juliaプログラミング』（進藤裕之・佐藤建太著）のご紹介です。www.coronasha.co.jp/np/isbn/978433...

タグ：

posted at 11:13:47

@rica_matsumoto3

電車の中で子供がギャンギャン泣いてテンパっいるお母さんが、絵本をだしたりオモチャをだしたり、あの手この手で一生懸命あやしてた。いつもちいちゃなピカチュウのキーホルダーを持ち歩いているようにしてる自分。それを渡しに行って、「泣かないで〜ホラ見て見て光るピカチュウ ♪」
そしたら

タグ：

posted at 11:17:39

桑満おさむ @kuwamitsuosamu

コロナウイルスが話題の今、「病原体」とは何かを再確認しよう！ www.cl20.jp/portal/pathogen/

わかりやすい。

タグ：

posted at 11:21:50

@rica_matsumoto3

パッと泣きやんだ。
「ピカチュウひかるね〜♪ありがとうございます」とお母さんが笑顔になってくれた。
多分子供は知らない人に急に話かけられてビックリして緊張のあまり泣きやんだと思われる。。😅笑笑
電車降りる時「子育て頑張って〜♪」とお別れしました。とさ…

タグ：

posted at 11:23:45

#統計

ベイズ統計に関するそれらの【とにかくわかりやすかった説明】は(厳密には具体的にどういう説明なのか特定する必要がありますが)、ことごとく間違っていたと判断する方が無難だと思いました。本当はすっきりしようがない話のはずなので。

「騙された！」と感じることは結構大事なことかも。 twitter.com/psych91415647/...

タグ：統計

posted at 11:25:39

@kcrt 電子版が出るとは聞いていますが，紙版より出版が遅くなるそうです。また，他のコロナ社の書籍を見る限り，Kindle版としては出ないのではないかと思います。

タグ：

posted at 11:28:18

#統計関連

モデルの何らかの適切な意味での精度を相対的に客観的に比較する方法の話。

モデルをテキトーに与えてデータにフィットするように調節して終わりだと、さすがにまずい。 twitter.com/genkuroki/stat...

タグ：統計

posted at 11:33:19

ベイズ統計たん @Bayesian_tan

ベイズの穴 arxiv.org/abs/2002.06467

タグ：

posted at 11:40:32

#統計自由エネルギー(または周辺尤度)をモデル選択に使うときには、本当は自由エネルギー自体を使うのではなく、標本分布に関するその期待値を使いたいのだが、それは不可能なので実現値の方を使う。

この辺について全然理解していない人達は非常に多いと思う。KL情報量のSanovの定理が基礎になる。

タグ：統計

posted at 11:42:08

このツイート、たくさんRTされていて「男性センターもあっていい」というリプも来るけど、一番の肝は絶望的な話の通じなさ、前提条件の無自覚さにあるわけで。この断絶が橋下との間だけでなく、男性社会そのものとの間にも存在していたのなら、めちゃ怖いなと。番組は20年前だけど、今はどうだろうか。

タグ：

posted at 11:44:20

原田　実 @gishigaku

@aoi_mokei @ROBOmito というか、それってもともとは右翼のような＞表現規制推進＋反ワクチン
twitter.com/aoi_mokei/stat...

タグ：

posted at 13:01:29

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

covthapapa @ArisatoMINAMI

@antiMulti @gishigaku @sinkurou @aoi_mokei 逆に言うと、左右とは別次元の対立軸が、確固として存在する、ということですね。

タグ：

posted at 13:20:22

@0QGIUhndmSpa8I2

渡辺澄夫先生のベイズ統計とか、このPDFとかが面白いです：
watanabe-www.math.dis.titech.ac.jp/users/swatanab...

タグ：

posted at 13:27:06

非公開

タグ：

posted at xx:xx:xx

Ishida the Brain Dam @tbs_i

この辺は「中学校くらいで『ムー』にはまって後からそこから抜ける」という経験をした方が偽科学に騙されにくい、みたいな話にならないかなと思ってみてる。 twitter.com/genkuroki/stat...

タグ：

posted at 14:13:23

非公開

タグ：

posted at xx:xx:xx

共立出版アリがと蟻 @1738310

左の本が、東京大学出版会さんの発行図書の中で、2019年にジュンク堂書店（丸善含む）で一番売れた本で、右の本が、共立出版の発行図書の中で、2019年に一番売れた本です。沢山のお客さまにお買い求めいただきました。誠にアリがとうございました。【撮影協力：ジュンク堂書店池袋本店7階】 pic.twitter.com/YDi8TPfX1K

タグ：

posted at 15:01:01

Nishiyama @Enutoesu

深層学習はガウス過程 - 作って遊ぶ機械学習。 machine-learning.hatenablog.com/entry/2018/01/...

タグ：

posted at 15:14:21

#統計

⭕️統計力学だけではなく、ベイズ統計にも分配函数や自由エネルギーが出て来る。

⭕️統計力学だけではなく、統計学でもエントロピー(正確には相対エントロピー)に当たる概念が基本的な役目を果たす。

❌統計力学におけるカノニカル分布をShannonエントロピーで基礎付けられる。

最後のは誤り。

タグ：統計

posted at 15:34:37

#統計例えば、ℝ上の分布p(x)で、∫H(x)p(x)dx = c = (一定)の条件のもとでShannonエントロピー ∫p(x) log p(x) dx を最大化するものは、p(x) = Z⁻¹(β) exp(-βH(x)) (βはcから決まる)の形になる。これで「Shannonエントロピーで統計力学のカノニカル分布を基礎付けられる！」と思うのは誤り。

タグ：統計

posted at 15:34:39

#統計一般化。事前分布φ(x)のもとで、分布p(x)で、

∫H(x)p(x)dx = c = (一定)

の条件のもとでKullback-Leibler情報量

D(p||φ) = ∫p(x) log(p(x)/φ(x)) dx

を最小化するものは、

p(x) = Z⁻¹(β)exp(-βH(x))φ(x) (βはcから決まる)

の形になる。統計学での指数型分布族はこのようn理解できる。

タグ：統計

posted at 15:34:40

#統計 KL情報量D(p||φ)は、各々が分布φ(x)に従う独立同分布確率変数列の経験分布として、分布p(x)に近い分布がどれだけ出てき難いかを意味する指標になっている(Sanovの定理)。

KL情報量最小化は「最も出てき易い分布を見つけること」に対応している。

タグ：統計

posted at 15:34:40

#統計物理的には、独立同分布の仮定は強過ぎるので、Kullback-Leibler情報量によるカノニカル分布の特徴付けは物理的には不完全な特徴付けになる。

Sanovの定理に類似のある種に漸近挙動を物理的に仮定するというのが、おそらく標準的な統計力学でのカノニカル分布の導出法である。

タグ：統計

posted at 15:34:41

#統計物理的厳密性にこだわればそういうことになるのだが、KL情報量を使った独立同分布確率変数列版のカノニカル分布として、統計学によく現れる指数型分布族を導出する議論は分かりやすい。

指数型分布族のそういう特徴付けは、統計力学のトイモデルになっていると思ってよい。

タグ：統計

posted at 15:34:42

#統計独立同分布確率変数列 X_1, X_2, …, X_n は熱浴のトイモデルだとみなせ、そのうちの1つ、例えば X_1 は熱浴に接している系のトイモデルだとみなされる。

独立同分布確率変数列の数学は比較的易しいので、易しいトイモデルが欲しい人はそのように考えればよい。

タグ：統計

posted at 15:34:43

#統計「熱浴」の概念はおそらく統計力学を超えて数学全般において重要な意味を持つ可能性のある(曖昧だが)基本的な概念だと思われる。

タグ：統計

posted at 15:34:43

#統計「熱浴を含む系の全体」X_1, …, X_n の分布を「エネルギー保存則」

(1/n)Σ_{k=1}^n H(X_k) = c = (一定) (正確には≦または≧とする)

によって条件付き確率分布に制限してやると、n→∞で「熱浴に接している系」X_1の分布は const. exp(-βH(x))φ(x) 型の分布に近付く(Sanovの定理の帰結)。 twitter.com/genkuroki/stat...

タグ：統計

posted at 15:41:52

#統計より詳しい解説↓

genkuroki.github.io/documents/2016...
Kullback-Leibler 情報量と Sanov の定理

タグ：統計

posted at 15:43:39

#統計統計学にKullback-Leibler情報量が自然に出て来る理由を理解するためには、「情報量」という用語の字面的意味に寄せた理解を目指すと損で、素直にSanovの定理を学ぶ方がよいと思う。

モデルの分布p*(x)による推測先の分布q(x)のシミュレーションの誤差がKL情報量D(q||p*)になる(Sanovの定理)。 twitter.com/yosuke_hash/st...

タグ：統計

posted at 15:58:12

#統計未知の確率分布q(x)が生成したサンプル(データ)から、未知のq(x)を近似的にシミュレートできそうな分布p*(x)を作ることは、統計学や機械学習の基本的な目標の1つなので、Sanovの定理によって、シミュレーションの誤差としてKullback-Leibler情報量が必然的に必要になるわけです。

タグ：統計

posted at 16:01:12

#統計ぶっちゃけ、Kullback-Leibler情報量を「情報量」として特徴付ける自然そうに条件をどんなに説明しても、それだけだと「どうして統計学でそれが必須になるのか」は分からないと思う。

統計学に適合するKL情報量の理解の仕方はSanovの定理だと思う。統計力学との関係もSanovの定理が本質的。

タグ：統計

posted at 16:05:03

#統計統計学入門で必要な確率論の三種の神器は

* 大数の法則
* 中心極限定理
* Kullback-Leibler情報量のSanovの定理

の3つ。Sanovの定理の易しい解説は少ない。

タグ：統計

posted at 16:06:58

岩永直子 Naoko Iwanaga『言 @nonbeepanda

#統計赤池弘次さんも1980年の論説で、統計学でエントロピーの類似物が自然に出て来る理由を本質的にSanovの定理を経由して説明しています。

Sanovの定理は、統計学にKullback-Leibler情報量が出て来る理由の「由緒正しい説明の仕方」だと思います。

twitter.com/genkuroki/stat...

タグ：統計

posted at 16:11:16

男子よ…。 pic.twitter.com/TkmH4pUC39

タグ：

posted at 16:16:24

@kado_judo0312 Juliaと4Kモニターイィゾォ

タグ：

posted at 16:46:10

Kei NAKAGAWA @keeeeei0315

あ〜る菊池誠(反緊縮)公式 @kikumaco

裾の厚い金融時系列によりフィットさせるため、ガウス過程潜在変数モデルをt過程潜在変数モデル(TPLVM)に拡張したワーキングペーパーをMAZINの内山さんと書きました。昨年SigFinで賞を頂いた内容の潜在変数バージョンです。 twitter.com/StatMLPapers/s...

タグ：

posted at 18:09:50

僕たちは、ボランティアで正しい情報を発信してくれている医療関係者を守らなくてはならないのです。絶対にです

タグ：

posted at 18:13:45

みずたま @go_tom_89

この手は詰めろです。放置すると…

７手詰

#おやすみ詰将棋 pic.twitter.com/8UDSvAe08N

タグ：おやすみ詰将棋

posted at 18:24:24

log(hoge(x))ではなく、loghoge(x) を使うだけではなく、他にも注意するべきことがあります。

統計学では x log y 型の計算がよく出て来ます。これは

xlogy(x, y) = iszero(x) ? x : x*log(y)

のように、xが0の場合を例外扱いしておかないとはまる場合がある。

twitter.com/genkuroki/stat...

タグ：

posted at 19:11:26

safemul(x, y) = iszero(x) ? x : x*y

としておいて、

safemul(x, log(y))

を使うと失敗します。なぜならば、x = 0 であっても、y < 0 のとき log(y) はエラーになるので、safemul(x, log(y)) はエラーになります。

xlogy(x, y) = iszero(x) ? x : x*log(y)

ならば大丈夫。

タグ：

posted at 19:13:50

他にも統計学の計算では logsumexp の実装で注意しないと失敗する。完璧な実装を目指したい人は以下のリンク先を参照。log(1+x) を直接計算してくれる函数があるなら、それを使う。

統計がらみの計算では、細かいところで注意を要することが結構あります。

twitter.com/genkuroki/stat...

タグ：

posted at 19:17:32

積分定数 @sekibunnteisuu

＞今の小学校ダメなの？

その通り。駄目なのは答案ではなく教え方の方。 twitter.com/kotochi_pa/sta...

タグ：

posted at 19:32:19

うちの家系は7代前からPythonパッケージのインストールを生業としていて、祖父と父は人間国宝にも指定されていた。

タグ：

posted at 20:15:37

TrueRoad @trueroad_jp

いくつか RT しましたが、拙作の原ノ味フォントが CTAN に収録 (haranoaji, haranoaji-extra) され、すぐに TeX Live にも収録されました。pTeX マップファイルも原ノ味に対応するものが入ったので、最新の TeX Live であれば、 pTeX, LuaTeX, XeTeX いずれも簡単に使えるようになったと思います。

タグ：

posted at 20:24:56

祖父のインストールしたパッケージは本当に見事で、子供の頃はそんな祖父に憧れて自分もインストール職人になりたいと本気で思った時期もあったが、中学2年のときに祖父は他界し、父への反発からその後はJuliaやRustのパッケージインストールに夢中になっていった。

タグ：

posted at 20:25:21

@bicycle1885 思ってたよりもPython歴長くてびっくりしました！

タグ：

posted at 20:35:52

@MathSorcerer みなさん驚かれるようです

タグ：

posted at 20:46:38

#統計 #Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
Bernoulli分布モデル

にLOOCV(1個抜き出し交差検証)とWBICを追加した。すでにAIC、WAIC、BIC、自由エネルギーは含まれていた。手書きの計算ノートも追加(添付画像)。

勉強用にどうぞ！ pic.twitter.com/d5lJYH5b3l

タグ： Julia言語統計

posted at 21:38:29

#統計ベルヌイ分布モデルの

AIC、WAIC、LOOCV、BIC、WBIC、自由エネルギー

の計算は、ベータ函数やガンマ函数の基本性質と{ディ、トリ}ガンマ函数の定義を知っていれば、正規分布モデルの場合よりずっと楽です。

情報量規準の最も簡単な例として、教育的だと思いました。

タグ：統計

posted at 21:43:57

Hidekazu Shiozawa @shiozawa_h

float f,b;void setup(){size(480,480,P3D);colorMode(HSB,3);noStroke();}
void draw(){clear();perspective(1,1,390,440);translate(240,240);
pointLight(0,0,3,0,0,200);rotateY(f+=.01);rotate(PI/4);
for(b=1;b<9;){fill(sin(b/4-f*2)+2,3,3);rotateY(f/4);box(b++*60);}}
#つぶやきProcessing pic.twitter.com/G2522Zp0FR

タグ：つぶやきProcessing

posted at 21:49:17

EARLの医学ツイート @EARL_med_tw

#統計 #Julia言語ベルヌイ分布モデルのノートに

AIC、WAIC、BIC、自由エネルギー

だけではなく、

LOOCV、WBIC

を付け加えました。

twitter.com/genkuroki/stat...

タグ： Julia言語統計

posted at 21:49:54

検査の種類も様々な感染症の経過も鑑別も分からないなら黙っててもらえますか？いくら専門外でもあなたクラスがそういうこと言うのはリスコミぶち壊しでむしろ外来診療を壊滅させるんですよ twitter.com/hatoyamayukio/...

タグ：

posted at 21:51:52

堀　成美 @narumita

感染しても症状を自覚しないほどの軽症者、病院に行こうと思わないレベルの軽症者がたくさんいるタイプの感染症では、地域流行してから「検査して見つけで隔離だ！」というアイデアじたいがズレてます。
高齢者や病気の人が重症化して苦しくならないよう、医療機関に負担がかからないことが重要です。 twitter.com/hatoyamayukio/...

タグ：

posted at 21:53:32

残り6件のツイートを見る

歩行者は右側通行 @JikanBae

職場トータルでMicrosoft Officeに契約してくれるより、みんなでLibreOffice使いましょうよってなるほうが嬉しいんだけどな。（僕はワープロとプレゼンにはLaTeXを、計算にはRとJuliaを使ってるけど）

タグ：

posted at 22:00:27

Yuki Nagai @cometscome_phys

このノートの34ページのn回ゲームしてその平均値がUを下回るまで何度もn回ゲームすると最終的にはカノニカル分布が出るって話、Juliaで実装してみたけどnが増えるにつれて平均値が期待値に近づくから数値計算的にはカノニカル分布が出ない。自分の理解が間違ってるかな
genkuroki.github.io/documents/2016...

タグ：

posted at 22:28:14

Yuki Nagai @cometscome_phys

何度戻ってもいつも平均値が期待値に近くなってて求めるUにならずまた戻ることになり、全然終わらなかった

タグ：

posted at 22:29:08

Yuki Nagai @cometscome_phys