黒木玄 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年10月27日(木)
新しい記事がZennにアップされました！#Julia言語 #Julia日本語記事
zenn.dev/ohno/articles/...👈
タグ： Julia日本語記事 Julia言語
posted at 22:52:48
ソーシャルメディアと陰謀論的信念の関係を調べたところ、ヤフコメの閲覧頻度が高いことと関連がある一方で、Twitterの利用頻度が高いことは関連が低いということです。
著者は（予想外に）と書いています。
私もTwitterは陰謀論の巣窟と思っていましたが、そう感じるメカニズムも書かれています。
タグ：
posted at 21:19:36
#統計 添付画像は既出の最近の統計学の適切な使用法に関する論文
journals.sagepub.com/doi/10.1177/02...
が引用している図の孫引きです。
統計ソフトの十分な支援の下で、この図のイメージで考えることを楽にできるようになっていれば、検定と区間推定の間の整合性をなくすことは難しくなります。 pic.twitter.com/6Y17v0uBKL
タグ： 統計
posted at 21:06:27
#統計 「有意差決戦主義」時代に整備されて流行した検定法の中には、「違いがない」の型の帰無仮説専用の検定法になっていて、「パラメータθの値はθ₀である」の型の仮説群を扱うようになっていなかったり、扱えるように拡張されていても標準的に教えられていないものが沢山ある。続き twitter.com/genkuroki/stat...
タグ： 統計
posted at 20:59:26
葛屋カツキ＠7/20-23日初個展 原宿 @S9992020
鹿児島空港のツイートがバズったので、鹿児島県のさつま揚げが美味しかったことをお伝えします！
1.「揚立屋」さんトウモロコシのさつま揚げ、超うんまい。
2.指宿の民宿「たかよし」さんの手作りのさつま揚げ最高！
3.揚げたてのサツマアゲを食べるタヌキ（ふわふわ） pic.twitter.com/tSmXes34oq
タグ：
posted at 20:48:03
@muuumin20 @Hyrodium #Julia言語
import Foo: g してしまうと、g(x)=～ でFoo.gにメソッドを追加できてしまうのがちょっと怖いです。
using Foo: gならばそうできなくなり、どのモジュールのgなのかが分かるようにFoo.g(x)=～と書くことになります。
この手のことで失敗しがちな私にとってusingはありがたいです。 pic.twitter.com/O1fvl3kyJd
タグ： Julia言語
posted at 19:32:34
これはまた面倒なことに。
リプライから規格制定時の議事録が辿れるんだけど、
「日本のみその製造業者さんで、この規格を取得できない業者さんっているのでしょうか」
「基本的にはできるという形で規格を詰めてきました。」
という問答があり、大豆不使用の味噌については全く考慮されなかった模様。 twitter.com/iitom3/status/...
タグ：
posted at 18:48:51
ハフポスト日本版 / 会話を生み出す国際 @HuffPostJapan
公立校の教員に残業代を支払わないことを定める法律（給特法）の抜本的な見直しなどを求める署名6万4000筆を自民に提出した現職教員、西村祐二さん（@kimamanigo0815）
「教員が望んでいるのはお金ではなく、業務を減らすことだ」と訴えています。 www.huffingtonpost.jp/entry/story_jp...
タグ：
posted at 18:35:09
#統計 WMW検定に関する多くの解説で、
* 正規分布を仮定しなくてよい。
* t検定との比較で検出力があまり落ちない。
というようなことばかり強調して、
* 「2群の母集団分布が等しい」という強烈に強いがゆえに否定し易い帰無仮説を採用している。
という事実を強調しないから、誤用が止まらない。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 18:18:09
青線はBrunner-Munzelでのシフトθの95%信頼区間。
橙線はWelchのt検定に付随する母平均の差の95%信頼区間。
Welchのt検定は外れ値に弱いので自由度2のt分布が生成する極端な値によって信頼区間の幅が急激に大きくなる場合があります。
BMのシフトθではそうならない。
github.com/genkuroki/publ... pic.twitter.com/tztgAGuUHr
タグ：
posted at 18:05:59
#統計 添付画像は、2つの群の母集団分布がどちらも自由度2のt分布のとき、サイズm,nの標本をランダムに生成して、その各々ごとに信頼区間をプロットしたものです。
青線はBrunner-Munzelでのシフトθの95%信頼区間。
橙線はWelchのt検定に付随する母平均の差の95%信頼区間。
github.com/genkuroki/publ... pic.twitter.com/kG0SR91UsV
タグ： 統計
posted at 18:03:07
脱線続き(別スレッドでWilcoxonの順位和検定=Mann-WhitneyのU検定について詳しく説明した)
↓ twitter.com/genkuroki/stat...
タグ：
posted at 17:58:47
#統計 Brunner-Munzel検定でのシフトパラメータθの信頼区間は、2群の母集団分布が正規分布なら(不等分散可)、Welchのt検定に付随する信頼区間と似た結果を与え(添付画像①)、外れ値が出るケースでも区間の幅が急激に広くなったりしません(添付画像②)。
青線=BM
橙線=Welch t
github.com/genkuroki/publ... pic.twitter.com/r4FO58TyWs
タグ： 統計
posted at 17:53:24
#統計 Brunner-Munzel検定についても、Wilcoxon-Mann-Whitney検定の場合のように
P(X<Y+θ)+P(X=Y+θ)/2=1/2
の意味での効果量θの信頼区間を計算してくれると便利だと思うのですが、そういう既存のRのパッケージの存在を知らない。
私による #Julia言語 での実装例→ github.com/genkuroki/publ...
posted at 17:42:41
#統計 Wilcoxon-Mann-Whitney検定とBrunner-Munzel検定については私のツイログ検索を参照
↓
twilog.org/genkuroki/sear...
タグ： 統計
posted at 17:37:30
#統計 この仕組みを使うと、X=Y+θの意味での効果量θの信頼区間を計算する方法が得られます。
例えばRのwilcox.testをconf. int=Tで使えば、X=Y+θの意味での効果量θの信頼区間も表示してくれます。
github.com/genkuroki/publ... pic.twitter.com/UUNBmgs3HU
タグ： 統計
posted at 17:36:02
#統計 Wilcoxon-Mann-Whitney検定は、帰無仮説として、「2群の母平均が等しい」のような弱い条件ではなく、「2群の母集団分布が等しい」という強烈に強いがゆえに否定し易い条件を採用する代わりに、データの数値の順位情報しか使わない検定。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 17:35:52
@iitom3 @uwajima_iimiso 今回の農文協のJAS規格制定の際の議事録をこちらに貼ります。味噌の定義について海外のオリジナル味噌と差別化するために「大豆と麹菌の使用」を含めているようです。調査会委員もまさか日本に大豆を使わない味噌があるとは認識していなかったと推測されます。
twitter.com/nakaji_minami/...
タグ：
posted at 16:31:22
Michal Sojka (@wenta @wentasah
New release of EmacsVterm.jl is out. It improves #JuliaLang integration with Emacs. New features: Can display images 🖼️ in Emacs buffers. Emacs part is now upstreamed in julia-repl package. github.com/wentasah/Emacs...
タグ： JuliaLang
posted at 16:03:55
葛屋カツキ＠7/20-23日初個展 原宿 @S9992020
【心からおすすめ】鹿児島から飛行機で絶景の島「沖永良部島」に行きました。
「もぉりもぉり」に私のサイン色紙を置いていただいてます😂女性１人でも入れる素敵な居酒屋さんです。（予約推奨）
島のラーメン屋「しん真」の味噌味も絶品でした！
3月頃から海に入れるらしい、また行きたい✨ pic.twitter.com/ok73HkUU6y
タグ：
posted at 15:37:56
藤本 衡 / Kou Fujimoto @huzymotto
経験的には学生実験の所要時間はガンマ分布的な傾向があって、たまーにすごく時間のかかる班がある。そういう班でも現実的な時間で実験が終わるように内容設計すると、結果的にかなり多くの班は予定時間よりかなり早く終わることになるんだよね…
ja.wikipedia.org/wiki/%E3%82%AC...
タグ：
posted at 15:20:10
#統計 脱線続き。既出であることを示すリンク
Wilcoxon-Mann-Whitney検定の使用可能条件の厳しさと、WMW検定の代替案としてのBrunner-Munzel検定の紹介については私のツイログを参照。
↓
twilog.org/genkuroki/sear...
タグ： 統計
posted at 14:50:10
#統計 脱線。Wilcoxonの順位和検定=Mann-WhitneyのU検定は、何の前提も無しに使うと、「２群の母集団分布は等しい」という仮説の検定になります。本当に扱いたい問題が「片方がもう一方より優れているか否か」であるときにWilcoxonの順位和検定を使いたい場合には非常に強い前提が必要になる(既出)。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:41:50
#統計 統計学的ツールをその詳細を完璧に理解してから使うというのは至難の技で、常にそのように要求するのは非現実的だと思います。
統計ソフトによる支援が重要です。
しかし、以下のリンク先で紹介した論文にある統計学の使い方へのソフト側の準拠は非常に遅れています。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:33:11
#統計 あと、以下のツイートで批判対象を広げてしまいましたが、「批判」は「問題を指摘することによって、建設的な提案をすること」を意味し、「攻撃」を意味しないことにも注意してください。(本当はこういうことは言う必要がない。) twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:21:56
#統計 「有意差決戦主義」の影響下で統計分析の報告をすると、たとえ効果の推定も同時に報告していても、相当に味の悪いことになりがち。
そういう障碍を取り除くための概念的にシンプルな枠組みが以下のリンク先で紹介した前者の論文に書いてあります。
理解者が増えると良いと思います。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:15:18
#統計 P値と区間推定の結果に整合性がないように見える結果が得られた場合にも、正直に両方報告する責任が生じています。
自分にとって都合の良い側のみを報告することは、自分自身を科学的でない人間にしてしまうことになります。
正直に両方報告して言い訳を書かなければいけない。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 14:10:35
#統計 最近の論文
journals.sagepub.com/doi/10.1177/02...
stat.columbia.edu/~gelman/resear...
を読むと、最近の統計学の正しい使用法に関する議論は、まさにこのスレッドが批判している対象の問題を扱っているように見える。
前者の論文のP値函数の使用の勧めはこのスレッドで扱っている問題への1つの対策になっています。
タグ： 統計
posted at 14:00:23
#統計 一般論として、以下のリンク先の疑問を解消するためには、効果の指標θの信頼区間が自然に対応するのは
数値θ₀を任意に動かせる場合の仮説θ=θ₀の検定法
であることに注意する必要があります。その場合には、
仮説θ=θ₀が棄却されない数値θ₀の範囲
として信頼区間が自然に決まります。続く twitter.com/tendy10tendy/s...
タグ： 統計
posted at 14:00:16
#統計 以上によって、
www.shionogi.com/content/dam/sh...
の内容に沿った質問への回答を終えます。添付画像はこのPDFより。手書きの青字は私によるコメント。
この具体的な場合についてではなく、一般論としての回答を以下で行います。 pic.twitter.com/aUIIhnmhpi
タグ： 統計
posted at 14:00:15
#統計 リンク先スレッドの疑問に答えます。長いです。
我々が相対的により興味ある問題は、
(W) 症状消失までの時間の真の分布がプラセボ群と等しいか否か
ではなく、薬の効き目に関する
(M) 症状消失までの時間の中央値がプラセボ群より真に短いか否か
の方です。これが件の問題の状況。続く twitter.com/tendy10tendy/s...
タグ： 統計
posted at 14:00:10
葛屋カツキ＠7/20-23日初個展 原宿 @S9992020
思いもよらぬ大サービスに機内で笑ってしまった！
鹿児島空港ーーまた来るよ！ pic.twitter.com/WJCMvm24ex
タグ：
posted at 13:03:04
赤ベコもらったので猫パンチしたら揺れてかわいいかなって思ってやってみたんだけどそこまでやれとは言ってない…😂 pic.twitter.com/Q3Mkw4XFu7
タグ：
posted at 10:52:01
EMの比嘉照夫氏が叙勲されたとのことで知ったんですけど。紺綬褒章は500万円～。
比嘉氏は「瑞宝章 教育研究功労」です。
www8.cao.go.jp/shokun/hatsure...
＞EM菌開発者は叙勲されるし、文部科学省のホームページにはEM菌教育のページまであります
twitter.com/breathingpower...
タグ：
posted at 10:05:54
いや。素の pharo.org でも動かせるように and/or ふつうに素直に こうか…
| a |
a := (2 to: 100) asArray.
a difference: (a gather: [:x | x * a])
でも、Smalltalkで書くと、同じ方針でも途端にメモリがもったいなく感じて、もうちょっとどうにかしたくなってしまう不思議^^; pic.twitter.com/Lhu3Wi8OOH
タグ：
posted at 09:38:40
@GreatDemon1701 @tendy10tendy @golgo_sardine @sekibunnteisuu #超算数 順序を変えただけで【あれもこれも正しければ公式ですと言ってしまうのは、ただ混乱を招くだけ】とは、あたかも公式が公式球のようにそれを使わないと正当なゲームが成立しないものであるかのよう。また、正方形は一辺×一辺じゃないと公式ではないという見解もtwitter.com/genkuroki/stat...。
タグ： 超算数
posted at 08:22:43
面白かったので、Squeak Smalltalk でも似た方針で^^;
| a |
a := (2 to: 100) asArray.
a difference: a +* (Matrix rows: 1 columns: a size contents: a)
"=> #(2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97) "
タグ：
posted at 07:32:46
APLで書くとこんな感じ？
{(~⍵∊⍵∘.×⍵)/⍵}1↓⍳100
tryapl.org/?clear&q=%7B(~... twitter.com/genkuroki/stat...
タグ：
posted at 07:18:19
ソ連が自由でないことは誰の目にも明らかで、それが科学におよぼす結果もまたはっきりしています。最高の例の一つは、植物が後天的に得た特性を次世代に伝えうるという、遺伝の理論を唱えたルイセンコでしょう。
タグ：
posted at 00:02:19
