黒木玄 Gen Kuroki(@genkuroki)/2020年12月17日

厳密に言えば、C言語で書けない方法でC言語を超えることが出来れば、より速い言語と言えるけどC言語で書けないことがないので表現としては誤りかな🧑🏻‍💻
ただ速いコードに行き着くまでのスピードはJuliaの方が速いと思う。 twitter.com/mattn_jp/statu...

タグ：

posted at 00:18:10

非公開

タグ：

posted at xx:xx:xx

時系列や空間系列などのi.i.d.仮定をおけない場合についての先生のコメントはここで読んで、この辺りは研究になる話なんだなという認識。

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

タグ：

posted at 01:06:49

エコ鉄 @eco_tetsu

私は，別の方法（別解）を見つけるのを楽しんでいた．

教えた通りにしろと言う人．
それもありだという人．

うちの親は後者なので平和に過ごせました． twitter.com/genkuroki/stat...

タグ：

posted at 01:10:57

K. Yam @yam_cpp

ちょっとpython bindingをいじるたびにC++のビルドを数分待つのがなかなかつらい．．Julia使いたい．．

タグ：

posted at 01:18:48

とはいえMTBS2018にそれぞれの時点での実現値が得られたもとでの条件付独立での適用の例は挙げられています。

タグ：

posted at 01:23:37

Danner @danner_2009

ほんこれ。
小3,4年の担任が創意工夫、自発性を褒めてくれる先生で授業が楽しかった。この時期に運悪くハズレのパターンマッチング教師にぶつかっていたら、勉強が嫌いになっていたと思う。 twitter.com/genkuroki/stat...

タグ：

posted at 01:25:20

統計たん @stattan

カッコ良すぎる。。。「神経回路網を理解するために、ただそのためだけに、具体的な計算のくりかえしの果てに、それらの方法こそが、この問題の自然であるということに到達したのです。もちろん、自然な理解に到達したあと、その基礎の上でさらに新しい具体的計算をくりかえしています。」 twitter.com/katzkagaya/sta...

タグ：

posted at 01:27:40

それを「予測誤差」というのは上のページでも言われているように、僕も腑に落ちないし、別時点でのデータ点を評価に使うのが一番無難に思える。

タグ：

posted at 01:28:17

統計たん @stattan

完全に感謝の正拳突きの境地。やばい。熱いぜ。。。 twitter.com/stattan/status...

タグ：

posted at 01:28:30

非公開

タグ：

posted at xx:xx:xx

#統計

答えはこうです！

帰無仮説「表の出る確率はθ以上である」について、「コインをn回投げた場合」と「ちょうどk回表が出るまでコインを投げた場合」の片側検定のP値は、それぞれ事前分布がBeta(1,0), Beta(0,0)の場合の事後分布で帰無仮説が成立する確率に一致する。

nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat... pic.twitter.com/vzptQGnvxw

タグ：統計

posted at 05:55:29

#統計つまり、「止め方の違い」は「事前分布の違い」にちょうど対応している！

ただし、ベータ分布Beta(a,b)はa,bがともに正の場合にのみ定義されているので、Beta(1,0), Beta(0,0)は所謂improper priorsです。a,b>0として事後分布を計算した後に、(a,b)→(1,0), (0,0)と極限を取ると考えてもよい。

タグ：統計

posted at 06:01:07

#統計要するに、コイン投げのデータの取得法の違い(実験の止め方の違い)を知っている人が、事前分布に実験の止め方の情報を組み込めば、通常の仮説検定におけるexactなP値と完全に同じ数値(その数値は止め方によって違う)を事後分布において帰無仮説が成立する確率として得られるわけです。

タグ：統計

posted at 06:06:07

#統計添付画像のOut[7], Out[8]に0が並んでいるのは、2種類の片側検定におけるexactなP値が、異なるimproper事前分布に関する事後分布内で帰無仮説が成立する確率にぴったり一致していることの確認になっています。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/uNFWeIW9TT

タグ：統計

posted at 06:08:57

#統計帰無仮説「表の出る確率はθ以上である」の片側検定について

(1)n回投げたら表がk回の場合のP値
(2)ちょうどk回表が出るまでn回の場合のP値
(3)フラット事前分布の事後分布で帰無仮説が成立する確率

を比較。上で述べたように(1),(2)も事後分布で解釈可能。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/z0Dj9OeuKA

タグ：統計

posted at 06:20:06

#統計「コインを投げる回数nを固定した場合」と「固定したkについてk回表が出るまでコインを投げ続けた場合」で、尤度函数は定数倍の違いしかないので、同じ事前分布のベイズ統計の結果は同じになってしまうが、事前分布に実験の止め方の情報を組み込めば2つの場合に違う結果を得ることができる。

タグ：統計

posted at 06:25:15

#統計 20世紀の主観主義のベイズ主義者の有力者達は、尤度函数に定数倍の違いしか無ければ結論は同じにならなければいけないという説(所謂「尤度原理」というデタラメ)を信じてしまったようだ。

しかし、素直にベイズ法を眺めれば、事前分布に「止め方の違い」の情報を適切に組み込むことができる。

タグ：統計

posted at 06:29:59

#統計ベイズ統計ではモデルの設定に「事前分布」の自由度が増えており、その増えた自由度を使って、それ抜きだと不可能だったことをできるようになっている。

私がベイズ統計の枠組みが優れていると思うのはそういう理由である。

「主義」によって事前分布の使用法を制限するのは極めて有害。

タグ：統計

posted at 06:36:57

#統計ベルヌイ分布モデルにおいて「止め方の違い」は等しい事前分布によるベイズ統計に結果に現れない。そのままだと、P値を計算する統計分析にベイズ統計は劣っていることになってしまう。

しかし、このスレッドで示したように、事前分布の取り方を適切に変えることによって、異なるP値を再現可能。

タグ：統計

posted at 06:36:58

#統計ところが、「主義」に基くベイズ主義のベイズ統計の信者達は、仮説検定ではできていたデータ取得法の区別をベイズ統計ではできないと勘違いし、「尤度原理」を満たすベイズ統計でできないことを仮説検定ができてしまうのは仮説検定が間違った統計学だからだと主張した。

単純に馬鹿だと思う。

タグ：統計

posted at 06:41:07

#統計補足：単純なベルヌイ分布モデルでベイズ統計を使う必要は全然ないと思います。「平方根を計算する人」になってどんぶり勘定すれば良い。一般に総和のばらつきの大きさは√nに比例するので、何はともあれ√nを概算したくなる。

目的ごとに適切な手抜きをすることが、数学の応用ではよい。

タグ：統計

posted at 06:49:55

名目換算とか購買力平価換算で国際比較するのはやめた方が良いよね。

タグ：

posted at 06:53:52

#統計以下のリンク先に書いたことについては、「何を言っているのかわからない」と感じる人が多いと思う。

なぜならば「コインを投げる回数nを固定」「表がちょうどk回出るまでコインを投げた場合」で統計分析の結論が同じにならなければいけない、と本気で言う人の存在を普通は想定できないからだ。 twitter.com/genkuroki/stat...

タグ：統計

posted at 07:03:09

#統計これは、掛算順序問題に似ている。

「2×8ならタコが2本足になる」とか「x円のノート8冊の文脈でも、8×xは8円のノートx冊の意味になる」と本気で言っている人がいることを普通は想定できない。

だから、掛算順序問題について保護者達に注意を喚起しても以前はなかなか信じてもらえなかった。

タグ：統計

posted at 07:07:35

Clément Canonne @ccanonne_

Stuff I wish I had known sooner: "Pinsker's inequality is cancelled," a thread. 🧵

If you want to relate total variation (TV) and Kullback-Leibler divergence (KL), then everyone, from textbooks to Google, will point you to Pinsker's inequality: TV(p,q) ≤ √(½ KL(p||q) )

1/ pic.twitter.com/NXm03yUg2f

タグ：

posted at 07:15:31

#統計続き↓ twitter.com/genkuroki/stat...

タグ：統計

posted at 07:23:45

Han Zhao @hanzhao_ml

Also, there are related discussions in Section 5 of Prof. Yihong Wu's lecture notes: www.stat.yale.edu/~yw562/teachin... twitter.com/ccanonne_/stat...

タグ：

posted at 07:36:25

#統計事例：所謂豊田『瀕死本』では、「尤度原理」に基いて、データの取得法が違っていても気にする必要がないと強調している。

確かに(定数倍の違いを無視して)尤度函数の情報だけを使う統計学がまともなら、気にする必要はない。

もちろん、そうではないのだが。

研究不正し放題の世界(笑) pic.twitter.com/SStIZEaZbL

タグ：統計

posted at 07:56:30

#統計補足↓ twitter.com/genkuroki/stat...

タグ：統計

posted at 07:59:08

非公開

タグ：

posted at xx:xx:xx

るこら。 @rucola_

「おかあさん、どうしてこれいちいち書かなきゃいけないの。かかなくても答えがわかるのに」
と言って息子に泣かれた小2のときに色々思い切った対策取るべきだったと後悔しています。 twitter.com/genkuroki/stat...

タグ：

posted at 08:34:19

#数楽線形代数ネタ

What Is the Matrix Sign Function? – Nick Higham
nhigham.com/2020/12/15/wha...

タグ：数楽

posted at 08:35:00

ランダム分割した領域でモデルのテストができて的中率が云々という論文が批判されるということ近年ありましたね…そのとき言及はしませんでしたが…

タグ：

posted at 10:08:03

Inazuma110 @Inazuma_110

講義でJuliaでPythonにマウンティングするためみたいなコードが要求された pic.twitter.com/VbnyjjqRBN

タグ：

posted at 10:08:11

大須賀覚 @SatoruO

峰先生@minesohの本。超絶おすすめ

誰もが抱くコロナの疑問に徹底的に答えてます。網羅的で深い内容なのに、対話形式ですごくわかりやすい。深い知識に裏打ちされているからこそできる明確な回答が気持ちいい

コロナ本の決定版。これ読んでおけば、ひとまず安心だと思いますamzn.to/2KC5Zw8 pic.twitter.com/iMPec4bq7d

タグ：

posted at 11:11:50

仙台弁こけし @jugokeshi

宮城県さ大雪警報でたっちゃ〜 pic.twitter.com/RmzIhKX84R

タグ：

posted at 13:10:40

音を置き去りにした twitter.com/stattan/status...

タグ：

posted at 13:20:33

yudai.jl @physics303

Julia言語内でPyhonを呼びたかったら，using PyCall を使えば良い？

タグ：

posted at 14:17:41

Yuki Nagai @cometscome_phys

@physics303 この記事の後半参考になるかもしれません。

qiita.com/cometscome_phy...

タグ：

posted at 14:23:31

yudai.jl @physics303

@cometscome_phys ありがとうございます！top levelだとpy""" """ は不要なんですかね．

タグ：

posted at 14:25:44

Yuki Nagai @cometscome_phys

@physics303 一番下がJuliaのコード、下から2番目が等価なPythonのコードですので、見比べれば使い方がわかると思います

タグ：

posted at 14:30:59

河合祐介 @tkawai18_tkawai

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

こういうのを「工夫」だと思っている人がいて寒さで固まっているｗ pic.twitter.com/hCajZcF6ZB

タグ：

posted at 14:35:18

非公開

タグ：

posted at xx:xx:xx

放射能デマが甲状腺検査を食い物にしようとしていますが、学校での甲状腺検査がもつ半強制性こそが今いちばん批判されているところです。
甲状腺検査そのものを中止するのが望ましいと僕は考えます。しかし、とにかく学校検診は即刻中止すべきです。
放射能デマの言うことを聞いてはだめです twitter.com/makomelo/statu...

タグ：

posted at 15:13:39

非公開

タグ：

posted at xx:xx:xx

Tom Kwong @tomkwong

川端裕人『ドードー鳥と孤独鳥』（国書刊 @Rsider

#AdventOfCode Day 17 animation! This is also my first 3D plot 🤟
Powered by: #julialang pic.twitter.com/pwactrvdr8

タグ： AdventOfCode julialang

posted at 16:10:32

川端裕人『ドードー鳥と孤独鳥』（国書刊 @Rsider

【2x2表で考えるヒント】
眼科医の人が、色覚検査をしないかったために、就職の時にひどいめにあうなどの被害が出ており、色覚検査をすべきというふうに言う時、問わなければならないのは、色覚検査が「ひどいめ」を抑止する効果があるかどうか。そのためには下の表を埋めて考える必要がある。 pic.twitter.com/vC8kRsIkG5

タグ：

posted at 16:33:24

川端裕人『ドードー鳥と孤独鳥』（国書刊 @Rsider

でも、多くの場合、先の図のaにあたる人がいたということが示されるだけ。それだけでは、色覚検査と関連しているのかわからない。実際、色覚検査をしていたけどひどい目にあったcの人も結構たくさんいる。色覚検査をしなかったけどひどい目にあわなかったbの人はもっとたくさんいる。

タグ：

posted at 16:34:56

結局、ad/bc（オッズ比）を計算してみないといけないよと某所で説明するはめになった。聞く耳をもっていただけただけ、マシであった。

タグ：

posted at 16:35:23

テキーラ氏 @kiu_x

@yusai00 今年１月、僕の車のガラスルーフに出たやつです pic.twitter.com/wMYYa9v6Kz

タグ：

posted at 17:00:52

takeruiso @owl8032

python/C APIを使ったモジュール作りかけてるけど、juliaに移行しようかなという気になってきた。

タグ：

posted at 17:10:45

ひびき @hibikiw

「理論疫学者・西浦博の挑戦」を読み終えたのだけど、三歩先を照らす数理モデル、伝播の動きを大きく変える行動変容、一定以下の感染数に押さえ込まれた時に威力を発揮するクラスター対策など優れた道具は既に手に入れているのに、その使いどころは社会のコンセンサス次第なのだと改めて感じた。

タグ：

posted at 17:14:27

Yuki Igarashi @bonprosoft

新青葉山キャンパス、だいぶ雰囲気が変わった気がする
ライトがもう少し明るくなれば青葉山ゲレンデでナイターいけそう pic.twitter.com/qFecZta8hb

タグ：

posted at 17:18:46

ひびき @hibikiw

最初の流行時はコンセンサスを醸成している時間がなかったので学者自ら語りかけたけど、これはやはりおかしいよね？と然るべき相手にバトンタッチしたら、逆にコンセンサスどころか半歩先の未来しか「見ない」状況に陥ってしまったのが勿体ない。でもこれはバトンを受けた側、そして我々の責任だよね。

タグ：

posted at 17:27:40

がっくん @gifted_te_nani

頑なに書くことを拒否。 pic.twitter.com/lnfYKwAZMJ

タグ：

posted at 17:32:11

M1-Macに乗り換え @kabukauna

M1 MacのRosetta2でインストールしたJulia1.53で、
RCall.jlがR_HOME is not a directory.と言われてビルド失敗した

Apple siliconが原因かもと心配だったけど、
ENV["R_HOME"]="*"
でビルド成功
JuliaのConda領域にインストールされた
#Julia言語 #M1Mac pic.twitter.com/6wpB9V7fmZ

タグ： Julia言語 M1Mac

posted at 18:05:09

akiaki5516 @doradora09

とりあえずJuliaは欲張りな言語というのは伝わるw pic.twitter.com/wGAdyezCxe

タグ：

posted at 19:16:44

猫のチャコとチー坊 @chakochiibou

#統計ない情報を勝手に読み取るのはまずい。

例えば、このスレッドでは「表が出る確率はθ以上である」という仮説を扱う場合にしか計算例を示していない。

パラメータの推定の場合は別の話になる。

掛算順序問題でも「教科書の教師用指導書」を「教科書ガイド」と読み取った人がいた(笑)

タグ：統計

posted at 20:10:41

床暖地蔵しながらエクササイズの応援してくれてます🐾
＃ネコ　＃ステイホーム pic.twitter.com/g6QP6pUjjj

タグ：ステイホームネコ

posted at 20:11:16

TaKu @takusansu

twitter.com/cafe_trico/sta...
教師等に逆らわない、奴隷精神の人材をピックアップしたいのかな。

タグ：

posted at 20:24:40

@heart8255

「検査を増す。陰性だとわかれば安心して旅行できる。」

これが間違いのもとだと誰か教えてあげて下さい(-人-;)

. twitter.com/CDP2017/status...

タグ：

posted at 20:50:16

TaKu @takusansu

@temmusu_n 参考文献情報
海外における割合指導に関する先行研究の分析
熊倉啓之
www.jstage.jst.go.jp/article/jssep/...
「くもわ」的なツールは割合の理解を深める事には繋がらないとばっさり。
3用法の指導は、他の方法に比べて最も効果が小さい。
等、書かれています。

タグ：

posted at 21:02:13

k @musicisthebest_

新たな症例を観測しました。 pic.twitter.com/54ja2KrS8m

タグ：

posted at 21:12:38

Hiromitsu Takagi @HiromitsuTakagi

有線電気通信法を知らない者はこのように多い。
twitter.com/tanrakukairo/s...

タグ：

posted at 21:36:25

Akinori Ito @akinori_ito

日本では教育のIT化っていうとタブレットなのはなぜだろう。インフラが整ってない国でさえOLPCなのに

タグ：

posted at 22:20:53

@temmusu_n @yamazaksv2 元ツイートが消されているので、貼り付けておきます。 pic.twitter.com/OBrQ01139t

タグ：

posted at 22:20:55

@sekibunnteisuu @yamazaksv2 #超算数ありがとうございます。ちなみに元ツイはアーカイブ済みです。

タグ：超算数

posted at 22:28:28

@temmusu_n @yamazaksv2 了解しました。

タグ：

posted at 22:29:35

@takusansu #超算数ありがとうございます。
Teahan, Teresa. "How I Learned to Do Percents." _The Arithmetic Teacher._ Vol. 26, No. 5 (January 1979), pp. 16-17 www.jstor.org/stable/41187744
ティーハンさんがobnoxious triangleの提唱者だということですね。Teahan's triangleで検索すると数件がヒット。

タグ：超算数

posted at 22:29:46

S.mit.Crystal @mitstream

学生時代には数値計算と言えばFortran一択だった身ですが、最近そのスキルがjuliaで使えそうだと知って、俄然juliaに興味がわいています。

タグ：

posted at 22:31:21

@sekibunnteisuu @yamazaksv2 アカウントは、鍵かけちゃいましたね。

タグ：

posted at 22:32:23

@temmusu_n @yamazaksv2 @nao03134522　ほんとだ。

　しかし、ツイッターで遭遇する小学校教師アカ、首をかしげたくなる人が多い。

タグ：

posted at 22:35:53

@takusansu #超算数ジャストルはコロナ対応で毎月100本の論文が無料で読めるようになっています。ティーハンさん、他にもこういう図が開発できるのではと述べています。オームの法則とか、類例は多いですよね。彼女の独創ではないですが、National Council of Teachers of Mathematicsの真面目な雑誌に載った。

タグ：超算数

posted at 22:40:33

#超算数
高1が学校の採点に憤慨していた。

有理式の計算問題　分子をx+1とすべきところとを(x+1)と括弧を付けたがために誤答扱い。

配点3点で0点だそうだ。

「無視していい」と言っておいたけど、

テストで何を測りたいのか、何を理解させたいのか、甚だ疑問

タグ：超算数

posted at 22:54:08

バク @baku_reshi

RustとかJuliaとかも学びたい

タグ：

posted at 22:56:05

@takusansu #超算数 Teahan's triangleと掛算の順序だけが違う図がありました。
Boling, Bruce A. "A Different Method for Solving Percentage Problems." _The Mathematics Teacher._ Vol. 78, No. 7 (October 1985), pp. 523-524. pic.twitter.com/FJvmMy4YFK

タグ：超算数

posted at 23:18:57

@takusansu #超算数少し毛色が違うのですが問題文のisとofを図式に流し込む覚え方(mnemonics)として、
McGivney, R. J., John Nitschke. "Is — of, a mnemonic for percentage problems." _Mathematics Teacher._, Vol. 81, No. 6 (September 1988), pp. 455-456. www.jstor.org/stable/27965881 pic.twitter.com/d45UhXwG4m

タグ：超算数

posted at 23:22:27

OokuboTact　大久保中二病中年 @OokuboTact

@takusansu #超算数マギヴニさんとニチキさんは、isにあたる数が上でofにあたる数が下と。理由はアルファベット順という暗記法。
Students find the formula easy to remember because "i" comes before "o" in the alphabet, and we write the numerator of a fraction before the denominator. (455)

タグ：超算数

posted at 23:24:35

#超算数

謎の交換法則 pic.twitter.com/PJ0I0nUtow

タグ：超算数

posted at 23:26:44

非公開

タグ：

posted at xx:xx:xx

@takusansu #超算数彼らのis-ofは掲載誌の通信欄で批判されました(1989年、サウス・オレゴン州立大学のJohn Whitesittさん)。生徒に考えさせる目的のこの雑誌で頭を使わなくてよい求答法とはどういうことかという趣旨。同じ紙面で、彼らが反論www.jstor.org/stable/27966257。1、普通の解き方と平行して指導する。

タグ：超算数

posted at 23:31:53

峰　宗太郎 @minesoh

きゃっきゃ、ちょっとアレなタイトル👶💦 twitter.com/TOTB1984/statu...

タグ：

posted at 23:40:44

@takusansu #超算数 2、数学で暗記法を使うのは伝統的なやり方。ここでπ<θ<3π/2でtanθ>0を教えつつ語呂も使うという例を出した。志の低さに頭を抱えてしまう。3、is-ofでできない子にも成功の喜びを味わわせることができる。
批判も擁護も日本と同じ。きはじが確認された初期からですtwitter.com/temmusu_n/stat...。

タグ：超算数

posted at 23:46:48

少しでもいい影響が広がりますようにと願いを込めて描きました。もし、読んでよかったよと思ってくださるかたがいたら、リツイートしていただけると喜びます。
#統合失調症 pic.twitter.com/79lhahr7Pv

タグ：統合失調症

posted at 23:48:08

pic.twitter.com/9kreslD23w

タグ：

posted at 23:48:10

pic.twitter.com/Pt0J3WUKAI

タグ：

posted at 23:48:12

pic.twitter.com/g7zTMUyA4m

タグ：

posted at 23:48:13

pic.twitter.com/cmkdAIDi7w

タグ：

posted at 23:48:15

pic.twitter.com/ilP1glGj1D

タグ：

posted at 23:48:17

pic.twitter.com/hQMk57lU4n

タグ：

posted at 23:50:45

pic.twitter.com/yJxHWKXKDg

タグ：

posted at 23:50:47

おまけ pic.twitter.com/UjF0u5LPVN

タグ：

posted at 23:50:48

#Julia言語おお！RCall.jlのビルドがそういう仕様になっていることを知らなかった。

Juliaを使い始めて、PythonやRなどの便利さを再確認するようになりました。

JuliaからRが呼べるようになっていると、統計関係のコードを書くときの確認をやり易くなる。 twitter.com/kabukauna/stat...

タグ： Julia言語

posted at 23:50:53

絵本を出版する前に公開した作品を、また投稿しました。出版したものと文章が若干違う、絵は全て違うのでご容赦ください。(絵本は、長文が読みにくい方を考え、この作品の文章をまとめています。。)願いを込めた作品です。

タグ：

posted at 23:56:59

ときどき、「Juliaを使うなら全部Juliaでやらないと意味がない。PythonやRを呼ぶならJuliaを使う意味がない」のような発言を見かけるけど、ものすごく間違っていると思う。

PythonやRからもJuliaをどんどん使えば良いと思う。

タグ：

posted at 23:57:05