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

Another great #julialang post
biojulia.net/post/hardware/

タグ： julialang

posted at 00:14:30

カイヤン @389jan

某かるてっと，正規乱数な1次関数モデルである限りは最尤でもベイズでもほぼ同じ結果になるんかなこれ（なりそう）

タグ：

posted at 00:16:21

カイヤン @389jan

分散推定までやればいけるかな？

タグ：

posted at 00:16:55

カイヤン @389jan

んんん，どういう仮想シナリオを立てると面白いだろうこの問題．散布図が全然違うのに1次関数フィット結果がとても近い値になる例だが，そうなってしまったスタート地点からどう進めていくのが面白いのかな．
とりあえず2次項までいれたらもちろん全然違う結果になるんだけど

タグ：

posted at 00:25:57

#統計 Anscome's quartetの4つのデータセットの単純線形回帰モデル p(y|x,a,b,σ) = (1/√(2πσ²))exp(-(y-(a+bx))²/σ²) の尤度函数はほぼ同じになる。

だから、ベイズ版でも4つのデータセットはほぼ同じ推定結果を与えます。

#Julia言語
 nbviewer.jupyter.org/gist/genkuroki...

twitter.com/389jan/status/... pic.twitter.com/BjMl5ooJnI

タグ： Julia言語統計

posted at 07:01:49

#統計指数型分布族のモデルではデータの代表値だけで尤度函数が決まってしまいます。十分統計量や共役事前分布関連の話題。

#Julia言語のTuring.jlでWAICやLOOCVなどを計算したい人は

nbviewer.jupyter.org/gist/genkuroki...
TuringでWAICやLOOCVを計算する方法

を参照。

タグ： Julia言語統計

posted at 07:04:13

#Julia言語 #統計

 nbviewer.jupyter.org/gist/genkuroki...
Baysian version of the plot of Anscombe's quartet

では、Rのデータセットの1つであるanscombeをJuliaで読み込んで、Turing.jlでベイズ版単純線形回帰を行い、予測分布を計算してプロットするということをやっています。

結構、参考になるはず。

タグ： Julia言語統計

posted at 07:08:26

#統計 #数楽非常に困ったことに多くの人が「ベイズ推定」だとか言っているPCRの{陽性,陰性}的中率の話は実は小学生レベルの割合の計算に過ぎません。

実際、本質的に同じ大学入試問題を小学生にやらせてみたら、確率などの用語を割合の言葉で言い換えてあげたら、比較的容易に解けました。続く

タグ：数楽統計

posted at 07:15:26

#統計 #数楽上の方で「単純線形回帰の尤度函数」というような専門用語をバリバリ使った言い方をして説明した事柄は、本質的に中学校で習う平方完成

ax²- 2bx = a(x-b/a)² - b²/a

の多変数(もしくは行列とベクトル)への一般化に過ぎません。

義務教育レベルの数学に理解は非常に重要！

タグ：数楽統計

posted at 07:19:08

#統計 #数楽「単純線形回帰の尤度函数」の計算は次の平方完成の問題に帰着する。

問題: n次元縦ベクトルy, n×m行列X, m次元縦ベクトルβについて

||y - Xβ||² = (y - Xβ)^T (y - Xβ)

をβの函数として平方完成せよ。ただし、X^T は X の転置を表すとし、X^T X の逆行列の存在を仮定する。

タグ：数楽統計

posted at 07:27:47

おり久 @OriqMen331

鈴木エイトジャーナリスト/作家『「山上 @cult_and_fraud

おはようございます。落合、木町通り本日も11:00〜営業致します。木町通り限定持ち帰りです。 pic.twitter.com/y7IqMNwMaY

タグ：

posted at 07:31:33

この件、リュック・モンタニエが（ヒト免疫不全ウイルス／エイズウイルス／HIV発見で2008年にノーベル生理学・医学賞を共同受賞）が #新型コロナウイルス人造説＆武漢の研究所漏出説を唱えている。モンタニエは #統一教会主催のトンデモ科学者会議にも参加するなど晩節を汚す科学者の一人 twitter.com/yahoonewstopic...

タグ：新型コロナウイルス統一教会

posted at 07:33:14

#統計 #数楽

小学生レベルの単なる割合計算についてわざわざ「ベイズの定理」やら「ベイズ推定」やらの用語を持ち出して、何か特別に高尚な事柄が必要だと思ってしまうのは、小学生の算数レベルの基礎ができていないからです。

続く

タグ：数楽統計

posted at 07:34:35

#統計 #数楽

同様に「単純線形回帰の尤度函数」の話を聞いて、義務教育レベルの数学でやったことの一般化に過ぎないことを全然認識できない人達は、中学校レベルの2次式の平方完成を十分に理解できていないことになると思う。

この手のことはどんなに高級な数学を知っていても起こり得る。

タグ：数楽統計

posted at 07:37:36

#統計 #数楽使っている用語だけを見ると高級そうな話題に見えていても、実際には、小中学校レベルの算数や数学の基本的な考え方に還元できてしまうことは少ないと思う。

数学では基礎が重要だと言われるゆえん！

タグ：数楽統計

posted at 07:40:48

以下のリンク先の話題は

yu-cho-f.jp/wp-content/upl...…
偏差値45未満の【当該私立大学文系学部卒者の期待生涯賃金の平均値は高卒平均値(1億9,203万円)を上回っているということである。その額もおおよそ3,700万円程度となっており、大学進学に伴う経済的効果が確認される】

とも関連付けられるべき。 twitter.com/daigaku_diary2...

タグ：

posted at 07:46:44

人道的な観点からだけではなく、偏差値45未満の大学であっても進学が経済的に合理的な場合があるという事実からも、「大学に行く人を減らすべきである」という類の言説は長期的には日本をダメにする可能性が高い意見だと考えられます。

twitter.com/genkuroki/stat...

タグ：

posted at 07:51:03

以下のスレッドも参照

twitter.com/ropcb08/status...

タグ：

posted at 07:51:49

Y Tambe @y_tambe

識者の反応「ま　　た　　モ　　ン　　タ　　ニ　　エ　　か」

タグ：

posted at 07:54:22

やばい、さっきいれたのでコーヒー豆が尽きた。

大好きな自作コーヒー牛乳がこれ以上飲めなくなるではないか。

タグ：

posted at 07:59:54

Nyoho @NeXTSTEP2OSX

Teamsが

• やっと9人まで同時表示できるように(今まで最大4)
• 挙手ができるように
• 背景が付けられるように(今まではぼかすだけ)

なった。
< What’s New in Microsoft Teams | April 2020 - Microsoft Tech Community - 1347929 techcommunity.microsoft.com/t5/microsoft-t...

タグ：

posted at 08:37:01

Nyoho @NeXTSTEP2OSX

ちなみにわしのではぼかす機能すら出てきません。

タグ：

posted at 09:15:09

Nyoho @NeXTSTEP2OSX

TeamsはAVX2.0 (Advanced Vector Extension 2) を使っているからというのが理由だと先日伺いました。しかしZoomはバーチャル背景もできているのでZoomはAVX2.0を使わずにやっていることがわかります。

タグ：

posted at 09:17:58

#Julia言語 LibSndFile.jlの使用例の追加

nbviewer.jupyter.org/gist/genkuroki...

タグ： Julia言語

posted at 12:20:45

nao @parasite2006

@samakikaku @pinpon_2011 @konamih @syoyuri つぶやきを使わせていただきました。問題ありましたら対応しますのでご連絡ください。togetter.com/li/695815?page=2

タグ：

posted at 13:40:08

yuri @syoyuri

まとめありがとうございます。
放射能が消えるかのよううたい批判された「創生水」に関係した樺山繁氏が、今は日本トリムにいて、新型コロナ“研究”に動きをみせていること、消費者庁も厚生労働省も注視して欲しいですね。 twitter.com/parasite2006/s...

タグ：

posted at 13:43:54

Keisuke @keisukefukuda

「1から始める Juliaプログラミング」届いたので連休中に読む。「１から始める」ってのが気が効いている（Juliaは配列のインデックスが1-origin）

タグ：

posted at 13:44:45

tsuchita @thnms_tct

横からみるとこんな感じです。 pic.twitter.com/G1dm6BnnxM

タグ：

posted at 14:14:50

完全理解者爆誕　#Julia言語 twitter.com/karaage0703/st...

タグ： Julia言語

posted at 14:45:38

１０万円給付申請受け付け　宮城の全市町村が月内開始 | 河北新報オンラインニュース / ONLINE NEWS sp.kahoku.co.jp/special/spe121... pic.twitter.com/e9bcRSo28l

タグ：

posted at 16:44:31

#超算数まず長方形をわかりやすい言葉で定義し、次に正方形を長方形の一種であることが分かる定義をするスタイルは、教科書会社による工夫かもしれません。当時有効だった指導要領(には詳しいことが書いていないのでその)解説では、正方形が先に定義されています(www.mext.go.jp/component/a_me...: 93頁)。 twitter.com/genkuroki/stat... pic.twitter.com/b4Opg7CzKZ

タグ：超算数

posted at 17:19:08

#超算数解説算数の定義はどちらかといえば、正方形は長方形と別のものであるという主張として読まれることが多いのかもしれません。また、算数教育界の見解では、1998年の指導要領で【正方形と長方形とは別の形として指導するように決められ[た]】twitter.com/OokuboTact/sta...そう。

タグ：超算数

posted at 17:19:10

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

#超算数 98年指導要領には詳しいことはなんにも書かれていませんwww.nier.go.jp/guideline/h10e...。算数教育指導用語辞典は業界のお約束どおり、指導要領を誤って定義しています。そこで指導要領解説算数編(2007 [1999]: 100)を参照すると、上掲の2009年版とほぼ同じ文言がありました。 pic.twitter.com/h7yAzOOEOC

タグ：超算数

posted at 17:19:11

Stay home, write Julia.

タグ：

posted at 17:25:10

#統計 #数楽使っている用語だけを見ると高級そうな話題に見えていても、実際には、小中学校レベルの算数や数学の基本的な考え方に還元できてしまうことは多いと思う。

数学では基礎が重要だと言われるゆえん！

【以上は訂正版。「少ない」→「多い」(笑)】

タグ：数楽統計

posted at 17:38:38

#超算数本当は、更に一つ前の指導要領(にはなんにも書いてないのでそ)の解説本である
文部省『小学校指導書算数編』東京、東洋館出版、1989年。MESC 1-8906
を見て、そこで正方形と長方形がどう定義されているかを確認したいところですが、暫く実行できません。間接的な検証法を取ります。

タグ：超算数

posted at 17:50:30

#超算数図形の包摂関係を小学校で教えることを唱道したのは、中島健三です。その『算数・数学教育と数学的な考え方』(1981)によれば、77年指導要領www.nier.go.jp/guideline/s52e...にある6年の【基本的な図形の相互関係について調べること】は、【趣旨が多少あいまいになっている】(16頁)が、包摂関係の指導。

タグ：超算数

posted at 17:50:31

Atsushi Sakai @Atsushi_twi

Bye-bye Python. Hello Julia! towardsdatascience.com/bye-bye-python...

タグ：

posted at 18:18:46

Matt Henderson @matthen2

certain Julia sets can be 'mated' - glued together so they cover the whole complex plane (shown here on a sphere). The resulting picture is a characterization of the Julia set of (z²+ c₁) / (z²+ c₂) pic.twitter.com/456wLLcPlP

タグ：

posted at 18:59:55

Matt Henderson @matthen2

made by adapting code from Claude Heiland-Allen mathr.co.uk/blog/2020-01-1...

タグ：

posted at 19:00:28

Matt Henderson @matthen2

I was inspired by Sarah Koch's talk www.youtube.com/watch?v=YoibJg...

タグ：

posted at 19:01:15

#統計さらに訂正

twitter.com/genkuroki/stat...

❌ p(y|x,a,b,σ) = (1/√(2πσ²))exp(-(y-(a+bx))²/σ²)

⭕️ p(y|x,a,b,σ) = (1/√(2πσ²))exp(-(y-(a+bx))²/(2σ²))

exp内の分母の2を忘れていた。実際に使うのは

-2log p(y|x,a,b,σ) = (y-(a+bx))²/σ² + lod(2πσ²)

の方。確率はその対数を扱うと便利。

タグ：統計

posted at 19:37:11

#統計

問題: Anscome's quartetの4つのデータセットがほぼ同一の単純線形回帰の与える理由は何か？単純線形回帰のベイズ版でも通用する説明を述べよ。

という問題は、線形回帰について理解しているかどうかを確認するためには非常に良い演習問題になると思いました。

タグ：統計

posted at 19:42:12

#統計データセット (y_i, x_{i1}, …, x_{im}) (i=1,…,n) の線形回帰

y_i = x_{i1}β_1 + … + x_{im}β_m + ε_i, ε_i ～ Normal(0,σ)

の対数尤度函数の-2倍は、適当に縦ベクトルy, βと行列Xを定めると

L(β) = ||y - Xβ||²/σ² + n log(2πσ²)

と書けます。||・||はユークリッドノルム。続く

タグ：統計

posted at 19:52:06

#統計続き～これは、転置を ' で表すと

L(β) = β'X'Xβ - 2(X'y)'β + y'y

と書き直されます。だから、モデルの有度函数は X'X, X'y, y'y で決まる。それらが等しければ、同じ尤度函数が得られます。

尤度函数が同じなら、最尤法(この場合には最小二乗法と一致)やベイズ推測の結果も同じになる。続く

タグ：統計

posted at 19:52:07

#統計 Anscome's quartetの4つのデータセットは実際にほぼそうなるように作ってあるわけです。ただし、データセット (x_i, y_i) から x_{i1}, x_{i2} を x_{i1}=1, x_{i2}=x_i で作る(m=2)。

ここまで解説を読めば

nbviewer.jupyter.org/gist/genkuroki...

のスクショである添付画像1の内容を理解できると思う。 pic.twitter.com/yU8csVqDrM

タグ：統計

posted at 19:56:11

#統計 #数楽上の方では「平方完成」を強調する説明をしてしまいましたが、Anscome's quartetが単純線形回帰で同じ推定結果を与えることの理解には「平方完成」さえいらず、本質的に

(ax-b)² = a²x² - 2abx + b²

の一般化に過ぎない

||Xβ-y||² = β'X'Xβ - 2(X'y)'β + y'y

だけで十分。訂正に続く

タグ：数楽統計

posted at 20:02:17

#統計訂正

❌ L(β) = β'X'Xβ - 2(X'y)'β + y'y

⭕️ ||y - Xβ||² = β'X'Xβ - 2(X'y)'β + y'y

左辺を間違った。 twitter.com/genkuroki/stat...

タグ：統計

posted at 20:04:28

#数楽多分、

(ax-b)² = a²x² - 2abx + b²

の計算なら「易しい」と思える人は結構いると思う。行列やベクトルの計算にも同様に慣れることができれば、同じ感覚で

||Xβ-y||² = β'X'Xβ - 2(X'y)'β + y'y

も「易しい」と感じられるようになります。

その瞬間にAnscome's quartetも易しくなる。

タグ：数楽

posted at 20:08:26

#数楽正しく学べば、個々の具体的な状況ごとに「やり方を暗記」するのではなく、基本的な計算や推論に習熟していることを使って、複雑で高級に見える事柄であっても、義務教育レベルの算数や数学と同じ感覚で扱える基本的な事柄に帰着してしまうことができるようになります。

これは実際に便利です。

タグ：数楽

posted at 20:12:18

optical_frog @optical_frog

スワン先生の解説が少し不十分だったので，the の発音について補足してたときに，動画の音声がでなかったので，こちらでフォローしてみますかね．まず，大きく分けて，「強い発音」と「弱い発音」があります．そんで，中学校英語で習う「母音の前 vs. 子音の前」の区別は弱い発音のなかの区別です． pic.twitter.com/JCkrHV92Sw

タグ：

posted at 20:22:58

マナマの生マナマコ @chor0sK

直線上への射影 #julialang #Emacs pic.twitter.com/4g6AbuWKaf

タグ： Emacs julialang

posted at 20:31:42

マナマの生マナマコ @chor0sK

Emacs Org-ModeはJupyter Notebookよりスクショの見た目が良い(過言の滝)

タグ：

posted at 20:34:54

TaKu @takusansu

#超算数 @sekibunnteisuu 「「試行錯誤」をさせようとしない人達がいるようだ」という題名で掲示板に書き込みました。
8254.teacup.com/kakezannojunjo...

タグ：超算数

posted at 20:51:16

Akinori Ito @akinori_ito

Google Meetに参加者を一人ずつ呼び出すためのWebアプリを書きました。α版ですが暇な人は遊んでみてください
aitoweb.world.coocan.jp/WaitingRoom/re...

タグ：

posted at 21:06:20

マナマの生マナマコ @chor0sK

combinedims([v1,v2,v3])よりもhcat(v1,v2,v3)のほうがタイプ数少ない件について #julialang

タグ： julialang

posted at 21:13:34

Marco Godano @GodanoMarco

Let me briefly explain my thought process when I have to translate a line from Japanese into English.
The line in question is "・・・うひいいいいいいいい～！！".
(thread: 1/47)

タグ：

posted at 21:21:29

岸政彦 @sociologbook

稲葉振一郎 @shinichiroinaba が踊ります

aidn.jp/wowa/26573512

タグ：

posted at 21:24:44

Prof. Feynman @ProfFeynman

OokuboTact　大久保中二病中年 @OokuboTact

Your degrees aren't as valuable as your skills. 🧠

タグ：

posted at 21:40:16

OokuboTact　大久保中二病中年 @OokuboTact

@temmusu_n 参考　中島健三の図形の包摂関係 pic.twitter.com/moLdNQ7Z0R

タグ：

posted at 21:56:27

@temmusu_n （続き）

中島健三の図形 pic.twitter.com/dL5X94ukuF

タグ：

posted at 22:00:27

Interface ６月号の #Julia言語特集は５日間で書いたので読む側もGWで一通り読めるようになるはず（ほんまか？）

続く：

タグ： Julia言語

posted at 22:16:00

例えば

- REPLを触る
- REPL の上でコードをコピペする
- Jupyter Notebook でペチペチ遊ぶ
- ローカル環境でJuliaのProjectを作る
- 適当なアプリを作る

で5日で行けるかな（白目）

続く：

タグ：

posted at 22:16:01

堅苦しい文法というセクションは完全に排除してサンプルコードで必要そうな文法が登場するようにした．

素材例は素朴に，素材を固定しそれを活用する使う道具を複数紹介し「自分でコードを書くと行った応用場面であの本で触った概念をそのまま使えそう」と誘導ができるように書いたつもりではある．

タグ：

posted at 22:16:01

ゴマ「ダメだったら飼い主が悪かったということで」
飼い主「(´・ω・｀)ウゥゥ」

タグ：

posted at 22:16:01

怠惰に暮らしたいトド @tsukuru_ouu

ちなみにJuliaのアプリケーションを作る章での最終成果物は

github.com/terasakisatosh...

で公開しています．完成に至るまでのプロセスは本で書いている通り．

タグ：

posted at 22:16:01

OokuboTact　大久保中二病中年 @OokuboTact

思考力に優れている方でも、自身の専門分野の「作法」にはとらわれてしまいます。

理論物理学の押川教授は、不十分な情報で方針決定するのを仕事とする医者の判断を、後知恵で「外しまくって」と呼びます。

一方、医者は与えられたリソースで最善を尽くすので、リソース自体の拡充の考えは弱いです。 twitter.com/masakioshikawa...

タグ：

posted at 22:41:29

@temmusu_n @takusansu

（続き）　目次 pic.twitter.com/Q6XhpmTjak

タグ：

posted at 22:43:45

押川さんが特別に変なだけで、専門分野の「作法」にとらわれるという問題ではないと思います。

「分野A vs. 分野B」型の話題にすることはこのタイミングでは我々の社会にとって非常に有害。

押川さん個人を「この件では他の理論物理学者よりも劣っている」と批判する方が良いと思いました。 twitter.com/tsukuru_ouu/st...

タグ：

posted at 22:55:14

@OokuboTact @takusansu #超算数ということは、ここtwitter.com/ookubotact/sta...で、包摂関係への深入りをちょっと批判がましく云々しているのは中島健三ということなんですね。

タグ：超算数

posted at 22:55:53