黒木玄 Gen Kuroki(@genkuroki)/2020年05月20日

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

本日の作品。
軌跡を残しながら端でバウンドするムーバーを描画している部分が回転しているように見える見せ方で...。
#p5js pic.twitter.com/XrFt4pAptG

タグ： p5js

posted at 00:09:36

その話はビデオでも講義ノートでも演義の解答例でも何度も説明したと思うぞ twitter.com/monokage9999/s...

タグ：

posted at 00:14:07

ロボ太 @kaityo256

小学校でシャボン玉を作る授業があった。みんながシャボン玉で遊びだすなか、黙々と石鹸を定規で削って水に溶かす奴がいた。授業時間も終わる頃、皆が見守る中、満を持して超特濃シャボン液でシャボン玉を作ろうとした彼だが、液が濃すぎてシャボン玉が全く作れなかったの、なんか示唆に富むよね……

タグ：

posted at 00:15:21

#数楽「対数を取ってから比較すれば、相対誤差は小さい」という話。

log (N+1)! = log N! + log(N+1)

で、Nが巨大なら、右辺の第2項は第1項より大幅に小さくなる。

N! ≈ N^N の類の大胆な「近似」も同じ理屈で納得できる。

現実に出会う量は対数の側により直結している場合が少なくない。 twitter.com/monokage9999/s...

タグ：数楽

posted at 00:29:54

#数楽確率や場合の数およびそれらに類似の量は対数を取って扱うことが非常に多い。

タグ：数楽

posted at 00:31:22

Re:RT #Julia言語の偉い人から直接コメントをもらう最も簡単な方法の例。「Juliaでやってみたけど遅かったです」と書く。「さらにPythonより遅い」と書くと効果的。

「みんなでコードの問題点を見つけて、Juliaのコードの最適化を行おうぜ！」みたいな話になって行く(笑)

タグ： Julia言語

posted at 00:46:38

むろ @nii_tsuku

@narita_yusuke 自分もプロスペクト理論派閥に含まれるのかもですが、この実験は確率加重関数のノンパラメトリックな検証であり、損失回避「自体」についての検証は同程度には強くないというのが自分の理解です。もちろん、多くの論文に使用されている確率加重関数が検証できなかったこと自体、非常に重い知見ですが。

タグ：

posted at 02:42:58

成田悠輔 @narita_yusuke

@nii_tsuku 勉強になります。そうですよね、それが理由で論文の導入もだいぶ入り組んだ説明になってる印象です。

タグ：

posted at 02:51:52

The Unexplained @Unexplained

Children in Hokkaido, Japan are taught from a young age: NEVER STEP IN PUDDLES pic.twitter.com/HpezzbArsP

タグ：

posted at 06:34:36

Dr.こば @KobaKobauro

こういう意見を見るといつもワンパンマンのこのシーンを思い出します twitter.com/lullymiura/sta... pic.twitter.com/SPythmJnK1

タグ：

posted at 07:46:05

Michael P.H. Stumpf @theosysbio

This looks very exciting: open source computer algebra using #JuliaLang: OSCAR Computer Algebra System

#computeralgebra is one area where I still require commercial software, despite my love for ⁦@SymPy⁩.
OSCAR could ring in change. oscar.computeralgebra.de

タグ： computeralgebra JuliaLang

posted at 08:04:47

あおじるPPPP @kale_aojiru

給付や助成は施しではないので感謝を強要される謂れはないし、その金は税金でもない。 twitter.com/aac022419/stat...

タグ：

posted at 08:19:19

#統計増強版

超一般論 ↔ 俗な呼び方

観測可能な変数で実際に値が観測されたもの ↔ データ

観測可能な変数で値が得られていないもの ↔ (予測の対象)

それら以外の確率変数 ↔ パラメータ

確率変数でない変数 ↔ パラメータもしくはハイパーパラメータ

タグ：統計

posted at 08:41:44

#統計超一般論：p(x_1,…,x_N,w|z)はx_j, wについては確率密度函数でzについては単なる函数であるとする。x_j, wは確率変数とみなせ、zはみなせない。変数x_j達は観測可能変数扱い。

x_1,…,x_n, n≦Nの観測値X_1,…,X_n (データ)が得られたとき、x_{n+1},…,x_N,wを予測・推測したいならば～続く

タグ：統計

posted at 08:41:46

#統計続き～、モデル内で観測値と同じ値が生成されたという条件で制限した条件付き確率分布

p*(x_{n+1},…,x_N,w|z)
=p(X_1,…,X_n, x_{n+1},…,x_N,w|z)/Z(z),

Z(z)=∫dx_{n+1}…∫dx_N∫dw p(X_1,…,X_n, x_{n+1},…,x_N,w|z)

を考え、Z(z)を最大化するように(ハイパー)パラメータzを調節する。

タグ：統計

posted at 08:41:46

舞台裏的な話。この論文で数学的事実を事実だけでもはっきり書くことになったのは「最尤法は事後分布を正規分布で近似できないとき」と書いたら「事後分布はベイズ法の言葉で最尤法とは哲学が違う。意味が分からない」というレビュアに応答するためでした。peerj.com/articles/9036/

タグ：

posted at 08:44:09

タグ：統計

posted at 08:54:47

哲学が違う、つまり、解析者はまず哲学が違うものをどちらか選択して「正しい」推測をするという考えではなく、我々のアプローチが選択→結果の吟味の型になっていることを明示する必要が生じたのです。なので考察では基礎的枠組みから記述することにしました。これでレビュアは納得されたみたいです。

タグ：

posted at 08:55:00

ぼんてんぴょん(Bontenpøn) @y_bonten

こういうやりとりは自分の不勉強さもあって大変時間がかかりましたが、貴重ですね。林さんの参画を経てDiscussionはじめ大幅に改訂し出版までこぎつけられました。

タグ：

posted at 09:02:11

「猛勉強した」より「もう勉強した」のほうが強そう。勉強済かよ。

タグ：

posted at 09:03:23

#統計

例：p(x_1,…,x_N,w|z)=p(x_1|z)…p(x_N|z)φ(z)の場合から正則化された最尤法(MAP法)が出る。(注意：パラメータzは確率変数扱いされていない)

例：p(x,y,a,b|z)=p(a|z)p(x|a)p(b|a)p(y|b)は有効グラフ

z→a → b
　 ↓　↓
　 x 　 y

で記述されるモデルになっている。

タグ：統計

posted at 09:05:44

#統計モデルp(x_1,…,x_N,w|z)でx_j達は観測可能な変数で、実際にx_1,…,x_n, n≦Nの観測値が得られたとき、他の変数の予測をそのモデル内でするには、そのモデル内でも同じ値が観測されたという設定を採用することが自然だろう。

そのように考えた瞬間に条件付き確率分布を考えていることになる。

タグ：統計

posted at 09:10:26

#統計要するに、確率モデルを使った予測・推測では、最尤法やベイズ統計などの「俗な呼び名」と無関係に、現実での観測値と同じ値がモデル内でも観測されたという条件で制限した条件付き確率分布を考えることは必然的であるということである。(実際、最尤法の予測分布もそのようにして出せる。)

タグ：統計

posted at 09:14:28

#統計「条件付き確率分布の概念が実は神だった」という話。

最尤法の予測分布もモデル内条件付き確率分布の一種。

ガウス過程回帰は条件付き確率分布を易しく作れる場合。

統計学でMCMC法は主に条件付き確率分布のサンプルを作るために使われている。

とにかく、条件付き確率分布の構成が神。

タグ：統計

posted at 09:19:28

Daisuke Tano @tanosensei

感謝する必要はまったくない。バンバン要求するべし。 twitter.com/aac022419/stat...

タグ：

posted at 09:24:23

【ディープラーニングは「最小二乗法」】【全国民PCR検査推進】とかだけではない！(笑)

ytb.hatenablog.com/entry/20071021...
小林慶一郎氏曰く【「リスクに対処するための金融商品自身が新たなリスクになる」という論文のロジックは～「ゲーデルの不完全性定理」と関係しているのではないか、と私は考えてきた】

タグ：

posted at 09:53:18

ToshijiKawagoe @ToshijiKawagoe

累積的プロスペクト理論は、そもそもカーネマン＆トヴェルスキーの確率重み付け関数自体に問題があってその後他人が修正版を出していますし、また損失回避性の測度も人によって色々定義があって、理論としては穴だらけというのが本当のところですよね＾＾； twitter.com/narita_yusuke/...

タグ：

posted at 09:55:39

経済論壇でのシバキアゲ言説流布で御活躍するだけではなく、多くの分野において噴飯物のトンデモ発言を定期的に繰り返して笑いモノになって来た人物を政府が重用し、その最近の発言に共産党トップまで共感を示す。

これが現実。ちと苦し過ぎ。

タグ：

posted at 09:57:17

posi_posi @posi_posi8

凄い(凄い) twitter.com/genkuroki/stat...

タグ：

posted at 09:58:42

ToshijiKawagoe @ToshijiKawagoe

Masa Yamamoto予測誤差が大き @mshero_y

こうしたプロスペクト理論の問題点についても、拙著『意思決定の科学』（講談社ブルーバックス、９月刊行予定）に書きました。入門書ですが、結構つっこんだことを書いていますので、お楽しみに！ twitter.com/ToshijiKawagoe...

タグ：

posted at 09:59:03

普通の常識があれば… twitter.com/genkuroki/stat...

タグ：

posted at 10:02:58

成田悠輔 @narita_yusuke

@ToshijiKawagoe 穴だらけの方が人を惹きつけるのかもしれませんね笑　この論文では修正版の方を棄却してるようです。

タグ：

posted at 10:12:03

非公開

タグ：

posted at xx:xx:xx

仮装した　#カイカムリ
 twitter.com/numagasa/statu... twitter.com/numagasa/statu...

タグ：カイカムリ

posted at 11:13:26

歩行者は右側通行 @JikanBae

OokuboTact　大久保中二病中年 @OokuboTact

いつのまにか #Julia言語がThe Computer Language Benchmark GameのTier 1に入ってる。
benchmarksgame-team.pages.debian.net/benchmarksgame...

タグ： Julia言語

posted at 11:16:35

＃超算数

藤沢利喜太郎（1861年〜1933年）の小学校算数の掛け算の考えは高木貞治と同じだと思う。
黒表紙教科書は時代は、掛け算の順序固定指導はなかったと思う。

タグ：超算数

posted at 11:52:44

清　史弘 @f_sei

これ、やはり私以外の方からも指摘が的確だという意見があるので、再度あげておきます。
お時間のある方はどうぞお読みください。
こういうことを元にして高校生、大学生の将来が決められようとしています。
enomoton2011.hatenadiary.org/entry/20121001...

タグ：

posted at 11:58:58

news.yahoo.co.jp/articles/00230...
小林慶一郎氏曰く【ラインハート教授と～ロゴフ教授は～「債務比率が90％を超えると、その国の経済成長率が１％程度低下する」という傾向があることを報告】

これはExcelのコードの誤りが原因の間違いであったことは有名。経済の専門家扱いしてはいけないことが分かる。 pic.twitter.com/DIWJLwYSpQ

タグ：

posted at 13:25:32

ラインハルト&ロゴフの報告の結論は、実際にはExcelのコードの誤りに基く間違いであったという話はかなり面白いです。詳しい情報は私のツイログ検索経由で確認できるはずです。

twilog.org/genkuroki/sear...

彼らのトンデモ論文を潰した論文は↓

twitter.com/genkuroki/stat...

タグ：

posted at 13:29:37

以上の引用から分かるように、小林慶一郎氏は専門のマクロ経済学において、Excelのコードミスの結果の報告を自身の緊縮シバキアゲ言説を正当化するために引用するへまをしでかすほどの無知無能な人物。

このような人物に経済政策について相談する奴もまた日本を害する人物として扱われるべきです。

タグ：

posted at 13:32:56

Nori（情報ブロイラー） @bgnori

専門とされる分野においても無知無能な人物が社会的に重用されているケースは特に経済において重大な問題を引き起こして来たと思います。

タグ：

posted at 13:34:41

あー、そんなのもあったなぁ。同じ人だったんだ。 twitter.com/genkuroki/stat...

タグ：

posted at 13:45:36

Juliaの利用者・開発者アンケート調査があるんだけど、今年は英語に加えてスペイン語・中国語・日本語でも調査するみたい。日本語版の翻訳をした。

タグ：

posted at 14:16:09

やっぱ日本語圏のJuliaユーザ多くなってきてるのかな。

タグ：

posted at 14:25:12

#Julia言語

行列の成分がほとんど0なら

using LinearAlgebra, SparseArrays
n = 10
C = sparse(1:n-1, 2:n, 1, n, n)
↓
1 0×10 SparseMatrixCSC{Int64,Int64} with 9 stored entries:
[1, 2] = 1
[2, 3] = 1
[3, 4] = 1
...

とか？ twitter.com/misaki70533995... pic.twitter.com/EMwhB9j9qV

タグ： Julia言語

posted at 14:29:21

#Julia言語を使い出してから、「無駄にメモリを使うと遅くなること」がものすごく気になるようになってしまった。

無駄にメモリを使っても多くの場合に十分速いのですが、ときどき「これはまずい」というケースに打ち当たる。

タグ： Julia言語

posted at 14:33:57

#統計訂正

以下のリンク先の後者の例

❌有効グラフ
⭕️有向グラフ

正確にはループを持たない有向グラフ。

twitter.com/genkuroki/stat...

タグ：統計

posted at 14:55:38

#統計注意：超一般論の枠組みで計算した予測分布の優劣を観測データだけに基いて確率的に判定する方法の開発は重要。単体のモデルで予測しっぱなしだと、どこまで妥当なことをやっているのか不明になる。

安易に証拠を示さずに「交差検証をやればいい」と言っている人達の無責任な態度は要注意！

タグ：統計

posted at 14:59:50

#統計注意：以上のような確率モデリングの話と統計学入門で習う「仮説検定」の話が全然違う話に見えてしまっているとしたら誤解している。

「モデルKとそのパラメータ空間を制限して得られるモデルHの最尤法の予測精度の推定値の比較」と「帰無仮説Hの仮説検定」はもろに関係している。

タグ：統計

posted at 15:04:52

積分定数 @sekibunnteisuu

福岡県、人権啓発番組に削除要請　手紙1通のクレームで：朝日新聞デジタル www.asahi.com/articles/ASN5M...

タグ：

posted at 15:26:42

Yuki Nagai @cometscome_phys

fortranでjuliaのpmapみたいな関数を作りたい

タグ：

posted at 15:27:39

ceptree @ceptree

午後というかもう4時半

タグ：

posted at 16:28:42

@ceptree それはどうかな？

タグ：

posted at 16:30:14

ceptree @ceptree

@bicycle1885 パパパパパパパパパーパー pic.twitter.com/QL8K5LaaGK

タグ：

posted at 16:31:36

@ceptree 🥳

タグ：

posted at 16:35:05

非公開

タグ：

posted at xx:xx:xx

Tsuyoshi Miyakawa @tsuyomiyakawa

結核感染についてこういったご意見があるようなので少し説明。
人口あたり年間の結核新規感染者数の国ごとの値は長期的に安定していて、2000年の値と2018年の値の相関は0.87です。つまり国別ランクは20年近くあまり大きくはかわりません。
twitter.com/mixingale/stat...

タグ：

posted at 16:53:23

非公開

タグ：

posted at xx:xx:xx

Tsuyoshi Miyakawa @tsuyomiyakawa

また、ここで重要なのはある時点での「結核保菌率」ではなく、結核に一度でも感染したことのある人の率ですね。結核蔓延国は昔から現在までかなり一貫しているので、我々の論文で使用している指標（人口あたり年間の結核新規感染者数の20年弱の平均）でおおむね代替できると考えられます。

タグ：

posted at 16:56:16

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

ちくわ @pavlov469

終わり続けるデータサイエンスブームと始まり続けるJuliaブーム

タグ：

posted at 17:25:24

ちくわ @pavlov469

もちろん，今年はJuliaが始まると思ってますよ

タグ：

posted at 17:29:11

BuzzFeed Japan Medic @BFJMedical

【New】自民党の「HPVワクチンの積極的勧奨再開を目指す議員連盟」が開かれました。

9価ワクチン承認後の速やかな定期接種化や、啓発用リーフレットを誤解を招かない表現にすることなどを厚労省に求めました。（岩永直子@nonbeepanda）
www.buzzfeed.com/jp/naokoiwanag...

タグ：

posted at 17:39:55

#統計私もlikelihoodを「尤度」と訳したのは分かり難いと思う。

それとは別に、「尤度最大化」は「一番あり得そうな」とか「最ももっともらしい」を意味しないことは強調されて然るべき(多くの教科書の説明がひどい)。

「尤度」の定義は「モデル内で観測データと同じ値が生じる確率(密度)」。続く twitter.com/yaneuraou/stat...

タグ：統計

posted at 17:51:39

#統計続き。「尤度」の定義は「モデル内で観測データと同じ値が生じる確率(密度)」に過ぎないので、「モデルの観測データへの適合度」の指標としては適切であっても、「モデルが真の法則にどれだけ適合しているか」の指標としては不適切です。

尤度が大きい方がもっともらしさが下がる場合がある！

タグ：統計

posted at 17:51:40

#統計続き。パラメータを沢山持つモデルでパラメータをモデルがデータにぴったりフィットするように(すなわち尤度を最大化するように)調節すると、もっともらしさの欠けた精度の低い推定結果が得られます(オーバーフィッティング)。

タグ：統計

posted at 17:55:32

お！discourse.julialang.org/t/ann-commonma...

タグ：

posted at 17:58:19

#統計尤度は「もっともらしさ」を意味しないので、英語のlikelihoodという呼び名は不適切であったと言ってよいと思います。

likelihoodはミスリーディングな名付け方なので、それに「尤度」という意味不明の訳語を割り振ったのは次善の策としては有効になっていた可能性もありました。続く

タグ：統計

posted at 17:58:24

#統計続き。しかし、多くの教科書が「尤度」について正確な説明をせずに、「尤度はもっともらしさを表す」というような言い方で不適切な用語のlikelihoodの直訳になっている「もっともらしさ」という言い方を導入して台無しにしている。

タグ：統計

posted at 18:00:27

非公開

タグ：

posted at xx:xx:xx

岩田健太郎 K Iwata, MD, M @georgebest1969

#統計試案：尤度はlikelihood(尤もらしさ(もっともらしさ))の訳語だが、もっともらしさの指標ではないので、あえて意味不明に感じられる尤度(ゆうど)という用語を用いる。「尤度」の定義は「モデル内で観測データと同じ値が生成される確率(密度)」であり、「モデルの観測データへの適合度～略

タグ：統計

posted at 18:04:38

間違ってます。今の日本の現実があるおかげは、少なくとも部分的には西浦先生のおかげです。ただし、人材育成が必要なのは事実。逆に西浦先生がいなかったら、日本は終わっていたかも、思う、医療現場は本当にギリギリでした。 twitter.com/lullymiura/sta...

タグ：

posted at 18:04:45

Julia製のCommonMark実装っすよ皆さん！

タグ：

posted at 18:13:09

これらはとてもいい記事 RTs：

中学生から始める適応進化理論
lambtani.hatenablog.jp/entry/2020/05/...
中学生から始める適応進化理論2：よくある誤解
lambtani.hatenablog.jp/entry/2020/05/...

タグ：

posted at 19:03:48

子供の頃にPythonばかり書かせてJuliaを禁止すれば、その反動で大人になったときJuliaばかり書くようになるのでは🤔

タグ：

posted at 19:10:09

天むす名古屋 Temmus @temmusu_n

義務教育ではPythonを採用して下さい！

タグ：

posted at 19:12:30

天むす名古屋 Temmus @temmusu_n

@m50185 教育理論には割とひどいものもありますが、研究者の相互批判が機能している面はあります。そういう指導書をノンノンさんはお読みになったことでしょう。しかし教師が最も参考にするのは教科書の教師用指導書でなんです(そういう質問調査あり)。現場の理論的支柱は一般人の目に届きにくいところに存在。

タグ：

posted at 19:51:41

@m50185 教師用指導書の購入は一般人には不可能です。自治体の中央図書館にもないし、なぜか古本も流通しません。教科書の記述の意図が指導書を見て初めて分かることも多く、公開文書である指導要領、文科省解説、論文なんかでは、なんでこのせんせいはこんなこと言うの？という疑問に答えは得られないことも。

タグ：

posted at 20:35:43

Sogachin23 @sogachin2

今日は図書館で、「StanとRでベイズ統計モデリング　松浦健太郎著　共立出版」を借りた。ベイズ統計をRにて実装できるように勉強する。

タグ：

posted at 21:06:21

hidesan @waverider_hide

こういう仕事はせっせとやってほしい。
「9価ワクチン承認後は速やかに定期接種に」　自民党のHPV議連が訴える www.buzzfeed.com/jp/naokoiwanag... @nonbeepandaさんから

タグ：

posted at 21:06:46

RochejacMonmo @RochejacMonmo

大学入試のあり方に関する検討会議（第7回）
新井紀子氏の資料
「大学入学共通テストにおける記述式問題のあり方について」
www.mext.go.jp/content/202005...
記述式に対する提案の数IA部分。
実は具体性に乏しいこととなぜ数IAという制約を十分に考慮していないと感じる。 pic.twitter.com/Pu9TDzAju6

タグ：

posted at 22:58:42

RochejacMonmo @RochejacMonmo

『問題文をきちんと読解しないと解けない「読ませる問題」』の意味がはっきりしない。
例えば試行調査の問題は、日常との関わりや問題と必ずしも関係しない会話文などが増えた結果、そこから必要な情報を取り出す力が求められるとは言っても、問題そのものの数学的な内容が減退してしまっていないか。

タグ：

posted at 23:02:04

RochejacMonmo @RochejacMonmo