黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2019年11月04日(月)
ごまふあざらし(GomahuAzaras @MathSorcerer
PILでよくあるのが () いるんだっけ?問題があるけれど
#Julia言語 だと multi dispatch でよしなにしてくれる.
```
function goma(xs)
return sum(x for x in xs)
end
function goma(x::T, y::T...) where {T<:Integer}
return goma((x, y...)) # 上のゴマを呼ぶ
end
``` pic.twitter.com/E3CU9PQ2Z6
タグ: Julia言語
posted at 23:57:15

#統計
nbviewer.jupyter.org/github/genkuro...
r×cの分割表におけるPearsonのχ²統計量
も更新。期待度数が
1.0 2.0 3.0 4.0
1.5 3.0 4.5 6.0
2.5 5.0 7.5 10.0
と5未満のセルが8個もあるのに、ピアソンのχ²統計量の分布はχ²分布で非常によく近似されています。もちろん周辺度数の全固定はしていない。 pic.twitter.com/miIZzJsBGD
タグ: 統計
posted at 23:46:58

#統計 添付画像は
nbviewer.jupyter.org/gist/genkuroki...
より。すべてのセルの期待値が等しい場合(離散性の問題が最も緩和される場合)
無印Fisher検定さん曰く「Poisson, multinomial, 2 binomialsで45度線に近付かないのはワイだけ」www
n=4×{50,100}ではhypergeometricでχ²,G,mid-p Fisherがほぼ一致。 pic.twitter.com/2yKaxqdNv1
タグ: 統計
posted at 23:35:29

明日の資料として平成30年度試行調査の記述問題部分を整理したファイルを作ったので掲示しておきます。ヒアリングの予習用にお使いください。
www.dropbox.com/s/2xva5fqbsez4...
www.dropbox.com/s/bt0967sxgs0o...
タグ:
posted at 23:30:30

#統計 そして、そのYatesさんの誤解はFisher vs. Pearson論争におけるFisherさんの間違った議論が由来になっているものと思われます。
現代の数学科のセミナーだとFisherさんはフルボッコにされていたと思う。もとの確率分布を勝手に変えて証明しても無意味。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 23:20:01


#統計
nbviewer.jupyter.org/gist/genkuroki...
2×2の分割表の独立性検定
を更新した。そこにあるプロット群を見れば、Fisher検定が採用している超幾何分布モデルが離散性の問題をひどく悪化させることも見えます。2×2の右下の部分がサンプルを超幾何分布で生成した場合。 pic.twitter.com/JxWzrM8tso
タグ: 統計
posted at 23:10:28

@kale_aojiru #超算数 2005年のある調査では、千人の3年生が学年末に、掛順2割正解、割り算の文章題8割という成績を残していますtwitter.com/temmusu_n/stat...。学年が上がるにつれ、掛順はバラツキが出、割り算は習熟すると仮定すれば、相関は逆相関であると予想できます。
タグ: 超算数
posted at 23:07:32

#統計 Yatesさんの論文(1984)からの引用
mathfaculty.fullerton.edu/sbehseta/Yates...
多項分布や2つの二項分布でも離散性の問題があるのに、それを勝手に超幾何分布で置き換えて離散性の問題を大きく悪化させることによって誤差を大きくしたP値の値を報告することにどういう意味があるのやら。 pic.twitter.com/1MB847IEaw
タグ: 統計
posted at 23:05:37
ごまふあざらし(GomahuAzaras @MathSorcerer
Flux が Julia 1.4 で動かない問題がJulia slackで認知されたようなので将来的に解決するでしょう. pic.twitter.com/lH4dediKvC
タグ:
posted at 23:02:45






一応言っておくと、遠山啓自身は「かけ算の順序でバツ」を批判していた。遠山啓自身はそれほど露骨な超算数ではない。
しかし、彼の薫陶を受けた銀林浩など数教協の面々は完全なトンデモ算数に陥っている。
タグ:
posted at 22:20:45

超算数の起源の1人である遠山啓は、「氷が溶けたら何になる?」で「春になる」と答えたらバツになったという真偽不明のエピソードに「ケシカラン採点だ」とdisっていた。
彼らの中では矛盾はないのだろう。彼らには、4人に3個ずつを4×3とするのは、1+1を3とするのと同様の「間違い」なんだろう。 twitter.com/devilGod510/st...
タグ:
posted at 22:19:02

マスメディアの報道だけではない。
文部科学大臣補佐官だった鈴木寛氏ですら、そうした言説を無批判に提示していることは見逃せない。
twitter.com/RochejacMonmo/...
タグ:
posted at 22:12:58

あたかも今の大学入試が記述式を十分に実施していないかのようなミスリードが行われていたことは繰り返し指摘しなければならない。個別二次試験での数理社の記述式問題を丸ごと無視する議論。
英語民間試験導入推進側の「今の大学入試は実質1.5技能」というミスリードと相通じる酷い印象操作だ。 twitter.com/RochejacMonmo/...
タグ:
posted at 22:09:18


非公開
タグ:
posted at xx:xx:xx


非公開
タグ:
posted at xx:xx:xx


`@test_broken` で broken test cases 作っておくと、直した時こんな感じで教えてくれてちょっと嬉しくなるのでオススメです。これとかcodecovとか、ある種のgamificationがあるとまだテスト書く気になる。
#julialang pic.twitter.com/C7i148278c
タグ: julialang
posted at 21:20:16

名著「CPUの創りかた」のTD4のバグが取れて、ようやく完成した。動画は
1. リセットボタン押す
2. 手動クロック(ボタン連打)
3. 自動クロック(1Hz)
4. 自動クロック(10Hz)
の順に実行したところ
* 左の縦4黄LEDがAレジスタ
* 左の縦4白LEDがプログラムカウンタ
* 右の縦白8LEDが機械語
* 下の畑がROM pic.twitter.com/8aoUJ51Wej
タグ:
posted at 20:09:41


ニセ科学商売が巧妙に法的規制から逃れて、野放しとなっている状況について経験をもとに書きました。
消費者が勉強して賢くなることで自衛しなければなりませんが、このままで良いのか…。 twitter.com/kumikokatase/s...
タグ:
posted at 19:41:12

『消費者法ニュース』に寄稿しました。
行政の対応が難しいニセ科学商品の問題…片瀬久美子(サイエンスライター)
www.clnn.net/number/news121...
タグ:
posted at 17:00:16

Non-Newtonian fluids turn into solids when pressure is applied. It's called the process of impact-activated solidification that occurs when compressive forces are applied to fluid-grain suspensions buff.ly/2NB5r7r pic.twitter.com/HbdRhygqXm
タグ:
posted at 16:53:26
ごまふあざらし(GomahuAzaras @MathSorcerer
#Julia言語
プロットなんだけれど
julia -e 'using Makie; plot([1,2,3]) |> display; sleep(10)'
のようにすればREPLを明示的に起動する必要ないことを知った.
display と sleep がいるけれど.
タグ: Julia言語
posted at 16:18:36

@sshojiro_jp Makefileも添付してありますので、Linux環境でのビルドも可能です(Boost C++ Librariesと、Eigenは必要です)。
タグ:
posted at 15:54:32

@sshojiro_jp 動きます。ただしビルドには、Boost C++ Libraries ( www.boost.org )と、Eigen ( eigen.tuxfamily.org/index.php?titl... )が必要です。
タグ:
posted at 15:52:51

非公開
タグ:
posted at xx:xx:xx

水素原子に対するシュレーディンガー方程式のFEMのプログラムが完成しました。以下URLでソースコードを公開しています。下の画像の通り、厳密な波動関数と、極めて近い結果が得られました。基底状態のエネルギー固有値も、厳密な値の99.997%を与えることが分かりました。
github.com/dc1394/hydroge... pic.twitter.com/TMNB7o5GNg
タグ:
posted at 15:46:15

Initial conditions:
Pendulum 1: position: (0.43, 1.50), mass: 3.24
Pendulum 2: position: (-0.85, 2.69), mass: 1.44 pic.twitter.com/Mc619a2zrY
タグ:
posted at 15:20:21

#統計 しまった。上で紹介したFisherさんの間違った議論については以下のリンク先のYatesさんの誤解の指摘の続きとして書いたつもりだった。投稿するスレッドを間違えました。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 14:02:32

#統計 しまった。以下のリンク先以降のスレッドはこちらに書いたつもりだったものです。間違ってしまった。
Yatesさんの誤りが、Fisherさんの誤りに由来していることの確認です。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 14:00:57








#統計 実際にFisherさんが間違った議論(勝手に全周辺度数を固定する議論)をしている証拠を論文
zenodo.org/record/1449484
から引用しておきます(添付画像)。
r×cの分割表の独立性検定のχ²統計量が漸近的に自由度(r-1)(c-1)のχ²分布に従うことを証明するためには、全周辺度数固定は必要ありません。 pic.twitter.com/iPhdw69fxy
タグ: 統計
posted at 13:15:14



動的に統計グラフを作りたいという用途には良いのかも。Hans RoslingのTEDトークを思い出した。 www.ted.com/talks/hans_ros...
タグ:
posted at 12:57:37







ごまふあざらし(GomahuAzaras @MathSorcerer
RT元うちのあざらし多めでワラッタ twitter.com/u_ribo/status/...
タグ:
posted at 12:26:24







学校に入り込むニセ科学 (925) (平凡社新書) 左巻 健男 www.amazon.co.jp/dp/4582859259/... @amazonJPさんから
↑「ただいま予約受付中です。」
*ぜひ、次も一緒に。
暮らしのなかのニセ科学 (平凡社新書) 左巻 健男 www.amazon.co.jp/dp/4582858473/...
タグ:
posted at 12:00:24



#統計 せっかくなので、Fisher vs. Pearson の有名な論争に関するFisherの1922年の論文も確認してみました。
zenodo.org/record/1449484
からダウンロードできます。案の定、Fisherさんも全周辺度数固定という不適切な方法で分割表のχ²検定のχ²分布の自由度を求めていました。
タグ: 統計
posted at 11:51:18





津田和俊/急激に進行した網膜剥離と闘って @kaztsuda
そういや、おれのところにも風疹の抗体検査とワクチン接種のクーポン届いたな。いや、頑張ってる感は伝わっては来るのだけど、風疹が流行れば風疹だけ、麻疹が流行れば麻疹だけ、それも毎回特定の年代や性別に限って、と場当たり的な政策を続けていても、モグラ叩きなんじゃないかなぁ。。。
タグ:
posted at 07:50:16


When you mix a top with a yoyo you get a spin top and when you launch a spin top you are a spin top sniper, a trick shot expert who plays with angular momentum buff.ly/2IDJ8dK pic.twitter.com/cyqF6raGKr
タグ:
posted at 01:00:25


#統計 仮に根拠があったとしても、私の計算結果を見れば分かるように、その根拠は間違っていると思われます。
nbviewer.jupyter.org/gist/genkuroki...
2×2の分割表の独立性検定
独立性の帰無仮説を満たす乱数でサンプルを10万個生成して、各検定でP値がx以下になる割合(確率)をx=0~0.05でプロット。 pic.twitter.com/R3Xxzt3NeV
タグ: 統計
posted at 00:12:42


#統計 #Julia言語
期待値が
a 10-a
2a 20-2a
の場合。
G-test → G検定(対数尤度比検定)
Fisher → Fisher検定
mid-p F → mid-p版Fisher検定
全般的にχ²検定がrobustでいい感じだと思う。
度数が5以下のセルがあってもこの程度の誤差ならいいんじゃないかな。Fisher検定よりはずっとよい。 pic.twitter.com/EqPgGhmeun
posted at 00:05:17