黒木玄 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月26日(金)
#超算数 詳細編は、合併(167)と、増加(168)とについて、具体的な操作とは絵の上にブロックやおはじきをおいて操作することであると考えているようである。「<参考>たし算について」(168)にある通り増加と合併は明確に異なるものと考えられているが、両者が持つ具体的な操作が相違するとは思われない。 pic.twitter.com/uYed59ctah
タグ: 超算数
posted at 00:37:19
#超算数 ちなみに奥付。現代の教師用指導書には作成に携わった人々のリストが載りません。ただし、この【編修】リストも、教科書のそれと違うのかどうなのか、調べていません。 pic.twitter.com/tyP7fzXejC
タグ: 超算数
posted at 00:40:38
@lotzun_DeuPol ああ、よっぽど自分の発言がまちがってるんじゃないかと不安だったんですねー。で、どっかの変な発言が自分を擁護してくれてると思って、大喜びしたんですねー。かわいそうに。
タグ:
posted at 03:10:26
どうしても趣味がJulia言語の方に行ってしまう(´・ω・`)大きな欠陥があるところ、不完全なところがまた愛おしいのですよ。#ロボット開発
www.mish.co.jp/applications/r...
タグ: ロボット開発
posted at 03:11:51
New package: AbstractPDEInterfaces v0.1.0 announced #JuliaLang
AbstractPDEInterfaces: Ecosystem for writing partial differential equation solvers
Registration: github.com/JuliaRegistrie...
Repository: github.com/vpuri3/Abstrac...
タグ: JuliaLang
posted at 05:02:27
@genkuroki こんなやり方があるんですね。追加したパッケージの依存関係で、ほかのパッケージのバージョンが変わるから、と思ってやり方を探ってたら、カーネル追加する方法がでてたので、それをやってました。
黒木さんの方が効率的ですね。ありがとうございます。
タグ:
posted at 05:19:03
#統計
xcelab.net/rm/statistical...
Statistical Rethinking 公式ウェブサイト
↓
xcelab.net/rmpubs/sr2/sta...
Chapters 1 and 2
Chapter 2. Small worlds and large worlds の全文を読めます。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 09:22:39
Julia言語は速いという話を聞いたので調べてみたら、この記事の実験ではfor文比較でpythonよりも100倍速いらしい。どうしてもpythonで長時間のfor文処理が発生してしまう場合にはその部分だけJuliaにやらせるのはありかも🤔 (並列処理もできる)
qiita.com/Ken-Kuroki/ite...
タグ:
posted at 09:48:39
このあたりのアドバイスを厳守した結果,爆速になりました!!! twitter.com/genkuroki/stat... pic.twitter.com/opAz4dSkLy
タグ:
posted at 10:31:31
#Julia言語 うぉっ!アロケーションの減らし具合がハンパない!くっそ速くなっている!すごすぎ。 twitter.com/physics303/sta...
タグ: Julia言語
posted at 10:48:38
不完全ガンマ関数を微分するってどうやるんだっけ?って調べていたらヤバい扉を開いてしまった。
おいおい、なんじゃこりゃ〜
mathworld.wolfram.com/MeijerG-Functi...
タグ:
posted at 11:45:14
#統計 Statistical Rethinkingより。これも私がいつも言っていたこと。
ベイズの定理を使っているかでベイズ統計か否かを区別するのは間違いです。
例えば、有病率、感度、特異度にベイズの定理を適用することを、ベイズ統計扱いすることは間違っています。
xcelab.net/rmpubs/sr2/sta... pic.twitter.com/7JpKdtDcq3
タグ: 統計
posted at 15:44:08
#統計 ベイズの定理は以下の2つのステップに分解されます。
リンク先の(2) p(x|y), p(y)が与えられているとき、p(x,y)=p(x|y)p(y)によってx,yの同時分布を作る。
リンク先の(1) x,yの同時分布p(x,y)から、p(x,y)=p(y|x)p(x) を満たすように、条件付き分布p(y|x)と周辺分布p(x)を作る。
続く twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:15:08
#統計 積分を使う場合は小学5年生レベルとは言えないのですが、最も易しい有病率、感度、特異度から陽性的中率を計算する問題は小5レベルの割合の問題に過ぎず、小5で割合について習って理解できた子に自信を付けさせるために利用できます。
(実際に利用しました😊) twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:22:54
#統計 Statistical Rethinkingでも、ベイズ統計と非ベイズ統計の類似性を強調していますね。
「頻度主義」と言わずに「非ベイジアン」と書いている所もよいと思いました。
xcelab.net/rmpubs/sr2/sta... のp.39より
↓ pic.twitter.com/aoWQXSQNhT
タグ: 統計
posted at 16:43:05
#統計 ベイズ統計の枠組みには、統計分析的に適切な数学的操作を作ることが易しい(確率論的に最も自然な操作が正解になることが多い)という利点があります。
データの欠損値の取り扱いもベイズ統計の枠組みでは概念的に非常に易しくなります。
xcelab.net/rmpubs/sr2/sta... のp.35より
↓ pic.twitter.com/PYfLy41g12
タグ: 統計
posted at 17:14:09
#統計 事前分布の任意性をベイズ統計を否定する理由に挙げることは、馬鹿げた態度であることも説明されています。
ほぼ同じことが、渡辺澄夫『ベイズ統計の理論と方法』にも書かれています。
この辺は「常識」扱いして良いことだと個人的には思います。
xcelab.net/rmpubs/sr2/sta... p.36より
↓ pic.twitter.com/Ab9upEg28f
タグ: 統計
posted at 17:25:33
julia言語 でこーゆベクトルを得たいのですが,手っ取り早い方法ありますか?
いまはJのサイズが5ですけど,Jのサイズが7とか9とかでも大丈夫なような実装をしようとすると,[1,1,...,1,k,1,...,1] の部分をどうしたらよいか分からないです. pic.twitter.com/Z4HIyKvG5D
タグ:
posted at 18:35:47
一番ナイーブな方法としては,
① onesでD次のベクトルをつくって,
② d番目を1個ずつ大きくしていって,
③ H[cum[d]-(d-2):cum[d+1]-d] に代入
だけど,もっと良い方法あったりしないかなと...
タグ:
posted at 18:41:51
記事を投稿しました! 固有ベクトル-固有値同一性(Eigenvector-eigenvalue identity)をJuliaで確認してみる on #Qiita qiita.com/cometscome_phy...
タグ: Qiita
posted at 18:57:37
新しい記事がQiitaにアップされました!#Julia言語 #Julia日本語記事
qiita.com/cometscome_phy...👈
タグ: Julia日本語記事 Julia言語
posted at 19:00:55
#数楽
nbviewer.org/github/genkuro...
genkuroki.github.io/documents/Calc...
の2.10節に、フルヴィッツのゼータ函数と対数ガンマ函数の関係に関する解説があります。
対数ガンマ函数の導函数は役に立つ基本特殊函数で、コンピュータでも効率的に実装されており、単に純粋数学的な興味を超えた余得があります。 twitter.com/10ios4/status/...
タグ: 数楽
posted at 19:19:40
#統計 一般化されたゼータ函数の概念は渡辺澄夫『ベイズ統計の理論と方法』の第4章で使われており(汚く計算すれば避けられる)、ゼータ函数の積分表示(Mellin変換表示)に慣れておくことも余得の1つになります。
関連の解説
↓
genkuroki.github.io/documents/2016...
nbviewer.org/github/genkuro...
タグ: 統計
posted at 19:24:50
Archerfish shoot their prey down with water droplets from their mouths. They are remarkably accurate and can bring down insects and other small animals up to 3 m above the water's surface
[📹 buff.ly/3RfpGqz]
[read more: buff.ly/3Kn98uh]
pic.twitter.com/QjOXiWX3H5
タグ:
posted at 20:38:33
@Rainmaker1973 Related species, one one land too and this is their favourite prey pic.twitter.com/JjpGeRY3GN
タグ:
posted at 20:51:29
ごまふあざらし(GomahuAzaras @MathSorcerer
既存のメソッドを上書きする type piracy なコードを書いてしまってある種の罪悪感を持つ一方で、自分自身が作ったパッケージだし使い捨てだしいいか📦♪といういい加減さをもつスキルを身につけた。
(設計はちゃんとしよう)
#Julia言語
タグ: Julia言語
posted at 21:01:23
「現在、鳥取県の条例改正により、鳥取県で有害図書に指定されるとAmazonで販売停止になることが発覚いたしました。…
有害図書の問題に関してここで1度考えておくべきではないでしょうか」
「ラジオライフ」2022年10月号の有害図書に関する記事 www.sansaibooks.co.jp/info/%e3%80%8c...
タグ:
posted at 21:37:25
@shigepong #統計 不完全ガンマ函数の難しい方の偏導函数の計算と #Julia言語 での実装例が
nbviewer.org/gist/genkuroki...
にあります。参考になれば幸いです。
一般化超幾何函数が必要ですが、計算自体は大学1年生レベルです。 twitter.com/genkuroki/stat... pic.twitter.com/k3Mo5XN2Dd
posted at 22:19:44
#統計 東大出版会の『統計学入門』を信頼できる教科書だと思ってしまうと、ミスリーディングな解説を受け入れてひどいことになるので要注意。
信頼区間、尤度、パラメトリック、…、基本概念達についてことごとくミスリーディングな解説が並んでいる。
信頼区間の解説のまずさの説明に続く。 twitter.com/tachnopolis25/...
タグ: 統計
posted at 23:05:06
#統計 添付画像は『統計学入門』での信頼区間の解説部分です。(11.36)に
P(L≦θ≦U)≧1-α
という不等式があります。実用的には出来るだけ
P(L≦θ≦U)≈1-α
になるようにします。例えばα=5%ならばこれは
P(L≦θ≦U)≈95%
という条件になります。続く pic.twitter.com/X7KUi6UZ4H
タグ: 統計
posted at 23:05:10
#統計 区間[L, U]は95%信頼区間なので、式
P(L≦θ≦U)≈95%
は
データを生成した確率分布を決めるパラメータθが
95%信頼区間に含まれる確率は約95%になる
ことを意味します。元の(11.36)では「95%以上になる」になりますが、P(~)が確率を表すことに変わりはありません。続く pic.twitter.com/vjykr6OOcg
タグ: 統計
posted at 23:05:19
#統計 (11.36)に確率を意味するP(~)が含まれているので、95%信頼区間の95%は区間にθが含まれる確率(の近似値や下限)を意味することは、確定しています。
まず、この点を押さえておく必要があります。
95%信頼区間の95%は確率です。
これを誤りだと言う人はどうかしています。続く pic.twitter.com/bRrUOKSsH0
タグ: 統計
posted at 23:05:22
#統計 受け入れてしまうと次のように考えてしまうだろう:
❌仮に現実の母集団から標本を何度も取り直して95%信頼区間を計算し直せば、その中で、現実の母集団に関する真の値が含まれる区間達の割合は95%になる。
これはひどい誤りなので、添付画像の説明はミスリーディングだと言わざるを得ません。 pic.twitter.com/E0ds9jCeIz
タグ: 統計
posted at 23:05:28
#統計 代替となる正しくてかつ実用的な信頼区間の解釈の仕方(confidence intervalではなく、compatibility intervalと解釈すること)については、最新の論文で滅茶苦茶わかりやすく解説されています。
Amrhein-Greenland June 3, 2022
journals.sagepub.com/doi/10.1177/02... pic.twitter.com/It7TTbecsr
タグ: 統計
posted at 23:11:17
#統計 より正確に言えば、現実の母集団から取得したデータから、信頼区間を計算すると、数学的フィクションである統計モデルと現実で取得したデータの整合性(consistency、両立性、compatibility)がある程度分かります。
journals.sagepub.com/doi/10.1177/02... pic.twitter.com/hLh6OdUajk
タグ: 統計
posted at 23:24:30
#統計 論文 journals.sagepub.com/doi/10.1177/02... では、以上で私が「統計学のモデル」とよんだものを background model とよんでいます。
統計モデルを無視せずにきちんと考えることが、現代では主流の考え方になっています。
タグ: 統計
posted at 23:31:00
日本におけるオンライン・ハラスメントの現状と対策:Twitterでの女性記者のツイート「炎上」を例に
f1000research.com/articles/10-11...
査読論文、著者が宣伝していたのを見かけ、興味深く読んだのですが…。
かなり問題のある論文に見えます、素人目には。
明らかな瑕疵が多数あり、少なくとも1つは致命的。
タグ:
posted at 23:46:52
#統計 博士号持ちなら、私が繰り返し紹介している短くてクリアな論文 journals.sagepub.com/doi/10.1177/02... くらい、すぐに読めるはずです。
タグ: 統計
posted at 23:50:33
問題はいろいろあるのですが、いちばん深刻なのは、用語定義の曖昧さ。
オンライン・ハラスメントについての論文なのに、「オンライン・ハラスメント」という用語の定義がブレていて、一貫していません。
ハラスメントの解釈や考え方より前に、論文としてきびしいです。
f1000research.com/articles/10-11...
タグ:
posted at 23:51:14
>オンライン・ハラスメントとは「インターネット上で行われる、扇情的なコメントやヘイトスピーチの繰り返しの投稿(「荒らし」「トロール」)、サイバーストーカー、身体的な脅迫、同意なしに性的に露骨な画像を公開すること(「リベンジポルノ」)、個人情報を公開する(略)などの行為」を意味する pic.twitter.com/9ml9OVTGw6
タグ:
posted at 23:55:18
これが本論文で宣言されている定義です。PEN Americaの定義を採用したものです。
身体的脅迫、リベンジポルノ、個人情報晒しなどの犯罪行為はいいとして、問題は「扇情的なコメントやヘイトスピーチの繰り返しの投稿」です。
そう、オンライン・ハラスメントの定義に「繰り返し」を入れています。 pic.twitter.com/CqvSx4MupY
タグ:
posted at 23:59:06