黒木玄 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 を見て下さい。
2020年06月17日(水)
#統計 R言語のbinom.testがデフォルトで表示するClopper-Pearsonの信頼区間は同時に表示されるp値と整合性がありません。
パラメーターpが95%信頼区間を外れているのに、p値が5%を超える場合が出て来てしまう。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki...
続く
twitter.com/kenmomd/status... pic.twitter.com/lB6wjNpat8
タグ: 統計
posted at 00:18:00
#統計 せっかくなので、データから決まるパラメーターの函数としてのP値函数もプロットしてみました。
パラ―メータを変えたときのP値の値がどうなるかが全部分かる。95%信頼区間はP値函数が5%以上になるパラメーターの範囲です。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/1i3uK5PAoz
タグ: 統計
posted at 00:27:06
#統計 #Julia言語
P値函数がプロットしてあれば、有意水準の横線を引くことによって、任意の信頼係数の信頼区間の様子がわかります。
特定のパラメーターのP値や信頼区間はP値函数が持つ豊富な情報の一部分を切り取ったものに過ぎません。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/BQgR5lwGFK
posted at 00:29:43
HPVワクチン「薬害」報道の経緯を解説した論文。朝日新聞から始まったことが明記されている。 academic.oup.com/cid/article/63...
タグ:
posted at 00:54:00
#超算数
妙な反応が出易い話題だと思いますが、12×3で文句無しに正解です。イチャモンつける奴らは馬鹿扱いで問題ない。
私のうちで似たようなことがあった場合には、「こういう場合には学校の先生からバツをつけられても気にしないこと!うちで全部マルをつけます」としています。 twitter.com/akkii_206/stat...
タグ: 超算数
posted at 01:15:22
julia> 1[1][1][1]
1
Found this on @JuliaLanguage slack!
speechless. <3
#julialang
タグ: julialang
posted at 01:19:09
#統計 #Julia言語 最初の話題に戻る。
P値の側の定義を「片側確率の2倍」に変えれば、P値とClopper-Pearsonの信頼区間が整合的であるようにできます。
function pval_dos(p, n, k)
bin = Binomial(n, p)
min(1, 2cdf(bin, k), 2ccdf(bin, k-1))
end
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/RX3RDLpZdo
posted at 02:15:46
"Maybe this could use a Web UI"
Regular Person: Learn JavaScript
Me: I wonder how hard it would be to hack #julialang into a browser engine as an alternative language until wasm dom bindings get good enough.
タグ: julialang
posted at 06:07:52
業者団体もどんどん信用なくなると、「次亜塩素酸水」自体だめになるよ👶まともに「次亜塩素酸水」に取り組んでいる人は声を上げて噴霧などはやめさせるべきでしょうね。全部だめになるよ本当に。業界ごと。 twitter.com/minesoh/status...
タグ:
posted at 07:29:27
#統計 #Julia言語 最初の話題に戻る。
P値の側の定義を「片側確率の2倍」に変えれば、P値とClopper-Pearsonの信頼区間が整合的であるようにできます。
function pval_dos(p, n, k)
bin = Binomial(n, p)
min(1, 2cdf(bin, k), 2ccdf(bin, k-1))
end pic.twitter.com/vjYE8n4Fkz
posted at 09:43:15
#統計 #RStats
すでに指摘したように、R言語のbinom.testが表示するP値と信頼区間のあいだには整合性がない。
実はさらに、fisher.testが表示するP値と信頼区間のあいだにも整合性がないです。詳しくは以下のリンク先を参照。
twitter.com/search?q=from%...
twitter.com/genkuroki/stat...
posted at 09:53:49
この辺りよく解ってなくて解りたい感がある。Rでの信頼区間とp>0.05の区間が合わないのって統計量(二項分布)がsymmetricじゃないのに片側確率の2倍で計算してるからってことなのかしら。 twitter.com/genkuroki/stat...
タグ:
posted at 10:11:56
#統計 Rのbinom.testやfisher.testでは、P値と信頼区間の計算で、異なるP値函数を採用しているので、表示されるP値と信頼区間の整合性がなくなっています。
binom.testやfisher.testが表示する信頼区間はP値を「片側確率の2倍」で定義した場合の信頼区間になっています。
twitter.com/Isa_rentacs/st... pic.twitter.com/uTCy8qYGu6
タグ: 統計
posted at 10:35:43
ポインタのインクリメントの増分が値型のサイズに依存するのは、ポインタをイテレータとして考えると自然だと思う。 twitter.com/dqneo/status/1...
タグ:
posted at 11:17:52
Mark Kittisopikul ht @markkitti
@Sov_trotter @JuliaLanguage julia> 1(2)(3)(4)
24
julia> 1/2π
0.15915494309189535
julia> exp(2(im)π)
1.0 - 2.4492935982947064e-16im
タグ:
posted at 11:38:11
Slack has a little "you're about to notify 6000 people in 24 time zones" warning. I feel like once you get above a certain threshold of timezones or people it should give you a special message. (This is the #julialang slack BTW - and yes we limit the power to rare announcements).
タグ: julialang
posted at 12:02:40
#JuliaLang ranks at 33 in the latest Tiobe. That is 15 spots higher than last month. While the index is quite noisy in that range, the rise is correlated with a huge influx of users into the #julialang community in the last few months.
www.tiobe.com/tiobe-index/
タグ: JuliaLang
posted at 12:34:24
非公開
タグ:
posted at xx:xx:xx
Mathematica で図がはげる.
PlotPointsをあげるとフリーズする.
さぁ,どうしよう. pic.twitter.com/v1HdKMEfW8
タグ:
posted at 12:49:37
未曾有の事態はビジネスチャンス
貴重なデータが取れそうですね
人命より読解力のほうが大切ですか
>今年継続してRSTを受検することは極めて重要
>格差が広がっていないかを検証できるのは今年をおいて他にない
>第二波
学校を休校すべきか
否かの判断をするひとつの指標
2020.6.17 #新井紀子 pic.twitter.com/3NdvBIrM07
タグ: 新井紀子
posted at 13:15:15
#統計 実際に、二項分布モデル(より正確にはベルヌイ分布モデル)の場合に、頻度論版のP値函数(正規分布近似で定義)とベイズ版のP値函数(事後分布(posterior)を使って定義)を重ねてプロットして作った動画が以下のリンク先にあります。何度見直してもよく一致している(笑)。
twitter.com/genkuroki/stat...
タグ: 統計
posted at 13:32:30
ごまふあざらし(GomahuAzaras @MathSorcerer
キューーーー(美味しそうな叫び声) twitter.com/genkuroki/stat...
タグ:
posted at 13:53:16
95%信頼区間と棄却域の話いまいち腑に落ちなかったので確率関数を書いてみた(縦線は平均)、黄と赤の外側の範囲がα=0.05での棄却域で、紫と茶が95%信頼区間を外側の累積*2=αとして計算したとき、のはず pic.twitter.com/eWOpcsHe9a
タグ:
posted at 13:58:43
このままだったら本当に手遅れになる
人間を賢くする方が、いまわたしが割ける力
リーディングスキルテストというのが
自分の研究人生で最後の仕事 2019.8.21(動画有)
www.asahi.com/articles/ASM7K...
ビジネスチャンスですね
>ドキュメントを自力で読めることは生き抜く上で必須
2020.6.17 #新井紀子 pic.twitter.com/z60UEUYJjE
タグ: 新井紀子
posted at 14:10:27
非公開
タグ:
posted at xx:xx:xx
この件は
国立情報学研究所の見解を聞きたいです
>新井紀子さんの良い側面だけを見て
新井さんを持ち上げている人達が
以上で指摘した問題を無視すれば
結果的に教育に破滅的な悪影響が生じる恐れがある
(黒木氏)
twitter.com/genkuroki/stat...
タグ:
posted at 14:39:03
答えは合っているのに、重大な間違いがあるかも知れないというのは、そもそも大きな矛盾。
試験の解答からわかることは限定的なので、その行間を読むような余計なことをしてはいけない。 twitter.com/HIRO_8_ORIH/st...
タグ:
posted at 15:46:29
#Julia言語 #JuliaLang
Mandelbrot sets by Pluto.jl
gist.github.com/genkuroki/cdad...
github.com/fonsp/Pluto.jl
github.com/fonsp/PlutoUI.jl pic.twitter.com/ObYlHx2vrp
posted at 16:01:23
@HIRO_8_ORIH 教えていないことでも、それに利点があれば評価すべきですね。教わっていないことをできるようにならなければ、社会や科学の発展はありえません。
それを評価できない会社はいずれ先細りでしょうね。
タグ:
posted at 16:28:37
#統計 関連
モデル bin=Binomial(n, p) とデータが「n回中k回成功」のP値を
sum(pdf(bin, j) for j in 0:n if pdf(bin, j) ≦ pdf(bin, k))
とするか
min(1, 2cdf(bin, k), 2ccdf(bin, k-1))
とするかの違い。それぞれが与える信頼区間は異なる。
twitter.com/isa_rentacs/st...
タグ: 統計
posted at 16:30:36
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
ガテン٩( ᐛ )و地道にマスクと手洗い @RrK9xKwM96pOWu6
@TOTB1984 HPVワクチンが出来たと勤務先のクリニックでMRさんから聞いた時、とても嬉しかった。
自分が細胞異形成でがんセンターに通ったから、これで若い娘さん達があの恐怖から逃れられると感謝した。
タグ:
posted at 16:54:58
非公開
タグ:
posted at xx:xx:xx
コロナでもそうでしたが、どうしてこう、医療を右と左に分けたがるんでしょうね。ちなみに私はコロナ禍の中、右翼・左翼両方のリストに入れられまくっております。何をもって右左認定してるのか意味不明です。
#コロナでイデオロギーはほんとアホらしい twitter.com/norichikahorie...
posted at 17:15:45
【オンライン集中講義】中島啓さん「3次元超対称性ゲージ理論のクーロン枝の数学的な定義」
2020年6月22~26日:Zoomによるオンライン授業
www.math.titech.ac.jp/~jimu/Lecture/...
www.math.titech.ac.jp/~jimu/Syllabus...
タグ:
posted at 17:56:05
@EARL_Med_Tw 安定の岩田健太郎
#コロナでイデオロギーはほんとアホらしい pic.twitter.com/BXHkulQP1n
posted at 17:59:52
@GodotOBrien @ROBOmito 「いろんな考え方があっていいから、かけ算に順番があってもいい」と考える人はけっこういるみたいですね。
でもこれ「教えた順番以外は×」という「いろんな考え方」とはまっこうから対立する自己矛盾を含んでいるわけですが…
タグ:
posted at 18:11:10
「ウサギ3羽の耳の数は全部でいくつですか?」と問題文に書いてあれば、それはその通りなので、かけ算の順番によって「3本耳のウサギが2羽いること」になることなど絶対にない。国語と算数の解釈が完全に間違っとる。
通常のウサギの耳は2本なので、それ以外の異常な状態をなぜ考える必要がある? twitter.com/miyuki_MathT/s...
タグ:
posted at 18:26:38
@osamu_takeuchi それだと手間が増えそうな予感。
最初に画像の表示サイズなどの試行錯誤が必要な感じ。
スクリプト一発でできて欲しいのです。
しかもWindowsで。(Cygwinはあり)
タグ:
posted at 19:02:50
@genkuroki @osamu_takeuchi 元が画像なら imagemagick の montage を使ってタイル画像を作ってからPDFに変換するとか. imagemagick.org/Usage/montage/...
タグ:
posted at 19:13:31
@genkuroki ブラウザだと画像の表示サイズやマージン等の調整結果をリロードだけで確認できるし、印刷時に倍率を変えることもできるしで、コマンドラインから直接 PDF を作成するよりずっと簡単に試行錯誤できますが、javascript + css の知識が必要になるところがネックですね。
タグ:
posted at 19:16:49
@akinori_ito @osamu_takeuchi helpx.adobe.com/jp/acrobat/kb/...
はまさに今日やっていた作業そのもの。
19回繰り返した(泣)
マウスを触らずに
jpegファイル達→1ページに2×4枚のPDF
の変換をしたい。
ファイルサイズは1ページあたり200KB程度に抑えたい。
タグ:
posted at 19:28:24
@genkuroki @osamu_takeuchi 元がPDFにできているなら,こういうのもあります polamjag.hatenablog.jp/entry/2015/02/...
タグ:
posted at 19:31:11
@genkuroki @osamu_takeuchi あとは imagemagick convert をつかうといちいちghostscript を呼ばなくても
% magick convert image.jpeg image.pdf
で変換できます
タグ:
posted at 19:33:50
試しにベネッセ全国実力診断テストをやってもらったのだが、算数のこの問題だけ間違えていた。あとは全部正解で国語は満点。
これが噂の「かける数とかけられる数」問題か。
②が正解なのは分かるが①は完全に間違いといえるのか? pic.twitter.com/YQwXrjIOsk
タグ:
posted at 19:51:34
これは誤り。順序擁護論者には、ゴリゴリのネトウヨもいるが、中核シンパもいる。
日教組の教研集会で報告者となることの多い、数学教育協議会は掛け算順序推進。森毅・銀林浩が典型。
掛け算に順序に左右はない。a×bもb×aも同じ。 twitter.com/phoquesan/stat...
タグ:
posted at 20:01:47
@akinori_ito @osamu_takeuchi いやあ、参考になりました。どうもありがとうございます。
ImageMagickで
magick montage hoge001.jpg … hoge008.jpg -tile -tile 2x4 -geometry 500x350+1+1 page001.jpg
のような感じで複数の写真を1枚にまとめることができ、
magick page001.jpg page002.jpg page003.jpg out.pdf
でOK.
タグ:
posted at 21:09:52
@HIRO_8_ORIH 先生や上司が言うことに疑問を持つことも、アリだと思います。大人の言うことが常に正しいわけではないので。
そもそも、自然科学の本質と自分にメリットがあるかどうかは関係がありません。
算数や理科の問題の解答が誰かの都合でころころ変わってはいけないんですよ。
タグ:
posted at 21:27:41
Sphericons are unique solids that roll in such a way that every point on their surface comes in contact with the plane- following the path shown here with white paper. Solids from the sphericon family all have one side and two edges buff.ly/2LHRh4x pic.twitter.com/8qxrofBto6
タグ:
posted at 21:30:16
反ワクチン勢力が存在するせいで日本政府が現在進行形でやっていることがコレ。ひどい怠慢。ウイルスが原因だと分かっていて、ワクチンで予防できることがわかっているのに、ワクチンを積極的に勧奨しないせいで、日本では大勢の若い女性の命も健康も守ることができないでいる。>RT pic.twitter.com/seHxRpWRK6
タグ:
posted at 21:41:20
#超算数 「小島順氏」という題名で掲示板に書き込みました。
8254.teacup.com/kakezannojunjo...
拗らせた人の大作です。
軽い気持ちで見てはいけません。
掛け算の順序について I ー 離散量からはじめる ー
math.oshirase.com/wp-content/upl...
掛け算の順序について II ー 連続量の場合 ー
math.oshirase.com/wp-content/upl...
タグ: 超算数
posted at 21:41:44
#オンライン授業
クォリティーを落としてpdfのサイズを下げたいならば、
magick page001.jpg page002.jpg page003.jpg -quality 70 out70.pdf
のようにすればよい。
twitter.com/genkuroki/stat...
タグ: オンライン授業
posted at 21:41:52
@HIRO_8_ORIH 授業態度なら、試験とは別に「平常点」として採点すればいいだけだと思います。
授業態度でテストの点数に手心を加えると、かえって問題になるでしょう。
タグ:
posted at 21:46:56
#Julia言語 Juliaの文字列はimmutable:
docs.julialang.org/en/v1/manual/s...
【As in Java, strings are immutable: the value of an AbstractString object cannot be changed. To construct a different string value, you construct a new string from parts of other strings.】
Arrayに変換すれば変更可。
タグ: Julia言語
posted at 21:48:31
#Julia言語
String と Array{Char,1} の相互変換
s::String -> collect(s)::Array{Char,1}
c::Array{Char,1} -> String(c)::String
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/ldusAqcEkS
タグ: Julia言語
posted at 21:48:32
OokuboTact 大久保中二病中年 @OokuboTact
@takusansu 掛け順関連ではないですが、小島順氏については、資料をアップしたことがあります。
8254.teacup.com/kakezannojunjo...
タグ:
posted at 21:49:06
@OokuboTact @takusansu #超算数 「算数は難しい」という言説を確固たるものにすることに貢献した人達の一人。どうして、「シンプルで易しい」としなかったのやら(呆)
小島順さんは谷村省吾さんが最近になって権威を引き上げてしまった感じで、谷村さんのファンが擁護する可能性を心にとめておいた方がよいかもしれません。
タグ: 超算数
posted at 22:06:02
@OokuboTact @takusansu #超算数 全部で68ページもある大作。今は読む余裕がないので、誰か解説して欲しい感じ。
math.oshirase.com/%e9%87%8f%e3%8...
↓
掛け算の順序について I ー 離散量からはじめる ー
math.oshirase.com/wp-content/upl...
↓
掛け算の順序について II ー 連続量の場合 ー
math.oshirase.com/wp-content/upl...
タグ: 超算数
posted at 22:08:10
OokuboTact 大久保中二病中年 @OokuboTact
@takusansu ざあっーと読みました。
銀林浩理論との類似を強く感じたけれど、
小島順氏は「小学生算数の順序固定」に賛成しているのか否定しているか全くわかりませんでした。
タグ:
posted at 22:09:41
@OokuboTact @takusansu #数楽 数学と物理の関連で幾何的に「量の次元」について理解することに興味があるならば、ツイッターで何度も紹介しているようにH.ワイル『空間・時間・物質』がお勧め。例えば、添付画像のような説明がある。今読むと結構分かり易い。
twitter.com/genkuroki/stat...
タグ: 数楽
posted at 22:12:26
OokuboTact 大久保中二病中年 @OokuboTact
@takusansu [ 銀林浩理論と小島順理論は似て非なるものか?】
参考 2017年の議論 pic.twitter.com/RICX1YwNOC
タグ:
posted at 22:15:50
非公開
タグ:
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
@takusansu (続き)
小島順理論を巡る説明(2017年)
#超算数 pic.twitter.com/UZu9M6meen
タグ: 超算数
posted at 22:35:17
文部科学省の「学校の新しい生活様式」に関するマニュアルが改定され、次亜塩素酸水を有人空間に噴霧することについて、「児童生徒等がいる空間で使用しないでください」から「メーカー等が工夫しているため学校側が使用の是非を吟味して判断せよ」に変わりました。
確認ですが、右側が新しい方です。 pic.twitter.com/3l11x7j11G
タグ:
posted at 22:41:08
有効性も安全性も確認されておらず、WHOも「いかなる状況でも推奨されない」としているにも関わらず、「推奨しない→吟味して使ってください」に判断を変えた文部科学省は、「科学」の名前を外した方がいいのではないでしょうか。
www.google.co.jp/amp/s/seniorgu...
タグ:
posted at 22:41:08
繰り返しますが、私は別に次亜塩素酸水の噴霧を毛嫌いしてるわけではありません。
特に個人的な思い入れはナシです。
ただ、人の健康のために使うのなら、しっかりとしたエビデンスを示してからやって欲しいと思うだけ。
「よくわかってないから使っていい」っておかしくないですか??
タグ:
posted at 22:49:04
ちなみに福岡県の飯塚市教委は、市立小中学校の全ての学級と児童クラブに1台ずつ配備できるよう450台を約1600万円で購入したそうです。
「児童生徒がいない放課後や始業前の活用を考えている」そうです。
安倍のマスクと同じく税金が使われている気がしてなりませんね😇
www.google.co.jp/amp/s/mainichi...
タグ:
posted at 22:53:53
OokuboTact 大久保中二病中年 @OokuboTact
@takusansu 参考
『圏論の歩き方』(2015年)という本の第12章の参考文献でも、小島順と銀林浩の本が参考文献にあがっている。
第12章の執筆は西郷甲矢人氏
#超算数 pic.twitter.com/s1BbLYMDk0
タグ: 超算数
posted at 23:07:55
教育「教えてない漢字や解き方は使ってはいけません」
教育「掛け算を逆にしてはいけません」
教育「読書したら必ず感動しなければなりません」
教育「英語を上手に発音してはいけません」
教育「道徳は大人の求める答えを書かなければいけません」
教育「日本のジ ##RT拡散
タグ:
posted at 23:26:00