黒木玄 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年07月07日(木)
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#Julia言語 添付画像の例を理解すれば多分やりたいことをできると思います。 twitter.com/nankyoku_dayo/... pic.twitter.com/iOoIzwFwtD
タグ： Julia言語
posted at 00:55:11
#Julia言語
for文を駆使して自力で何とかできるなら、これで問題は解決だと思います。
x,y,zのデータの形式によっては、必要な変換を楽に行う方法があることもあります。データの形式の詳細だけではなく、テストデータも用意すれば挑戦する人がきっと出て来ると思う。データの形式の情報が重要。 twitter.com/nankyoku_dayo/...
タグ： Julia言語
posted at 09:11:43
#Julia言語
入力が3つの同じ長さのベクトルxx,yy,zzのとき、座標(xx[i], yy[i])に色zz[i]でドットをプロットする方法→添付画像①②
データが添付画像①の形式(xx,yyがmeshgrid)ならfor文なしにheatmapで使える形式に変換できます→添付画像③
ソースコードはALTにある。 twitter.com/nankyoku_dayo/... pic.twitter.com/6x2tXTmanE
タグ： Julia言語
posted at 09:43:20
#統計 「わかる人にはわかる」ことは実に多い。
確率パラメータpはオッズu=p/(1-p)でp=u/(1-u)と表示されることが多い。x = log u のとき、p = logistic(x)=1/(1+e⁻ˣ)でロジスティックモデルの話になる。
添付画像中の1/θがpのオッズになっている。
ガンマ函数やベータ函数も「理解」が重要。 twitter.com/othinus1123581... pic.twitter.com/1yzFLhzikG
タグ： 統計
posted at 11:09:55
#統計 負の二項分布については以下が基本的。
* αが正の整数の負の二項分布はベルヌーイ試行で解釈できる。
* 負の二項分布の連続時間極限でガンマ分布が得られる。
* 期待値パラメータがガンマ分布に従うポアソン分布は負の二項分布そのものである。(ただしパラメータはオッズに変換される) pic.twitter.com/b4gslX7Oo5
タグ： 統計
posted at 11:10:04
#統計 成功確率パラメータpがベータ分布(共役事前分布)に従う二項分布には、「ベータ二項分布」と呼ばれています。
「負の超幾何分布」と呼ばれることもあり、ポリアの壺での解釈を持ちます。
en.wikipedia.org/wiki/Beta-bino....
タグ： 統計
posted at 11:24:31
#統計 cdf(Binomial(n,p),k)=ccdf(Beta(k+1,n-k),p)は具体的に書くと、
Σ_{i=0}^k binom(n,i) pⁱ(1-p)ⁿ⁻ⁱ
= ∫_p^1 tᵏ(1-t)ⁿ⁻ᵏ⁻¹/B(k+1,n-k).
これの左辺は基本特殊函数でコンピュータで効率的に実装されていて、左辺の和を地道に計算するより、圧倒的に速いです。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 11:41:33
#統計 二項分布の累積分布函数のコンピュータでの実装では、二項分布とベータ分布の累積分布函数の関係がもろに使われていたりします。
二項分布とベータ分布の累積分布函数の関係については以下のリンク先を参照。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 11:44:47
#統計 で、非常に役に立つ
Σ_{i=0}^k binom(n,i) pⁱ(1-p)ⁿ⁻ⁱ
= ∫_p^1 tᵏ(1-t)ⁿ⁻ᵏ⁻¹/B(k+1,n-k)
は、出目が1と0のベルヌーイ試行を、出目が0から1の連続的な一様分布に従う乱数の場合に持ち上げれば、確率論的直観で明らかな公式になります。
ベータ分布は一様分布の標本分布から出る。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 12:06:06
#統計 #Julia言語
負の二項分布→ガンマ分布
と
ベータ二項分布→負の二項分布
の数値的確認の視覚化。
こういうグラフは必ず描いておいた方がよい。 pic.twitter.com/JXDeQZOLuO
posted at 12:33:18
数学が絡むとその瞬間に教科書に書いてあることをフォローすれば理解が進むという考え方は通用しなくなると思う。
教科書に書かれていないこと(要するに試行錯誤)をノートにたくさん書き、教科書に描かれていないグラフを大量にプロットして様子を眺めてみる、というような地道さが必須になる。
タグ：
posted at 12:47:24
#統計 訂正：
❌p=u/(1-u)
⭕️p=u/(1+u)
符号はよく書き間違う。
この話題に興味ある読者達は勝手に訂正してくれるとわかっているので、どうしてもテキトーになる。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 12:58:27
#統計 X～Gamma(a,1), Y～Gamma(b,1)のとき、
P=X/(X+Y) ～ Beta(a, b)
となることは、本質的にベータ函数をガンマ函数で表す公式と同じ話。Pのオッズ
U=P/(1-P)=X/Y
はベータプライム分布 en.wikipedia.org/wiki/Beta_prim... に従います。続く
タグ： 統計
posted at 13:06:02
石橋湛山さんが草葉の陰で泣いている。 #経済 twitter.com/kuri_kurita/st... pic.twitter.com/Ztbs1UBgXe
タグ： 経済
posted at 13:30:32
New package: RealPolyhedralHomotopy v1.0.0 announced #JuliaLang
RealPolyhedralHomotopy: A package for finding real roots of systems of polynomial equations using polyhedral homotopy.
Registration: github.com/JuliaRegistrie...
Repository: github.com/klee669/RealPo...
タグ： JuliaLang
posted at 13:33:18
日銀に文句を言っているようでは、怖くて絶対に投票できません。
この手の人が国会議員になって、金融引き締めを推進してそれが実現してしまうと、金融引き締めによる経済崩壊によって、年間数千人単位で無駄に人が死ぬことになるでしょう。 twitter.com/obatakimiko_go...
タグ：
posted at 14:02:17
即効でひとを年間数千人オーダーで無駄に死なせる政策には、
* 日銀に金融引き締めを行わせる。
* 反ワクチン勢力に与する政策を実施する。
があり、差別的に人を苦しませることには
* 放射能デマに与すること
がある。どれか1つだけでも当てはまったら、その候補に投票するべきではないと思う。
タグ：
posted at 14:06:31
#統計 この辺も、色々読むと、妙なこだわりのある人達がいて、
❌Strene型ではなく、片側確率を2倍するClopper-Pearson型のP値を使うべきだ
のようなことを言う人達がいる。
⭕️トレードオフの問題でしかない
のに、妙なこだわりを強弁しても恥にならない世界が広がっている感じ。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:02:27
#統計 Rのbinom.test, fisher.test, poisson.testが表示するP値はどれもSterne型。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:03:56
#統計 RがSterne型のP値を採用していることについて、「この特殊な両側有意確率を求める方法に製薬企業が飛びつき、やがてそれがR等の統計ソフトにまで採用」と書いてあるのを見つけたことがある。😱
いやあ、びっくりした。😅 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:08:32
中央銀行による金融政策は経済政策のスーパーエース的な役割を果たしており、方針を変えると影響が巨大になります。
金融政策の雇用への影響は大きいので、より良い就職や転職を求めている人達は、金融引き締めを実施されてしまうと被害に遭う可能性が高いことを十分に認識して投票するべき。
タグ：
posted at 15:13:08
このタイミングで日銀に文句を言っている人達が苦しめることになる可能性があるのは、現在失業中もしくはこれから失業する可能性のある人達です。
こういう知識は義務教育の範囲内でも習っている可能性があって、「有権者」的には非常に大事。
タグ：
posted at 15:15:48
#統計 他にも変なこだわりに関する奇妙な議論を目にすることもあって、「条件付き推論」に関するNP vs. Fisher的な議論にも設定の詳細を曖昧にしたままの奇妙な議論がある。
「科学的に重要な条件については気にし、そうでないものは気にしない」で済むはずの議論がおかしなことになっていると思う。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:36:11
#統計 それらの変なこだわりのある検定に関する議論は、以下のリンク先で紹介した現代的なP値の使い方に関する議論からも大幅にずれており、時代遅れな感じがするものになっていると個人的には思います。 twitter.com/genkuroki/stat...
タグ： 統計
posted at 15:41:38
東野篤子 Atsuko Higashin @AtsukoHigashino
余談ですが、この戦争開始後の衝撃的な経験として、あるロシア専門家の方に「あんな汚職まみれの国、一度ロシアに倒して貰って再建し直した方がいいんじゃないの」と直接言われたことがあります。
はっきり申し上げますが、ウクライナに汚職が蔓延していることは、ロシアの侵攻を一切正当化しません。
タグ：
posted at 16:05:45
これはマジでこの世の人間全員によく覚えておいてほしいんだけど、頭を強くぶつけた後に嘔吐をした場合は直ちに病院を受診すべきです
最悪の場合、無知は死につながります
news.yahoo.co.jp/pickup/6431792
タグ：
posted at 20:27:48
こういう問題、医療従事者が「ヤバい」って判断できるのはまあ当たり前なんだけど、一般に知らない人がいるのは仕方ないことだと思う
なので常に啓発などをする人間が必要になるんだよね
そして致命的なイベントに繋がりかねない病態くらいはまじで義務教育とすべきだと思う(n回目
タグ：
posted at 21:54:12
比例については、アニメや漫画の将来を考えるとやっぱし赤松さんに入れることになりそう。
山田太郎さんと一緒になって頑張って欲しい。 twitter.com/kenakamatsu/st...
タグ：
posted at 23:36:29
惜しいな。第一文は正しいけど、第二文で緊縮経済政策だということが露呈している。
正しくは、国債発行で賄う。そのためには金融緩和を継続することは必須。これが共産党には理解できない。
消費税を完全に廃止しても年間25兆円くらいなので、国債で賄える。 twitter.com/yamazoejimusyo...
タグ：
posted at 23:46:05
