黒木玄 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年01月24日(月)
#統計 推定量に不偏性を課すとそのせいで平均二乗誤差が大きくなってしまうことがよくあります。
特に古典的に有名な話はStein推定の話。
nbviewer.org/github/genkuro...
より初等的には、正規母集団の不偏分散は平均二乗誤差が大きくなる簡単な例になっている。
不偏性にはこだわるのはよくないっぽい。
タグ： 統計
posted at 23:49:00
#統計 #Julia言語 iidの不偏分散と補正無しの標本分散ではそれらの推定先である母分散σ²に対する平均二乗誤差が小さいか？
母集団分布の尖度(正規分布で0になるやつ)が-3/2より大きいならば十分大きなnについて、
不偏分散の平均二乗誤差 ＞ 標本分散の平均二乗誤差
github.com/genkuroki/publ... pic.twitter.com/IjGZaZsJFu
posted at 23:35:32
「共通テスト科目にしなければ永久に改善しません」の一言には度肝を抜かれたが、この大人の都合丸出しな本音に「高大接続」「大学入試改革」の問題点がほぼ全て詰まっていると思う。
タグ：
posted at 22:49:59
ご紹介ありがとうございます、計算が複雑だったことで詳しい考察ができていませんでしたので、参考にします。 twitter.com/genkuroki/stat...
タグ：
posted at 21:40:08
マックイーン！知ってっか〜、天才ラマヌジャンを見出したことで有名なハーディもめちゃくちゃ業績がある数学者なんだぜ〜。ハーディは純粋数学は「役に立たないからこそ美しい」という美学の持ち主だったんだけど、後年彼の数論がコンピュータに、調和解析が制御に応用されていくのがおもしれーよな！
タグ：
posted at 21:32:29
新型コロナウイルスはオミクロンで『ただの風邪』になったのでは？
という意見も多いが、感染力(実効再生産数)もインフルエンザより高く、重症化もしやすい。
インフルエンザでも、小児で脳症やけいれん、心筋炎をきたす。新型コロナウイルスはまだインフルエンザまで管理可能な状態になってはいない
タグ：
posted at 21:30:31
なぜか宮沢孝幸や長尾和宏など、妊婦も脳出血高齢者も『診療しない』人達が『5類で解決』って言うんだけど、何類だろうが、感染者数は減らない
なおかつ5類にすればリアルタイムな患者数は反映されない。つまりピークアウトしたかどうか、現場が把握しづらくなる twitter.com/MIKITO_777/sta...
タグ：
posted at 21:24:54
あっという間に流れていってしまったが、ベーシックインカムは勤労意欲を削がないという実験結果を紹介する記事を見た。好感する論説だったが、BI論者にとってはまずい結果なんじゃない。企業が人件費を上げるインセンティブにならないということなのだから。
タグ：
posted at 21:24:41
タモリ倶楽部オープニングテーマ
ロイヤルティーンズのショートショーツ♬︎1958年
@retoro_mode pic.twitter.com/DPrN9x1xty
タグ：
posted at 21:24:27
入院中は、ネット環境がよろしくなかったので画質を落としてYouTubeで柳時熏さんの10秒碁をよく見ていました。
www.youtube.com/channel/UCB_68...
タグ：
posted at 21:13:39
入院しているときに気を付けたことは、自分の病気についてインターネットで検索しないこと。
ついさっき検索してみて、「うん。入院中に検索しなくて良かった」と思いました。
検索で見つかる内視鏡検査の写真を見ると、「あれ？俺の方がずっとひどかったんじゃね？」と感じられるものが多数。😱
タグ：
posted at 21:04:35
#統計 一様最強力不偏検定は不偏検定という特殊な検定の中でのみ一様最強力であるにすぎないので、不偏でない検定との比較では、パラメータ(横軸)と真の値(縦軸)の組み合わせによっては、強さで負ける場合も出て来ます。以下の左下のグラフの赤い部分では実際に負けている。 pic.twitter.com/mTgO2Fu7ht
タグ： 統計
posted at 20:51:18
#統計 しかし #Julia言語 の強みを活かしたお気楽プログラミングによって、Fisher検定(のオッズ比=1という条件を外した一般化)の一様最強力不偏版にも対応しているコードになっているので、興味がある人は自分で試せます。
nbviewer.org/github/genkuro... pic.twitter.com/8d14DgRYcS
posted at 20:42:06
#統計 以上のパワーのプロット比較はn=10の場合。
以下のグラフはn=100の場合です。
n=100にすると、一様最強力不偏二項検定と正確二項検定(Rのbinom.test)の違いはかなり小さくなります。
上段の 2つのパワーのグラフがそっくりになっている。 pic.twitter.com/czX96klltt
タグ： 統計
posted at 20:35:09
#統計 "exact" はRのbinom.testのデフォルトと同じ正確二項検定で、eqtailedは片側確率の2倍でP値を定義した場合の二項検定です。
左下のプロットが青いので、前者の方が強力です。
なぜか片側確率の2倍が好きな人達がいる不思議。
そんなにパワーを下げたいのか？ pic.twitter.com/S5BSwCJME6
タグ： 統計
posted at 20:29:04
#統計 パワーの比較もしています。
二項分布モデルの場合の一様最強力不偏(umpu)検定と、不偏検定ではない正確二項検定(exact)の比較。
左下のグラフの青い部分ではumpuが強く、赤い部分ではexactが強い。
横軸は帰無仮説の"p"で、縦軸は真の分布の"p"です。
nbviewer.org/github/genkuro... pic.twitter.com/Z7k7oJqBiX
タグ： 統計
posted at 20:23:40
実用性と無関係な統計学の話題は気楽にできて良い。
少々間違っていても「仕事で使って被害を受けた」の様なことは言われずに済む。
緊張するのは実際によく使われている道具の話をするとき。
タグ：
posted at 20:16:55
#Julia言語 umpu!(r, dist, α)は指数型分布族の1次元の有限離散分布distと有意水準αと配列rを受け取って、配列rの中に一様最強力不偏検定の棄却確率を格納する函数。
umpu(dist, α)はそれから作られる検定函数phiを返す函数。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 20:06:15
#統計 最初はさらに素朴に、モデルの確率函数のパラメータによる偏微分を使う方法を自動微分を使った実装も試してみました。
その辺の試行錯誤の様子も残しておいた方がいいように思われたので
github.com/genkuroki/publ...
で公開しておきました。
自動微分は普通に使える便利な道具です。
タグ： 統計
posted at 20:00:32
#統計 #Julia言語
添付画像が、指数型分布族の1次元有限離散分布の場合の、素朴な一様最強力不偏検定の実装です。
どこにも二項分布のような具体的な分布を指定するコードが入っていない点に注目。
素朴な実装ではありますが、その点においてはものすごく抽象化されたコードになっています。 pic.twitter.com/jXppjJ8mPY
posted at 19:55:49
@arts_lib #統計 なお、両側をα/2未満に分ける方法は、n=4, p=1/3, α=0.1の場合にはうまく行きますが、必ずしもうまく行くとは限らないので注意が必要です。
例えばn=4, p=0.2, α=0.1の場合にはうまく行きません。 pic.twitter.com/gKVIVfodju
タグ： 統計
posted at 19:48:40
@arts_lib #統計
#Julia言語 で指数型分布族の一次元有限離散分布モデルの両側検定一般で使える一様最強力不偏検定を実装しました。パラメータに関する微分を使わずに実装してあります。
手計算をする場合にもそうした方が簡単です。
twitter.com/genkuroki/stat...
posted at 19:35:40
#統計 一応、Rのumpパッケージのumpu.binom 函数の結果と私の実装による計算結果を数値的に比較して一致していることを確認しています。
www.rdocumentation.org/packages/ump/v...
私の実装はそのパッケージのソースコードを見ずに独立に書いたものなので、これで双方の信頼性が少し上がっていると思います。
タグ： 統計
posted at 19:25:03
#統計 #Julia言語 二項分布モデルやFisherの非心超幾何分布で使える一様最強力不偏(umpu)検定を実装してみました。
↓
nbviewer.org/github/genkuro...
何の工夫もない素朴な実装です。
まだ細かい誤りが残っている可能性があります。見つけたら教えて下さい。 pic.twitter.com/MguCCKC9JU
posted at 19:19:07
ジギット野崎 Taiga Nozaki @glocker_tanz_jp
北西コーカサスの先住民、アディゲ諸族に伝わる剣舞を披露させていただきました🦅
#Адыгэ #ヨールカ祭 pic.twitter.com/Nav6Eug8EJ
posted at 18:54:28
小学校の先生はなぜこうなるのか質問されてもなんとか対応してほしい。かけ算の順番とかこだわっていないで。
キーワードは循環小数です。 twitter.com/Harar1_/status...
タグ：
posted at 17:50:34
#Julia言語
自分でパッケージを開発していてLTSもサポートしている人以外は、LTSの存在は忘れた方が良くて、最新安定版を使うというので良いと思います。
新しいv1.*のJuliaでも過去のv1.*のコードは動くし、機能的にも計算速度的にも改善されているので、最新安定版への移行のメリットは大きい。 twitter.com/h_hiro_/status...
タグ： Julia言語
posted at 17:30:02
@MKS403 @MMBIZ4 twitter.com/Nijntje_Nijntj...
失礼します。これは冗談じゃなく本気でおっしゃっているのでしょうか？
この答案を書いた子は、縦と横が分かっていないのでしょうか？ pic.twitter.com/IAGha1GFdJ
タグ：
posted at 16:50:12
Tomokazu KASHIO （加塩朋 @Tomokazu_Kashio
重積分の積の部分をなぜかマイナスで計算しているレポートが多くて
∫∫f(x)g(y)dxdy=∫f(x)dx - ∫g(y)dy
なんでだろ、授業中変なこと教えちゃったかな、って思ってたんだけど。
最初は積を表す「・」だったのが、丸写しを繰り返すうちに横長になり、いつの間にか「ー」になったことに気が付いた👹
タグ：
posted at 16:17:36
これとか: github.com/JuliaLang/juli...
これとか: github.com/JuliaLang/juli...
これ: github.com/JuliaLang/juli...
タグ：
posted at 16:12:37
こういうのが、まさにポルポト的なんだと思うよ。
「政治的に正しい胸の大きさ」があるかのような主張はリベラルともフェミニズムとも反し、強いていえばファシズム twitter.com/tawarayasotats...
タグ：
posted at 15:40:32
トリビアの泉で衝撃的だったやつ選手権
入選
・「123456789」はどう入れ替えても3で割り切れる。全ての組み合わせが割り切れる。
・4/4と6/6と8/8と10/10と12/12は毎年すべて同じ曜日になる。永遠に。
※数学得意な人、本当ですか？
タグ：
posted at 15:37:44
新しい記事がQiitaにアップされました！#Julia言語 #Julia日本語記事
qiita.com/ekzemplaro/ite...👈
タグ： Julia日本語記事 Julia言語
posted at 15:03:35
今朝もいい感じに出た。1回目は中途半端で2回目に気持ちよく出た。ここ数日の具合は非常によい。
退院してからしばらくのあいだ、大腸が回復し切っていない感じが強くあって、正直「参ったな」「便秘怖すぎ」と思っていたのですが、ここ数日間は非常によい。
タグ：
posted at 14:42:32
私は社会学も重要な研究分野だと若いときからずっと思っているので、最近の私のTLに流れて来る否定対象が大き過ぎる非難には与したくない。
気にいらないやつがいるなら名指しで個人を直接的に批判しろ、と言いたい。周囲を巻き込むな。
最近の闇のapjさんの発言には20年以上前から共感している。
タグ：
posted at 12:52:34
非公開
タグ：
posted at xx:xx:xx
両論併記は両論が同程度の確からしさである時にするものであって、ひとつの論が一方的に正しいなら両論併記の必要はないし、反論がマイナーなものなら扱いは小さくするべきです。
地球球体説と地球平面説を同列に扱うような両論併記は有害
タグ：
posted at 11:36:10
(そもそも反原発に命がけのヤツなんかいるかね?安全地帯から騒いでいるだけの印象なんだけど。まあ学者関係なら専門家としての信頼性かけてる人はいるだろうけど)
タグ：
posted at 11:30:35
つーか、いち福島県民として言わせてもらえばこの人の方がよほど迷惑なんだけどな。主張もさることながら、勝手に福島県民の意向を決めつけちゃう所とか。
twitter.com/tswaxAatkyoSSS...
タグ：
posted at 11:27:17
きくまこ先生は命がけじゃないかもしれないけど、こういうデマ信者からの嫌がらせに対峙しているのは結構な労力だろうな、って思うよ。学者としての業績にもならないのにね。 twitter.com/tswaxAatkyoSSS...
タグ：
posted at 11:24:10
職場の方に「スイートポテト作るならスーパーカップがあればいい。芋以外の必要なもの全部あれに入ってる」と教えていただいて無事美味しいスイートポテトが完成しました。生クリームもバターも砂糖もいらんかったんや…
タグ：
posted at 11:18:07
非公開
タグ：
posted at xx:xx:xx
@genkuroki 特に猗窩座はかっこよかったのでみんな真似してますよね。「破壊殺！」たぶん次に出てくる、逃げ回って責任転嫁するのが得意技の鬼は、あまりみんな真似しないと思う(子供のうちは)
タグ：
posted at 10:01:04
自分ちの子は鬼派の幼児だというツイートもたまに流れて来る。
さすがに子供は鬼滅のどこが面白いかを鋭く見抜いているようだ。 twitter.com/hiyori13/statu...
タグ：
posted at 09:54:14
明らかな誹謗中傷、名誉毀損なら、めでたく（？）当事者の地位を獲得したので、訴えて、削除させる、賠償金をとる、削除するまで賠償金上乗せなど、いろいろなパターンで求めていけばよい。そこで出た判断は、下級審の裁判例であっても積み重ねることに意味がある。相場が決まってくるので。 twitter.com/apj/status/148...
タグ：
posted at 09:49:13
p値に関するASAの提言後、影響力の強い雑誌において、統計学的有意差のない結果がどう結論で扱われていたかを見た研究。「効果がない」や「統計学的に有意差がない」との記載が多く、95%信頼区間と臨床的に意味のある差について十分な注意が払われていませんでした。
www.jclinepi.com/article/S0895-...
タグ：
posted at 09:45:34
ある表現を気に入らないと思うが直接訴訟ができない場合は、そもそも紛争の当事者になれない関係なのだから、他のルートで（雇用等の）法的関係に影響を及ぼして実現しようというのは、権利が無いことをズルして通そうとしているのに他ならない。これをむやみに認めてはならない。 twitter.com/apj/status/148...
タグ：
posted at 09:44:27
#数楽 高木貞治『解析概論』関連
昔から理学部数学科では数学の勉強における具体例の計算の重要性が強調されて来たと思う。
新しく見える抽象的一般論の評価でも非自明で面白い新しい具体例を含むかどうかが常に問われる。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 09:41:14
これは逆に味噌も糞も一緒に考えないとならない場面なんだよね。正義の女神が目隠しをしているのは味噌と糞を区別しないため。 twitter.com/apj/status/148...
タグ：
posted at 09:40:21
言論を巡って解雇等が起きた場合、そっちも訴訟すればというのはただの無責任。その解雇等に関する訴訟は、最初から直接本人を訴えていれば起きなかった可能性が高い上、表現をけしからんと主張する側が手間を省いて誰かにやらせるというゲスい行動原理によって生じたことになる。 twitter.com/apj/status/148...
タグ：
posted at 09:37:02
大学生の生活保護が禁止されているため、虐待から逃れるなどで生活困窮したら手詰まりになってしまう問題。
横須賀市が独自に生活扶助相当と通学費の支給制度化した。
市独自で現金支援策を講じる横須賀市の先進性とともに困窮学生を助けようとしない国の後進性も明らか www3.nhk.or.jp/news/html/2022...
タグ：
posted at 09:34:10
訴訟であれば対等に争って線引きができる。雇用関係のある間柄で規制がかかると、力関係の差から対等に争うのは至難の業。しかも素人が介在するから、判断の基準がブレがちで不安定になる。 twitter.com/apj/status/148...
タグ：
posted at 09:34:03
例のオープンレターが、一応は女性差別的な文化を無くすことを目指しているのと同様、私は、個人の言論について職場にクレーム入れて当然という文化を無くす（本人の提訴一択）ことを目指している。誹謗中傷を容認するという意味ではなく、法によらない表現の制限をするなという意味。
タグ：
posted at 09:30:42
@h_okumura 潜伏期間の使い道はRight tailで95や99の分位点が多そうです（暴露者行動制限の日数根拠）。それで95や99の分位点信頼区間をBootstrapかベイズで出してくるためにRight tailのNon zero probabilityを計算可能なtwo parameter分布を出してくる、ということですね。欲を言えば対数正規やワイブルも見たい
タグ：
posted at 09:28:30
普段愛用しているVScodeのwav-previewというプラグインがパワーアップして波形やスペクトルまで表示されるようになった！再生前に波形が確認できるのは地味にうれしい marketplace.visualstudio.com/items?itemName... pic.twitter.com/sQK0HfMgZI
タグ：
posted at 09:27:34
「Pythonによる制御工学入門」のPythonプログラム（一部除く）をJulia言語で書き直しました．JuliaでもControlSystems.jlを使えば，ひととおり勉強できそうです
サポートページにPython, Matlab, Juliaのコードを掲載しています
y373.sakura.ne.jp/minami/pyctrl
タグ：
posted at 09:04:50
『県民などに対する臨時特別協力要請の案として、２回目までのワクチン接種を終えていない人への不要不急の外出自粛』
ただ、感染予防対策を徹底しないノーマスク・反ワクチンこそ、『自粛』しない人達なんで、効果があるとは思えない
www3.nhk.or.jp/lnews/kofu/202...
タグ：
posted at 08:43:04
『とっとと詰めよ』
『医療者なんて名前の人殺しどもだけが快適に暮らす社会が成り立つわけ無いだろが』
『さっさと崩壊しろ
archive.ph/pizVp』
さすがにひどいので、通報＆発信者情報開示 pic.twitter.com/iPvJZ06o3J
タグ：
posted at 08:27:34
でも鬼はそれを口実にして他人を喰うことが正当化されると思ってて、それを理由に自分の恨みを発散させていいと思ってて、かわいそうだけどそれはやっぱりいかんよ、という大きなテーマが出てくる。なんか、現実世界で展開されているある種のドラマをまさにそのまま戯画化したようで、すごいわ。
タグ：
posted at 08:25:49
本音はむしろ「いいなあ、おめえら、いい思いしてんなあ、ねたましいなあ」というルサンチマンの塊で、こいつらが交互に手を組んでくるので始末に負えず、倒すには両方の首を同時に斬らねばならない！鬼は鬼になる事情があり、それは遊郭や社会の仕組み批判としても機能している
タグ：
posted at 08:25:49
鬼滅、おもしろいなあ。いまのシーズンでは遊郭に潜んで「姫」を座敷名につける鬼が傍若無人に人を殺して回ってたのを、天元さんたちがきて守勢にまわると「みんながあたしをいじめるの〜」と騒々しく泣き叫ぶの。するとその分身の醜い兄鬼が出てきて「よくも頭の悪い妹をいじめたなー」と言うんだけど
タグ：
posted at 08:25:48
『被害のほとんどない世代に効果の怪しいワクチン打つ意味って本当にあります？』
『0歳からたくさん打つのに～って言う人居ますけど、なら無意味な注射は減らした方が良くないです？
archive.ph/unpbf』
『無意味』『効果の怪しい』ワクチンは市販されませんよ、反ワクチンさん
タグ：
posted at 08:21:30
『あなたは接種していない医療関係者も、多数いるのをご存じないのかな？』
『いや、医師ならご存じですよね。こりゃ、失礼
archive.ph/1tBWb』
2021年の6月時点で2回接種90％ですけど、新型コロナワクチン未接種者の医療従事者が『多数』とは？
こりゃ、失礼
www.mhlw.go.jp/content/109000... pic.twitter.com/KuxRij5wOi
タグ：
posted at 08:17:01
非公開
タグ：
posted at xx:xx:xx
引用『東北大は自ら高校に調査を行ったところ、英語民間試験利用の賛成が8％にすぎないことを重く見た。』
一度大失敗をしているのだから、共通テストへ「情報」を導入するにしても、事前にこういう調査をしっかりやってそれが公表されていないとおかしいんですよね。
タグ：
posted at 08:07:07
Buitengebieden @buitengebieden
“Sorry mom, no hug, it’s freezing cold outside here” 😂 pic.twitter.com/hcn8CdsvRk
タグ：
posted at 07:59:15
『ただ撲滅する事は今現在では不可能ですね。SARS MARSも現存します』
まあ、火星と中東呼吸器症候群を間違える時点で、リプライする必要はない人ですよ
twitter.com/saws_skmdrt/st... pic.twitter.com/EB6y6mztMK
タグ：
posted at 07:55:17
Buitengebieden @buitengebieden
The dead leaf butterfly (Kallima inachus) is a master of beauty and disguise.. pic.twitter.com/lV7LLgXIM2
タグ：
posted at 07:45:12
きょう24日付東京版１面です。
23日の将棋・ ＃王将戦(スポニチ主催)第２局第２日(大阪府高槻市・山水館)は、＃藤井聡太 竜王(19)=４冠=が ＃渡辺明 王将(37)=３冠=に98手で勝利。開幕２連勝としました。第３局は29、30日に栃木県大田原市の「ホテル花月」で行われます。 pic.twitter.com/7YHb1kmKtO
posted at 07:40:38
@hollyonthehill So beautiful! I had something similar on my deck back in December. pic.twitter.com/lPdyU94b1V
タグ：
posted at 05:05:18
それでも職場を巻き込もうとするな直接本人を訴えろ、と１０年以上前から繰り返し言ってるわけでして。引きずるも引きずらないもそんなもの今更関係ないんですよ。クレーム入って当たり前、などと思ってる人がいる間は改善しない。 twitter.com/onnawaturaiyo1...
タグ：
posted at 03:55:44
Why you should avoid imperative state in your programs
ht @c0nc0rdance/@mattparlmer pic.twitter.com/o8oW5N8bDy
タグ：
posted at 03:08:32
結論は変わらないと言っているのですよ。
東電原発事故で放出された放射性ヨウ素によって甲状腺がんは増えていないし、今後も増えない twitter.com/tswaxaatkyosss...
タグ：
posted at 01:02:55
ごまふあざらし(GomahuAzaras @MathSorcerer
やっぱり時代はJuliaっぽい twitter.com/karaage0703/st...
タグ：
posted at 00:34:04
あと役所の権威を借りに有識者風にすり寄ってくるあたりを便利に使うのもこける原因の一つ。1人去れば1人入ってくるように役人にとっての便利屋委員がいるんだよ。営業に使われているだけ。お国の施策が食いものにされるのだがお友達になると目が曇る。
タグ：
posted at 00:12:30
教育ログを本人同意で適法化できるという専門家を起用しないこと。かつての権威も今の施策に大丈夫とは限らない。役人のスキルが問われるところだ。一応記録のために言っておいた。わからんことは黙殺するあたりも知っている。置物歴含めて30年みてるからな。
タグ：
posted at 00:09:57
nを与えて確率を計算するプログラムを書け、みたいな感じでプログラミングの試験問題にしようと思ったけど、ちょっと短時間では厳しそうなので没に。競プロ勢なら一瞬なのかな。n=9が全然帰ってこないんだが、ワイのプログラムが悪いのか…
タグ：
posted at 00:04:43
