黒木玄 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 を見て下さい。
2022年08月05日(金)
#統計 Brunner-Munzel検定を使った2つの軍が互角になるようなシフトaの信頼区間は、平均の差の信頼区間のように使えてかつ外れ値に対する頑健性もあるので、もしかしたら便利に使える可能性があると思いました。
こんな感じのことをやってみました。
↓
nbviewer.org/github/genkuro...
タグ: 統計
posted at 23:59:18
#統計 #Julia言語
正規分布の標本については、
Brunner-Munzel検定を使った互角になるような分布の位置のシフトaの信頼区間
と
Welchのt検定に付随する平均の差の信頼区間
はほぼ一致します。
しかし、外れ値が出る分布では後者は大きく広がる傾向がある(当たり前)。
nbviewer.org/github/genkuro... pic.twitter.com/4XHgyHhfly
posted at 23:56:01
#統計 #Julia言語
添付画像中のグラフは、X と Y + a が互角になるようなaのP値函数です。
橙の破線がaの点推定で、aの95%信頼区間が confint_location として表示されている。
このaに関する信頼区間を、Welchのt検定での平均の差の信頼区間と比較してみました。続く
nbviewer.org/github/genkuro... pic.twitter.com/Bp218GjjjR
posted at 23:50:06
#統計 #Julia言語 BM検定は、X軍とY軍から一人ずつ無作為に選んで戦わせたときのY側の勝率pに関する検定です(戦闘力が高い方が必ず勝つ設定)。
だからpに関する信頼区間を構成できる。
それとは別に、X軍とY群を互角にするような戦闘力のシフトの信頼区間も定義可能。続く
nbviewer.org/github/genkuro... pic.twitter.com/FW0n0lQUwB
posted at 23:45:57
@asanebou13 この20局の結果から最尤法により推定した里見先生のレーティングはなんと1730(レーティング差で勝率が計算できますが、20局の結果の確率が最大となるRが1730に)。全棋士中25位に相当します。実際、レーティング上位澤田七段、出口六段に勝っていますし・・・。
タグ:
posted at 23:39:06
#統計 #Julia言語
github.com/toshi-ara/brun... で指摘されているように、set.seed(1290)後のデータ生成で、100個中4個の値が私の実装と一致しませんでした。
これは、Rのbrunnermunzelパッケージは ⪅ の実装に失敗しているのではないかと思いました。
nbviewer.org/github/genkuro... pic.twitter.com/EK4PauEkNd
posted at 23:38:02
#統計 #Julia言語 Rのbrunnermunzelパッケージのgithubリポジトリの方も見てみた。
すると、permutation版の実装にバグがある疑いがissueになっていて、ずっと放置されていました。
これも確認してみました。続く
github.com/toshi-ara/brun...
posted at 23:34:19
#統計 #Julia言語 Rのbrunnermunzelパッケージの
cran.r-project.org/package=brunne...
にアクセスすると
Package ‘brunnermunzel’ was removed from the CRAN repository.
となっていました! pic.twitter.com/pu5ROsFRfy
posted at 23:31:38
#統計 #Julia言語
okumuralab.org/~okumura/stat/... にある例と同じ数値を私の実装がアウトプットしてくれることの確認。
github.com/genkuroki/publ... pic.twitter.com/1qD8TF8Ayt
posted at 23:29:12
#統計 #Julia言語 Brunner-Munzel検定を私も試しに実装してみました。色々やっています。
* Y側の勝率の信頼区間
* Yと互角になるようなXのシフトの信頼区間
* permutation版
* Welchのt検定との比較
* 信頼区間の比較を含む
github.com/genkuroki/publ... twitter.com/genkuroki/stat...
posted at 23:25:10
アイスランド気象局によると8月3日、首都レイキャビク近郊のファグラダルスフィヤル山が噴火した。この山は昨年も、6カ月にわたり噴火していた。 pic.twitter.com/5ArKHev6HO
タグ:
posted at 23:00:00
ビタミンK欠乏症に由来する頭蓋内出血が減ったら、こんどはビタミンKシロップの意義を見失う人が増えたというあたりがねえ……なんともはや。
昔の、かわいい1~2か月の赤ちゃんを失う時代に戻りたいのかね
タグ:
posted at 21:40:54
@physics303 行列積が使えるのではないかと思いました。コードしか見ていないので、使う際にはよく確認してください。 #Julia言語
nbviewer.org/gist/Lirimy/26...
タグ: Julia言語
posted at 20:27:25
里見香奈女流の公式戦成績
2021/7/23/から
〇黒田✖村田顕〇山本真〇高田〇本田
〇西川✖大地✖星野✖狩山〇浦野
〇澤田〇池永〇冨田〇古森〇浦野
〇出口✖今泉〇柵木〇折田✖星野
直近20局は14-6(.700)
タグ:
posted at 19:27:44
@Peraia2 #統計 母分散はσ² = E[(Xᵢ - μ)²] (E[ ]は期待値を取る操作)なのですが、動画では母分散とは異なるものを母分散と呼んでいます。 pic.twitter.com/4udfgDL0j4
タグ: 統計
posted at 18:35:41
@LeonSatyros 実際に教科書に証明がしっかり書いてあっても、「教科書に書いてある証明が正しいと先生が言っていたから正しい」のように考えちゃダメで、自分自身で証明を再構成して細かいギャップを埋めて、理解に必要な例を自分で作ったりしないと、許してもらえないです。
タグ:
posted at 17:06:07
@banyanyanmi YouTubeの動画の(不当な方法による)バックアップの仕方へのリンクが少し前の ocw.kyoto-u.ac.jp/news/6/ にあったので、バックアップが必要になると理解していたが、「チャンネルは残るようですね」とのこと。
不安だ。 twitter.com/tiseda/status/...
タグ:
posted at 16:58:45
@banyanyanmi いやあ、これにはびっくりですよね。
ocw.kyoto-u.ac.jp/news/6/ を今見たら、動画のバックアップの仕方は後でアナウンスするとなっていました。
すでにYouTubeにあるものはそのまま放置してくれないのかな?
m.youtube.com/watch?v=vz9cZn... を運良く見付けて繰り返し宣伝して来たのに! twitter.com/tiseda/status/...
タグ:
posted at 16:55:35
これ,深さ2の2つのテンソルGとT間の Multilinear Product っていう積の計算なんですけど,こんなに for loop 重ねるしかないのでしょうか.なんかすごい遅くなるようなことをやっている気がして不安.
定義はこの論文のdef5.を参照した.ieeexplore.ieee.org/abstract/docum... pic.twitter.com/IqrlkkGYyZ
タグ:
posted at 16:06:52
@samakikaku @breathingpower ヤマダHDもヒノキヤも知らんぷりかしら
左巻先生の記事
>日本ホメオパシーは福島県内で採取した土を希釈したレメディをペットボトルに入れ河川に投入
放射能を消すとするパフォーマンス
創生水・EM菌・水素水について放射能除去作用検証
結果は水道水除染と何の差も無い
president.jp/articles/-/497...
タグ:
posted at 15:15:43
@samakikaku @breathingpower 2010年8月25日 朝日新聞
ホメオパシーは「荒唐無稽」学術会議が全面否定談話
archive.ph/tS5ZH
>比嘉照夫氏、DND出口俊一氏はホメオパシー支持者。
twitter.com/samakikaku/sta...
タグ:
posted at 15:05:24
@samakikaku @breathingpower DND出口俊一氏はいつからEM研究機構顧問なんでしょうか。
>DND出口俊一氏(EM研究機構顧問 ヒノキヤグループ社外取締役)の「日本学術会議のホメオパシー批判とEM菌潰しの構図は同じ」が当時公開のFacebookにあった。
2010年8月25日 朝日新聞
twitter.com/samakikaku/sta...
タグ:
posted at 14:56:22
@samakikaku @breathingpower DND出口俊一氏は自分で役員やめたんですよね🤔
役員報酬
2021年12月期の社外役員4名に対する報酬総額は174万円増額で3115万円
irbank.net/E00325/fee
irbank.net/E00325
EMの記載はなし
↓
2014.1.1~2014.12.31
株式会社桧家ホールディングス
有価証券報告書
www.hinokiya-group.jp/ir/wp-content/... pic.twitter.com/VQJHTIWLQn
タグ:
posted at 14:46:07
「Juliaで精度保証付き数値計算 」更新しました。
Clenshawのアルゴリズム
taklab-blog.blogspot.com/2021/01/rigoro...
チェビシェフ多項式を用いた補間(チェビシェフ補間)の関数値を計算する方法O(n)で計算できるので最速を誇るが、区間演算と組み合わせると注意が必要です。反復計算に区間演算は弱い…
タグ:
posted at 14:21:17
@sekibunnteisuu うちの子が小1のときの学校のテストです。4+1+4はバツらしいです笑 まぁ、これで良いよとフォローしたし、私は積分定数さんと違うスタンスなので、教員個人に対してはなにも求めないですが…(単純にスタンスの違いです) pic.twitter.com/esJhRBdkI3
タグ:
posted at 13:49:59
OCWを拡充し続けるMIT,
理工系無料講座を拡充し続ける中国大学MOOC(icource163)
理工講義中心に受講すべて無料、各所テストセンターで認定試験受ければIIT他の正式な単位として認めるインドswayam
一方
予算が出されないため消滅する日本。
一体どこで差がついたのか…増税、予算削減、財務省の存在
タグ:
posted at 12:27:12
#超算数
いきなりエレファントな方策を教え込むと、全くナンセンスな苦肉の策をやってしまうことになると思う。
わり算の筆算だと、「たてる」「かける」「ひく」「おろす」という手順を徹底させることが行われているらしい。
www.shinko-keirin.co.jp/keirinkan/sans...
タグ: 超算数
posted at 11:06:03
非公開
タグ:
posted at xx:xx:xx
9価の定期接種化も進みそうです。
「委員からは、海外の接種状況を参考に、対象を9~14歳へ引き下げ、2回接種とするべきだとの意見が出た。」
↑ここも重要。15歳になるまでに打ち始めたら3回も打たなくて良いです。
9価ワクチン「問題なし」 HPV定期接種化へ議論
www.sankei.com/article/202208...
タグ:
posted at 10:52:12
Katsusuke SHIGETA・重田 @shigejam
@horilab OCWとMOOCについては、私が代表幹事のOpen Education Japan(oejapan.org)を皮切りに、意見を表面する準備を進めております。
タグ:
posted at 10:31:06
Katsusuke SHIGETA・重田 @shigejam
京大はFDの先端的拠点であり、京大OCW・MOOCは学習者・教育者・研究者にとって大変高い価値のある教育資源で、取組が終了されることの損失は計り知れません。 twitter.com/horilab/status...
タグ:
posted at 10:30:44
#統計 その動画解説の作成者は母分散の概念を誤解しているので、何も分かっていません。
YouTubeで学ぼうとする場合にはデタラメに感心して受け入れてしまうリスクに注意が必要だと思いました。 twitter.com/peraia2/status... pic.twitter.com/3875vovsCx
タグ: 統計
posted at 09:56:17
@setuko1234 @takedayasuhiro @hhhira @mcranagalle なるほど。誰でも安易に専門外に口を出すべきではないね。
紀藤弁護士は統一協会の件では信頼絶大なので、そちらの追求に力を入れ続けてほしいね。
タグ:
posted at 08:36:29
ICYMI, Geometry2D.jl and BeltTransmission.jl have had their 0.1 release:
mechanomy.com/products/BeltT...
mechanomy.com/products/Geome...
#JuliaLang pic.twitter.com/7WRJ0225Bh
タグ: JuliaLang
posted at 06:49:28
@hhhira @mcranagalle 紀藤正樹弁護士は、フクシマの放射線で死ぬ死ぬ詐欺師の木下黄太、三田茂、西尾正道、井戸謙一などと組んでますね。 pic.twitter.com/9Qo8R4DAn4
タグ:
posted at 03:17:35
@genkuroki @sunchanuiguru @keiosfcimailab #超算数 今井むつみ氏がどう考えているかはともかく、確実に悪影響はある。 pic.twitter.com/pFbrCpR8oG
タグ: 超算数
posted at 01:50:15
@genkuroki @sunchanuiguru @keiosfcimailab #超算数
>子供はなるべく与えられた数字だけを使って答えを導く傾向にある、
そうしろと教わる。あるいは直接教わらなくても、そうしないとならないと思い込んでいる。
本当に子どもが元から持っている性質かどうかは疑わしい。 pic.twitter.com/V0tvoIn46X
タグ: 超算数
posted at 01:46:37
@genkuroki @sunchanuiguru @keiosfcimailab #超算数 妙な納得した人と、真っ当な突っ込みしている人 pic.twitter.com/yIKwLkkNqp
タグ: 超算数
posted at 01:40:46
@genkuroki @sunchanuiguru #超算数 今井むつみ氏@keiosfcimailab は、おかしな算数教育が行われていることを知っているのかな?
「子供はこうしがち」と言っているけど、教え方の影響は検証したのかな?
タグ: 超算数
posted at 01:29:21
@genkuroki @sunchanuiguru #超算数
子供にとって交換法則が直感的なものじゃないというが、
かけ算順序をしつこく教わればそうなっちゃうよね。 pic.twitter.com/LdUULvpcb5
タグ: 超算数
posted at 01:26:02
@sekibunnteisuu @sunchanuiguru これ、科学的根拠があるんですかね?
多くの一般人は肩書きや聞き慣れない用語に権威を感じて、一理あるのかなと思ってしまう。
きちんとコントロールを用意して比較したもの以外は相当に要注意。
タグ:
posted at 01:17:05
#超算数 5割増しだと正解することが多いらしい。「大体の値を直感で答える」なら、以外と正解に近い値を出すかも。「出てきた数値で何らかの計算処理をしないとならない」が邪魔になっているように思う。
twitter.com/sekibunnteisuu...
タグ: 超算数
posted at 01:13:18
#統計
私も渡辺澄夫さんの仕事にはb函数方面から興味を持ちました。
今思えば、その偶然は非常にラッキーで、渡辺澄夫さんが批判している「主義に基く統計学」の洗礼を完全に回避することができました。 twitter.com/_zak3/status/1...
タグ: 統計
posted at 01:07:38
@sunchanuiguru 全然知らない人だった。今井むつみ氏は最近よく見かけて、今読んでいる「誰が国語力を殺すのか」でも出てきます。
記念館は知らなった。メットとか飾っているのかな?w
タグ:
posted at 01:04:13
#統計 ここでやっている計算は数学的にはq二項係数
q-binom(200, 100)
のqの多項式としての計数をすべて正確に求めているのと同じことです。
github.com/genkuroki/publ...
タグ: 統計
posted at 01:02:03
#統計 素朴にすべての組み合わせを生成する方法でやると、200個から100個選ぶ組み合わせの個数だけループを回すことになって死ぬ。
₂₀₀C₁₀₀ ≈ 9×10⁵⁸ 😱
github.com/genkuroki/publ... pic.twitter.com/prlslVSHKM
タグ: 統計
posted at 01:00:25
@sekibunnteisuu 堀元見ですか…
胡散臭いやつです。
算数とは関係ないですが、堀元見が経営していた(いまは放棄)キャンプ場から至近距離に藤本敏夫(加藤登紀子の夫)記念館があります。
タグ:
posted at 00:58:10
#統計 添付画像
①場合の数をInt64で計算したら、こんなことに!
②Int128でもダメ!
③BigIntでやっと中心極限定理による正規分布近似を再現できた。
④場合の数をFloat64で計算してもほぼ同じグラフをプロットできる。
github.com/genkuroki/publ...
Exact Wilcoxon-Mann-Whitney without ties pic.twitter.com/oX7KTKHUVW
タグ: 統計
posted at 00:55:50
#統計
github.com/genkuroki/publ...
順位に重複がない場合のWilcoxon-Mann-Whitney検定で使う確率分布(正確には場合の数)の正確な計算
3通りに実装。最後のやつは事前割り当てを使えば新規メモリ割り当てゼロで動く。
まじめに効率的な実装はしていない。分かり易さを重視。 pic.twitter.com/PohfemJlnB
タグ: 統計
posted at 00:52:53
Metal.jl - A GPU backend for Apple hardware: This talk presents updates on the development of a #GPU backend for Apple hardware (specifically the M-series chipset) along with a brief showcase of current capabilities and interface. youtu.be/IARikXzRU7s
#julialang #juliacon
posted at 00:35:42
@phasetrbot ありがとうございます。
こちらにとりあえず3日まで動作確認したコードが置いてあります。残りは順次確認次第載せる予定です。
cometscome.github.io/YukiNagai/ja/b...
タグ:
posted at 00:03:40