黒木玄 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 を見て下さい。
2017年09月19日(火)
@h_okumura #数楽
RとJulia言語で
oku.edu.mie-u.ac.jp/~okumura/stat/...
Fisherの正確検定かカイ2乗検定か
の結果(のsubset)を再現してみました。再現結果は以下のリンク先にあります。
nbviewer.jupyter.org/gist/genkuroki...
タグ： 数楽
posted at 00:16:03
@h_okumura #数楽 Rによってプロットした
fisher.test(correct=TRUE)
の場合のグラフが
oku.edu.mie-u.ac.jp/~okumura/stat/...
と
nbviewer.jupyter.org/gist/genkuroki...
の一番下で全然違う形になってしまっています。続く
タグ： 数楽
posted at 00:18:27
@h_okumura #数楽 続き
Julia言語からRを呼び出してテストしたことが原因のトラブルでないことも確認してみました。
gist.github.com/genkuroki/5aa1...
もしかして、chisq.test函数の仕様が変わったのでしょうか？
タグ： 数楽
posted at 00:20:34
#数楽 直接 Gist の方を見たい人は次のリンク先を見て下さい。確実に最新版を見たい人はGistを直接見た方が安全です。
gist.github.com/genkuroki/7694...
2×2の分割表に関するカイ二乗検定とFisherの正確検定の比較
タグ： 数楽
posted at 00:22:22
#数楽 #JuliaLang 以下ではJulia言語で2×2の分割表のカイ二乗検定とFisher検定のp値と求める函数を書いて、RCallパッケージを利用してR言語の結果と比較してバグっていないことを確認したりしています。
nbviewer.jupyter.org/gist/genkuroki...
posted at 00:28:52
#数楽 #JuliaLang 既存の統計のパッケージがあるのにわざわざ自分で函数を書いたのは、自分の勉強のため。パッケージをブラックボックスとして使うのは怖い。一度、自分で実装してみて、感じをつかんでおきたい。既存のパッケージもどこまで信用できるのかは不明だし。
posted at 00:30:31
#数楽 #統計 #JuliaLang 私が作ったグラフは添付画像の通り。
oku.edu.mie-u.ac.jp/~okumura/stat/...
の結果のsubsetがきちんと再現されています。 pic.twitter.com/aSwNDIgcoR
posted at 00:33:03
@genkuroki 例えば先ほど数回呟いたポアンカレ円盤の作図の話とかが図版混じりでサクサクっと描けるでしょうか、、、まぁ何をもって「サクサクっと」かってところが問題だけど。こんな話はGeoGebraみたいなのだと話が早そうだけどあれはあれで何やら無理から特化してて、、、
タグ：
posted at 00:43:47
@genkuroki ええ、まぁそう意味でなら何でも描けるのだとは思います。作図的な操作のうち、点を取る、直線／線分が引ける、2点の長さを取れる、円が描ける、あとは直線や円の交点が取れるかな。サクサクっと出来そうですか？ #自分で試せ
タグ： 自分で試せ
posted at 00:51:20
@genkuroki それ、何だか僕なんかは堪え性や持続性が弱ってるのか振り回される事や気になる事を抱えすぎてるのか、なかなか行動に起こせないし、起こしても続かないので困ります。誰か早よ♬
タグ：
posted at 01:27:40
#数楽 #統計 #JuliaLang ついでにG検定の場合も付け加えてみました。
gist.github.com/genkuroki/6d90... pic.twitter.com/Y1RtcIFnC5
posted at 02:29:50
#数楽 #統計 #JuliaLang 2×2の分割表でのカイ二乗検定とG検定を比較したい人は多いと思うので、p値が小さい部分のグラフも作成しました。 pic.twitter.com/dsbYsS3SGS
posted at 02:30:56
#数楽 #統計 #JuliaLang このグラフを見る限りにおいては、2×2の分割表では古典的な(補正無しの)カイ二乗検定で十分なんじゃないかなと思いました。
「このグラフ以外を見たらそうとは言えない」ということを知っている人がいたら、計算例を見せて下さい。
posted at 02:33:14
#数楽 #統計 #JuliaLang 公開しているJupyter notebook
gist.github.com/genkuroki/6d90...
のRCallを使っている部分以外は juliabox.com で動くと思います。誰にでも再検証可能。
posted at 02:40:26
#数楽 #統計 #JuliaLang 私が公開しているJupyter notebookを手元のパソコンで使いたい場合には次のリンク先を参考にして下さい。
gist.github.com/genkuroki/81de...
posted at 02:41:34
#数楽 #統計 #JuliaLang もとの試行回数は10^5だったのですが、以下のリンク先では10^7に増えています。計算時間は約2分。最適化も並列化も全然していません。気軽に"n"を増やせることはJulia言語の利点。
gist.github.com/genkuroki/6d90...
posted at 02:45:27
@Arakur65536 @sekibunnteisuu 「掛け算順序固定肯定派の意見」であなたが尊重するべきだと思っているものを引用して頂けませんか？あなたの想像上の意見ではなく、引用であることが重要なポイントです。
私はこの件について過去数年間継続的に注意を払っているのですが、尊重するべき意見を述べた人は皆無でした。 #超算数
タグ： 超算数
posted at 03:14:44
@Arakur65536 @sekibunnteisuu 今までの経験では、掛算順序固定強制支持の立場で尊重するべき意見を述べた人は皆無であり、過去数十年間分の文献を調べても尊重するべき意見はどこにも載っていませんでした。
まだ発見していない尊重するべき意見があるのであれば是非とも教えて欲しいものです。 #超算数
タグ： 超算数
posted at 03:17:10
@Arakur65536 @sekibunnteisuu 過去数年間、掛算順序固定強制支持の立場でどのようなことが言われているかについてほとんど知らないままで、勝手な想像で意見を述べる人達が次から次へと現われました。
そうならないことを祈っています。 #超算数
あなたが知っていることの証拠を見せて下さい。
タグ： 超算数
posted at 03:19:56
あれま、見てる分には面白い人だったのに…
ノイホイこと菅野完氏、Twitterアカウントが永久凍結される。怒りと嘆きと擁護の反応をまとめてみた。 - Togetterまとめ togetter.com/li/1152295 @togetter_jpさんから
タグ：
posted at 06:18:28
Smithsonian Magazine @SmithsonianMag
A new movie sets its doomed entrepreneurs amidst 17th-century “tulipmania.” smithmag.co/xXDtQo
タグ：
posted at 08:45:10
V745 Scoの新星爆発(チャンドラX線天文台)
chandra.harvard.edu/photo/2017/v745/
1937/1989/2014に観測。爆発後のガスの動きからこの種の新星爆発は白色矮星の両極方向に広がるというモデル作ってる
www.youtube.com/watch?v=MJZU15...
タグ：
posted at 08:51:42
「17世紀欧州ではチューリップが流行してオランダでは屋敷の価格よりもチューリップが畑の地価が上がるほど投機的に市場は盛り上がったけど、バブルが弾けてトレーダー連中は破滅した」
と言われているが、
「これの話は事実ではない」
まじか twitter.com/SmithsonianMag...
タグ：
posted at 09:10:35
「チューリップマニア」という書籍の紹介
press.uchicago.edu/ucp/books/book...
当時の貴族や商人がどれくらいチューリップ投機に金使ったとか検証してるらしい。
考えてみれば資料が残っている時代の話だから、市場経済の規模とか分析できる
タグ：
posted at 09:22:28
「学術会議」と一括りにできないことは知っておいたほうがいい。これまでも「第一部」はむしろ放射能デマを積極的に擁護してきた twitter.com/nonbeekumasan/...
タグ：
posted at 10:19:17
「学術会議9.12報告」がstudy2007なるどこの誰とも分からない人物が書いた「見捨てられた初期被曝」という無査読の一般向け書籍を参考文献に挙げたのはこの報告書の信頼性を著しく損なっているということです。なんだそれ、です twitter.com/nonbeekumasan/...
タグ：
posted at 10:28:57
非公開
タグ：
posted at xx:xx:xx
広く読まれてほしいと思います。／【SYNODOS】「福島で次世代に放射線被曝の影響は考えられない」ということ――日本学術会議の「合意」を読みとく／服部美咲 / フリーライター synodos.jp/society/20472
タグ：
posted at 11:33:32
福島に住む人を悩ませてきたことの一つが「いろんなことを言う人（専門家）がいる」であったと思う。『「両論併記」はすでに公平ではない』ことをメディアはどう受け止めるのか。ウケる話題や即時性が求められることだけでなく、こうした事実をしっかりと伝えるための役割があるのではないでしょうか。 twitter.com/hangui_/status...
タグ：
posted at 11:37:49
@Arakur65536 @genkuroki 読んだ限り、この人自身が教えいているようでもないし、これまでの議論を踏めている風でもなく、頭の中で「こうじゃないか？」と想像しているだけですね。周回遅れで論ずるに値しません。
タグ：
posted at 12:47:38
@Arakur65536 @genkuroki この方が「こういう理由で順序指導をしているのでは」ということを、そもそも順序指導をしている人は述べていません。
タグ：
posted at 12:48:33
両論併記は公平ではないという話を僕がプレイボーイだかSightだかのインタビューで話したのは震災の1年後くらいだったと思うけど、事態はあまり改善していない twitter.com/hangui_/status...
タグ：
posted at 12:58:35
#超算数 sukyokensenrichuo.blog.fc2.com/blog-entry-127...
＞学校では、縦横のどちらか基準を決めるようですが、数教研では子どもが何を１あたりと考えたのかを説明できれば正解としています。というのも基準となるモノは変わることもあり、自分の考えを自らの言葉で説明できるなら理解していると
タグ： 超算数
posted at 14:13:22
ご紹介するのは「ナカトジール」という中綴じが簡単に出来る定規。この大きさでA3や新聞紙まで中綴じ出来るのには驚き
一緒にシンプルなホチキスも入荷。この山吹色はなんだかかわいい
身近にあるものを無駄に冊子化したくなりますよ
■ナカトジール ¥380税別
■ホチキス ¥360税別 pic.twitter.com/yB9MSiCotR
タグ：
posted at 14:24:13
#超算数 sukyokensenrichuo.blog.fc2.com/category5-9.html
＞しかし、問題に書いてある数字の順番どおりにかけばよいと思い込んでいるため、６×７という式を書いてしまうケースが少なくありません。結果として答えの４２個は合っていますが、式が間違えているということになります。
タグ： 超算数
posted at 14:24:23
#超算数 sukyokensenrichuo.blog.fc2.com/blog-entry-661...
タイルを長方形状に並べて、縦が１あたり、横がいくつ分らしい。それで、３×４と４×３で意味が違うとか言っているけど、長方形状に並べた段階で、１あたりもいくつ分も区別不要、となるのが普通では？
タグ： 超算数
posted at 14:34:24
@Ra_koyama １あたり と いくつ分 を区別できると思っていると言うことは、その人自身が理解していないんですよね。水道方式・数教協の場合、掛け算を累加ではなく１あたり×いくつ分、とするので、主流派以上に意味に拘り順序に拘りがちになる。
タグ：
posted at 14:42:17
メガハウス「カードではじめて将棋」、10月下旬発売
www.megahouse.co.jp/megatoy/produc...
盤には初期配置の駒が薄くプリント、コマには進み方が表示。紙製だから盤やコマが大きくても手頃価格（900円） pic.twitter.com/oT2hdbOHdF
タグ：
posted at 15:07:24
@sekibunnteisuu @Arakur65536 Arakur65536さんへ
論外。
引用するときには、著作者を明記しないとダメ。
単なる思い込みで言っているだけの発言なのか、そうでないかを第三者が検証可能なだけ十分な情報を示さない単なる転載を議論の材料にしてもらえると勘違いしないように。
タグ：
posted at 15:12:27
@sekibunnteisuu @Arakur65536 私や積分定数さんを含むこの話題への数年前からの参加者達の多くは、社会的に算数教育の専門家とみなされる人達が書いたものを多数引用して議論しています。ツイッターの過去ログにその記録がすべて残っています。
タグ：
posted at 15:17:27
@sekibunnteisuu @Arakur65536 インターネット上の発言を引用するときにはリンクもはっておけば、この話題に注目している人達にとって親切です。
毎回きちんとリンク先などの情報を示すのは大変なので手抜くことがあってもいいと思う。
しかし、著者名を塗り潰す手間をかけたのはナンセンス。論外。
#超算数
タグ： 超算数
posted at 15:29:28
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#超算数 sukyokensenrichuo.blog.fc2.com/blog-entry-457...
＞答えはあっているのですが、解き方がしっかりと書ききれていません。
まず、お金は全部でいくらもっているのかという式を書かなくてはいけません。
タグ： 超算数
posted at 16:50:35
#超算数 水道方式を推奨しているこのブログ主は、１／３個で１００ｇだと1個では？を１００×３とするのは誤りで、１００÷１／３が正しい、と言っている。
selfyoji.blog28.fc2.com/blog-entry-476...
タグ： 超算数
posted at 16:54:10
１本a円の鉛筆２０本分の値段も１本２０円の鉛筆a本の値段もどっちも20aである。その文字式のルール（数字が前で文字は後で辞書順、＊は省略する）を適用する前の式はそれぞれ20*aとa*20で違うはず、なんてことは意識に上がりません。
#超算数 #掛算 twitter.com/sekibunnteisuu...
posted at 16:59:12
#超算数 タイルの縦横で 1あたり と いくつ分 を区別することは、頻繁に教えるらしい。
sukyokensenrichuo.blog.fc2.com/blog-entry-419...
タグ： 超算数
posted at 16:59:19
#超算数 ずつ に着目
sukyokensenrichuo.blog.fc2.com/blog-entry-311...
＞４×５＝２０ としてしまうお子様も結構多くいらっしゃいます。
しかし、問題をよく読んでみると「1本に花が５つずつ」とあります。つまり１あたり量がかかれているわけです。このことから正しい式は ５×４＝２０
タグ： 超算数
posted at 17:06:43
@AS_Insects 「割合が難しい」というのは、難しく教えられているのが原因だと思います。基本的には比です。たまたま2つの量が同じ種類と言うだけのこと。
水４８０ｇと食塩２０ｇを混ぜて食塩水５００ｇを作る。この食塩水１００ｇ中には食塩水何ｇ解けているか？
これが分かる子は割合の概念を理解している
タグ：
posted at 17:15:22
@AS_Insects 速さの概念を理解している状態では、速さは単位時間あたりの距離、というのがおまけみたいな物。
同様に、１００ｇ中の食塩が４ｇと分かる子には、割とか％などはおまけみたいな物。
タグ：
posted at 17:19:41
いつものどんぎつねとは違うキツネがいますね。
なにやら”りょうり”の作り方に悩んでるようで……！？
#けものフレンズ #どん兵衛 pic.twitter.com/aOuMcpZ9Jz
posted at 18:05:17
@genkuroki @sekibunnteisuu 不適切な転載を行ってしまい申し訳ありませんでした。このような不貞者にもはやこの問題を議論する資格などないでしょう。全ての前言を撤回し、この無学のために多くの時間を費やして下さった方々に謝罪します。
タグ：
posted at 18:11:33
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
コーシー分布は不機嫌な分布であるらしく、シミュレーションから作った分布と全く同じ分布をMCMCでフィッティングしてもめちゃめちゃ発散する。無論ReParameterize込みでやっているんだけれど。。
タグ：
posted at 19:45:21
作らせていただきましたー、どん兵衛はマイ主力夜食なので楽しかったです！
（あと何気に日清さんアニメ関わり、10年ぶり２度目で嬉し…） twitter.com/donbei_jp/stat...
タグ：
posted at 19:50:17
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
学年別漢字配当表の「標準」字体について⇒〈「本当はどちらでもよいのだが、負担を減らすために標準を示しただけ」という経緯と趣旨が、まったくと言ってよいほど周知されなかった〉toyokeizai.net/articles/-/188...
タグ：
posted at 21:18:41
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@InfinityLoopGIF @genkuroki 黄色い球体以外が今一つ回転しているように見えないのは底面のナミナミが固定されているからでしょうね。
タグ：
posted at 22:15:52
クウェート大使館。丹下健三作品。ここも来るの初めて。来春建て替え工事始まる前に見ておこうと思って。空中庭園に架けられた階段がとてもいい。天空の城感を強調している。無くなってしまうのは本当に残念だ。 pic.twitter.com/2vyAQgaH9u
タグ：
posted at 22:27:48
非公開
タグ：
posted at xx:xx:xx
#数楽 #統計
2×2の分割表の各行各列の独立性に関するカイ二乗検定、G検定、Fisherの正確検定の比較をシステマティックにやってみました。 #JuliaLang のJupyter notebookがリンク先にあります。
gist.github.com/genkuroki/32a8...
posted at 22:34:05
#数楽 #統計
#JuliaLang を使って行った計算は以下の通りです。
(1) 各行各列が独立な確率分布を用意し、それを用いて10万個のサンプルを生成。
(2) 各々のサンプルについてp値を計算し、p値がx以下になる確率(割合)を求めてプロット。
続く
posted at 22:43:35
#数楽 #統計 #JuliaLang
(3) p値の求め方は、カイ二乗検定、G検定、Fisherの正確検定の3通りあるので、そのプロットによって、それらの3つを比較することができる。
(4) p値がx以下になる確率は理想的にはxに一致するべき。誤差は小さい方がよいと考える。続く
posted at 22:45:14
#数楽 #統計 #JuliaLang リンク先を見ればわかるようにたくさんの場合をシステマティックにプロットしているのですが、p値がx以下になる確率がxに最も近くなるのがカイ二乗検定になる場合が大部分でした。
gist.github.com/genkuroki/32a8...
posted at 22:47:52
#数楽 #統計 奥村さんによるconditionalでない場合の計算を見てもわかるように、Fisherの正確検定ではp値がx以下になる確率はxよりもずっと小さくなります。「誤差」は非常に大きい。
oku.edu.mie-u.ac.jp/~okumura/stat/...
posted at 22:50:01
#数楽 #統計 続き。リンク先のたくさんのグラフを眺めてみればわかることですが、G検定よりもカイ二乗検定の方が「p値がx以下になる確率とxの違い」の小ささでは圧倒しています(小さいほどよいはず)。
gist.github.com/genkuroki/32a8...
posted at 23:02:18
#数楽 #統計 #JuliaLang
gist.github.com/genkuroki/32a8...
2×2の分割表のカイ二乗検定とG検定とFisherの正確検定の比較
をついさっき更新しました。もうダウンロードしてしまった気の早い人は新しい版に入れ替えて下さい。
posted at 23:16:14
#数楽 #統計 #JuliaLang
再検証をRでやってくれる人がいると特にうれしいです(特にG検定)。
oku.edu.mie-u.ac.jp/~okumura/stat/...
のコードをそのまま使うと私によるJulia言語による計算よりも20倍ほど時間がかかるはずです。その点は注意と工夫が必要。
posted at 23:19:29
#数楽 #統計 #JuliaLang すでに述べたことですが、自分でG検定を実装する場合には、x log(x/y) は x→0 で 0 になることに注意しなければいけません。Julia言語では0*log(0)=0*(-Inf)=NaNとなってエラーが出ないのではまった。
posted at 23:23:10
#数楽 #統計 #JuliaLang Gistに直接アクセスすると読めるようになるまで少し時間がかかるので、それが嫌な人は
nbviewer.jupyter.org/gist/genkuroki...
2×2の分割表のカイ二乗検定とG検定とFisherの正確検定の比較
の方を見て下さい。
posted at 23:26:00
#数楽 #統計 #JuliaLang プロット結果の例
左半分は n, p に対する多項分布で生成されたサンプルのp値がx以下になる確率のプロットです。カイ二乗検定の太実線が最も45度線に近いです。プロットの範囲はx=0～0.1です。有意水準5%に興味がある人向け。 pic.twitter.com/8PZpIOAQ9o
posted at 23:30:07
#数楽 #統計 #JuliaLang 脱線：「検定」に対応する英語は test です。日本語でも「検定」と言わずに「テスト」と言えば権威的な響きが薄れていいかも。数学用語で片仮名語にした方がわかりやすい場合は少なくないと思う。「行列の階数」は「行列のランク」と言わないと通じ難い。
posted at 23:33:08
非公開
タグ：
posted at xx:xx:xx
#数楽 #統計 #JuliaLang 上で示したプロット結果のサンプルと同じ場合を0から1まで全部プロットすると、添付画像のようにカイ二乗検定とG検定の45度線の区別がつかなくなります。拡大して眺めないとカイ二乗検定とG検定の違いがわからない。 pic.twitter.com/liWBItPJkG
posted at 23:36:38
#数楽 #統計 別の例
添付画像1
0～1の全体のプロット。
左半分を見るとカイ二乗検定とG検定の区別がつかないが、右半分を見るとG検定が小さなxで45度線から大きくずれていることがわかる。
添付画像2
x=0～0.1でG検定のp値がx以下になる確率はxよりかなり大きい。 pic.twitter.com/m4fkEler8T
posted at 23:52:14
