黒木玄 Gen Kuroki(@genkuroki)/2022年02月06日

#Julia言語以下のリンク先の計算を並列化無しでやっていたら1時間近くかかったと思います。

計算がはやく終わることが分かっていると、

　遊びでやる計算

が増えます。遊びでやる計算は数学的な事柄を理解するために非常に役に立つことが多い。

計算速いは正義！

github.com/genkuroki/publ... twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 00:13:11

#Julia言語並列化のためのコードは添付画像の赤線の部分の3行だけです。ものすごくシンプル！

そこに出て来るThreads.{nthreads, @ threads, threadid}について公式ドキュメントで調べて、rand!を使うメモリ割り当ての節約法を理解すればすぐに使える方法です。

github.com/genkuroki/publ... pic.twitter.com/Y9HMajM3s8

タグ： Julia言語

posted at 00:19:02

Massimo @Rainmaker1973

Dendrosenecio kilimanjari is a giant groundsel found atop Mount Kilimanjaro in Africa. The giant grounsels of the genus Dendrosenecio evolved, about a million years ago, from a Senecio that established itself on Mount Kilimanjaro [read more: ow.ly/A2PL50Bdw1F] pic.twitter.com/FdClQDLra3

タグ：

posted at 00:19:28

#Julia言語確率分布distのサイズnのサンプルを、ループ内で

X = rand(dist, n)

で作ると、ループを回した回数分だけ、これにメモリが割り当てられてしまいます。そうならないようにするためには using Randomしておいて、rand! を使います。

これを知らないと計算時間が大幅に無駄に伸びます。

タグ： Julia言語

posted at 00:22:05

#Julia言語中身が独立に互いに無関係に干渉せずに実行可能なforループはその前に Threads.@ threads を追加するだけで、一瞬で並列化できます。

しかし、ループの内側でrand(dist, n)を使うとメモリ割り当ての増大で死ぬ。

そうならないようにする処方箋が上で紹介されているわけです。

タグ： Julia言語

posted at 00:25:44

#Julia言語 GPUを使えばこの手の単純な繰り返しはさらに速くできるのですが、コードを書くのは少し面倒になります。

CuArrayとそれへのブロードキャストで済む場合は簡単なのですが。

Juliaでは並列化やGPU計算が比較的楽にできます。

タグ： Julia言語

posted at 00:28:19

@yaizawa

@nekoruri 真似してOpenSCADでコーディングしてみましたw
twitter.com/yaizawa/status...

タグ：

posted at 00:54:14

まぼろし博覧会@_fushiginamachi に言ってきた。
到着早々、セーラちゃん　@maboroshimusume　に遭遇。

共産趣味的部分を中心に、写真をアップする pic.twitter.com/S02EVpi0uY

タグ：

posted at 07:31:37

一瞬、「中核と革マルがない」と思ったら、左に仲良く並んでいました😀 pic.twitter.com/nGFjAi7NP5

タグ：

posted at 07:32:54

pic.twitter.com/xTx3M5F8xS

タグ：

posted at 07:33:24

pic.twitter.com/Gub3ISWZy5

タグ：

posted at 07:33:41

なぜか上田哲 pic.twitter.com/aK1YpXpwDn

タグ：

posted at 07:34:32

pic.twitter.com/f0SmEOKpJq

タグ：

posted at 07:35:16

日大闘争は、三島キャンパスにも及んでいたようだ。 pic.twitter.com/WfQLUKoWLm

タグ：

posted at 07:36:25

pic.twitter.com/DCi0DatRL0

タグ：

posted at 07:37:27

文革時のポスターと思われる。
「美帝」は「アメリカ帝国主義」だろうけど、

打倒？修　って何？おそらく、党内の修正主義を打倒する、という意味だろうけど、読めない。 pic.twitter.com/LfZSUsFS4Q

タグ：

posted at 07:39:55

非公開

タグ：

posted at xx:xx:xx

eval parse @evalparse

It's interesting to contrast how R, Python approach this.

Julia does not have DataFrames not html tools built in but packages are composable. So TableScraper.scrape_tables(URL)[1] |> DataFrame.

1/3 #julialang twitter.com/benhamner/stat...

タグ： julialang

posted at 08:48:07

Jonathan Chen @jowch

"So how's your research going?" 🤣

Created with @carbon_app pic.twitter.com/JgPm2kLPSI

タグ：

posted at 09:10:28

Jonathan Chen @jowch

Code is now available at github.com/jowch/Wordle.jl. More fun stuff to come!

タグ：

posted at 09:39:08

ナオミ @nao73714

過剰診断は、必要はなかったはずの手術や治療を選ばざるを得なくなることも問題だけど、手術しなければよいというものでもないです。診断そのものの心身と生活への影響は決して小さくはないです。

タグ：

posted at 09:53:38

ナオミ @nao73714

診断自体が重大な害をもたらさないよう、疾患の概念や観念を変えていこうという向きもあるけれど。どの疾患でもどの世代でも器用に付きあっていくことを求めるのも酷で、社会もそこまで変わってもいない。

タグ：

posted at 09:57:06

飯島明子 @a_iijimaa1

これが想像できない人は少なくないのですよね…。 twitter.com/nao73714/statu...

タグ：

posted at 10:07:04

sekkai @sekkai

知念実希人　物語り @MIKITO_777

実効再生産数がピークアウトしました！！って、1を超えている限りは感染者数は増え続けるのですがそれは☺️ twitter.com/sf_satoshifuji...

タグ：

posted at 10:28:38

違うんだ。

実効再生産数は１を切らない限り、感染は拡大し続けるから、『実効再生産数のピークアウト』を安心材料にしちゃダメなんだ。

実効再生産数が1を切って初めて『ピークアウト』という言葉を使うべきなんだ。 twitter.com/SF_SatoshiFuji...

タグ：

posted at 11:01:50

Ishida the Brain Dam @tbs_i

知念実希人　物語り @MIKITO_777

マルチコアとマルチスレッディング対応言語の普及で、「全部のコア使用率を100%に張り付ける」という楽しみ方が出てきたw twitter.com/genkuroki/stat...

タグ：

posted at 11:24:18

セーラちゃん　まぼろし博覧会 @maboroshimusume

あと、この実効再生産数の1.23も実は間違っています。

これはたしか世代時間5日での計算ですので、世代時間が2日ほどのオミクロンの実効再生産数を正確に表した数値ではありません。

現在の東京の実効再生産数は1.1程度で、あともう少しで1を切ってピークアウトが見えるレベルまでは来ています。 twitter.com/SF_SatoshiFuji...

タグ：

posted at 11:32:30

@sekibunnteisuu @_fushiginamachi ありがとうございました🥰‼️
まぼろし博覧会は、これからも日々進化増殖してお待ちしてます🤗💓💕🎊

タグ：

posted at 11:35:55

上松正和 @Uematsu1987

小菅信子 Nobuko M KOSUG @nobuko_kosuge

明確な誤り。
実効再生産数は「1人の感染者が平均で何人にうつすか」がなので1より大きいと感染が広がり続けます。
データを見るのは大事だけど、データを正しく解釈する力はもっと大事。数字は嘘をつかないけど、数字で嘘をつく人は多いから。
奈良県は大丈夫かな...。 twitter.com/SF_SatoshiFuji...

タグ：

posted at 11:53:30

某氏のブログ削除の件について。
これ、某氏が「全国的リアル著名人」だから注目が集まっただけだよね。この種の経緯で「削除」したりされたりした人たちは、2011年3月11日以来だけでか・な・りいたよ。

タグ：

posted at 11:55:05

#統計

標本平均と標本分散の分散と共分散は母分散σ²と母歪度κ₃と母尖度κ₄で書けます。κ₃=E[((X-μ)/σ)³]、κ₄=E[((X-μ)/σ)⁴]-3で、E[exp(t(X-μ)/σ)]=t²/2+κ₃t³/3!+κ₄t⁴/4!+O(t⁵).

var(X̅)=σ²/n
cov(X̅, S²)=(κ₃/n)σ³
var(S²)=(κ₄/n + 2/(n-1))σ⁴

正規分布でκ₃=κ₄=…=0.

タグ：統計

posted at 11:57:36

小菅信子 Nobuko M KOSUG @nobuko_kosuge

#統計 Z=√n(X-μ)/σとS²/σ²については

E[Z] = 0
var(Z) = 1
cov(Z, S²/σ²) = κ₃/√n
var(S²/σ²) = κ₄/n + 2/(n-1) ～ (κ₄ + 2)/n

正規分布ではκ₃=κ₄=0.

タグ：統計

posted at 12:03:00

小菅信子 Nobuko M KOSUG @nobuko_kosuge

まあ、削除させたら、「削除しろ」側も削除しろって書いた「事実」を削除したほうがいいかもね。いや、「削除しろの理屈」で言うとだけどね。

タグ：

posted at 12:09:05

人の生きた「証」をそう簡単に削除しろしろ命じるのは、どうかと思うね。

タグ：

posted at 12:12:44

きっちー @rikeden_net

わわ、自分史上初めての100いいね行きました。動作原理とか気になる方はこちらの動画もぜひご覧ください。
只今絶賛後編製作中です。
nico.ms/sm39526753

タグ：

posted at 12:13:08

#Julia言語

CPUのコアが全部100%にはりつき、パソコンのファンがシューっとなるのは楽しい。😊 twitter.com/tbs_i/status/1... pic.twitter.com/f88vb6xvkb

タグ： Julia言語

posted at 12:20:37

#統計 #統計 Z=√n(X-μ)/σとW=√n(S²/σ²-1)については

E[Z] = 0
E[W] = 0
var(Z) = 1
cov(Z, W) = κ₃
var(W) = κ₄+ 2n/(n-1) ～ κ₄ + 2

(Z, W)の分布は、平均が(0,0)、分散共分散行列が

1 κ₃
κ₃ κ₄+2

の2変量正規分布に収束する。

このスレッド中の後半のグラフでこれがよく見えている。

タグ：統計

posted at 12:31:41

#統計 1つ前のツイートの内容は、母平均と母分散の不偏推定の漸近挙動というザ・統計学的な話題。

あまりにも基本的な話題なので、コンピュータでグラフを描いて確認しまくっておかなければいけなかったのですが、ダメすぎる私はそれをやっていなかった。

反面教師扱いして頂ければうれしいです。

タグ：統計

posted at 12:40:22

なっちょ @toshiNpapa

ネタわからない人用にコッチがレオンちゃんです。 pic.twitter.com/WbYwz6IbJB

タグ：

posted at 12:41:05

「万年筆マネー vs. 又貸し説」を唱えている学者先生達の発言が馬鹿にされるようになった方が良さそう。

マンキューに「又貸し説」とやらが書いてあるらしい(←皮肉)。

いやあ、正常性バイアスが強すぎるせいで、こんなに馬鹿な奴らが増えていることに最近になるまで気付かなかった。

タグ：

posted at 12:53:25

「普通に標準的な教科書を読め」でお終いの話題のはず。

タグ：

posted at 12:54:13

藤井聡教授が顔写真を連ねている人達が怖すぎる！

こういう人達と名前を連ねてしまうような人物には日本の政策に関わらせることをやめさせるべきだと私は思います。

www.werise.tokyo/declaration/ pic.twitter.com/nLnOiHh4Gl

タグ：

posted at 14:12:02

関連
↓ twitter.com/sekkai/status/...

タグ：

posted at 14:13:02

藤井聡教授のウィキペディアより

ja.wikipedia.org/wiki/%E8%97%A4...
【安倍内閣・内閣官房参与～などを歴任した。2020年設立された「新型」コロナウィルスの存在を否定する（従来の4種のコロナウイルスに追加される5種類目のコロナウイルスとして定義する）WERISEにて共同宣言者を務める[1][2]。】

タグ：

posted at 14:19:06

こまねこさん @koma2023

なるほどなー twitter.com/genkuroki/stat...

タグ：

posted at 14:19:53

安倍内閣・内閣官房参与だった人物が○海氏などのめちゃくちゃ怖い人達と名前を連ねていることの恐怖感！

仮にこういう人物が再度政府で重用されるようになり、その人脈で怖い人達の影響力が政府に及びようになったらどうなる？

全然洒落にならない怖い話。

タグ：

posted at 14:22:03

藤井聡氏関連
↓ twitter.com/mikito_777/sta...

タグ：

posted at 14:33:35

藤井聡氏は【2020年設立された「新型」コロナウィルスの存在を否定する（従来の4種のコロナウイルスに追加される5種類目のコロナウイルスとして定義する）WERISEにて共同宣言者を務める[1][2]。】 ja.wikipedia.org/wiki/%E8%97%A4...

めっちゃまずい人達と名前を連ねていて全然笑えない。

藤井聡氏関連
↓ twitter.com/uematsu1987/st... pic.twitter.com/ptfENi81QP

タグ：

posted at 14:37:48

@kuri_kurita

小学生の算数でしばしば大人が問題にする部分ですが、問題が「何人すわれますか」なので答は「人」になります。長いすが７つあるので7人、その隣にまた7人...と数えていくと7+7+7+...となります。7を9回たす→7x9と書けます。足し算と掛け算、数式に出会う大切な初期なので色んな考えを試しましょう♪ twitter.com/i_diction/stat...

タグ：

posted at 14:54:21

#統計このスレッドでは、まずは易しい1標本のStudentのt検定について理解するためのプロットをソースコード付きで紹介しています。

「検定やP値の使用は時代遅れ。信頼区間を使いましょう」のように思った人は信頼区間が検定やP値で作られる区間であることを理解していません。それは非常にまずい。 twitter.com/genkuroki/stat...

タグ：統計

posted at 15:14:01

#統計

(I) データxについて、パラメータθの信頼係数1-αの信頼区間を与えること

と

(T) 値θ₀について、仮説「θ=θ₀」を有意水準αで棄却するデータxの範囲を与えること

は数学的に同値なので、信頼区間の理論は仮説検定の理論と完全に同じだとみなされます。

タグ：統計

posted at 15:19:38

#統計検定と信頼区間の表裏一体性は以下の図に要約されます。

信頼区間はデータから検定で棄却されないモデルのパラメータの値の範囲です。

これは勉強量を節約したい人にとって朗報！

検定やP値について理解していれば、信頼区間についてもすでに理解していることになります！ twitter.com/genkuroki/stat... pic.twitter.com/itxj4F94hs

タグ：統計

posted at 15:25:45

#統計というわけで、1標本の場合のt検定について理解することは、そのまま平均の信頼区間について理解することにもなっているわけです！

平均値の信頼区間について理解を深めるためには、平均値に関する検定について理解を深めることに集中すればよい。

タグ：統計

posted at 15:28:18

#統計 P値に関するASA声明の翻訳 biometrics.gr.jp/news/all/ASA.pdf がどのような経緯で行われたかについても書いてある

www.jstage.jst.go.jp/article/jjb/38...
ASA声明と疫学研究におけるP値
佐藤俊哉
2017

の最後の段落でも、信頼区間と検定の表裏一体性の指摘が引用されています。 twitter.com/genkuroki/stat... pic.twitter.com/Gc4uSQ7wdy

タグ：統計

posted at 15:38:25

#統計ところが「P値以外のアプローチ」として「信頼区間」を挙げる言説が流れてしまっていて、「信頼区間もまたP値を使う方法そのものである」という情報に触れ難くなってしまっている。

P値を理解していれば信頼区間も理解したことになって勉強時間を大幅に節約できるメリットが潰されている。

タグ：統計

posted at 15:41:46

#統計 ○○の効果がないことは△△モデルの立場で「θ=0」と書けるとき、仮説「θ=0」のP値のみを考えるから、P値の道具としての有用性が下がっているのです。

0以外の任意の値θ₀についても「○○の効果はθ₀である」=「θ=θ₀」という仮説のP値を全部考えれば信頼区間の話に繋がって行きます。

タグ：統計

posted at 15:46:05

akiranya @akirana

うわ。私的に、これはないわ。 twitter.com/genkuroki/stat...

タグ：

posted at 15:46:19

aramaki @aramaki34

藤井って人、トンデモさんのお仲間だったのかあ、、、
肩書が京都大学大学院教授？！ twitter.com/genkuroki/stat...

タグ：

posted at 17:33:44

タクラミックス @takuramix

怒らない白井（仮想軒と須雷堕は除く） @tatsuvar

「哲学」に妙な期待を抱いている人を目にする。
今の社会とか科学とかから「哲学が失われている」などと嘆いて見せる人が居て、その言わんとするところは分かるし共感するところもあるけれども、失われたものも必要になるものも哲学ではない…というか、哲学に変な期待を抱くのは危険だと思う。

タグ：

posted at 17:49:05

Julia日本語記事まとめ @julia_kizi

.@Osakata6u パワハラかどうかはよく分かりませんが、無意識に生徒を見下していますね。「これくらい覚えられないの？」（バカじゃない.. togetter.com/li/1833142#c10...

「遊んでやろうと言うのでお言葉に甘えて「掛け算の順序は暗記しないといけないのか..」togetter.com/li/1833142 にコメントしました。

タグ：

posted at 19:02:40

Julia日本語記事まとめ @julia_kizi

新しい記事がQiitaにアップされました！#Julia言語 #Julia日本語記事
 qiita.com/RhT/items/b9c9...👈

タグ： Julia日本語記事 Julia言語

posted at 19:34:23

新しい記事がQiitaにアップされました！#Julia言語 #Julia日本語記事
 qiita.com/RhT/items/d28e...👈

タグ： Julia日本語記事 Julia言語

posted at 19:34:25

Open Plant Pathology @OpenPlantPath

New blog post on OPP! Have you heard about Julia and want to know if it might be useful in your work as a plant pathologist? We've a new post showing how it compares to R for creating an AUDPC function. Check it out! openplantpathology.org/posts/2022-02-...

#JuliaLang #RStats

タグ： JuliaLang RStats

posted at 19:38:40

#統計訂正：Zの定義中のXは本当は標本平均X̅です。

κ₃, κ₄を母歪度、母尖度とすると、添付画像に続けて、

cov(Z², W) = κ₄/√n
var(Z²) = κ₄/n + 2

期待値1のZ²=n(X̅ - μ)²/σ²と期待値0のW=√n(S²/σ² - 1)の共分散とZ²の分散が母尖度κ₄ (正規分布で0になるバージョン)で決まる。 twitter.com/genkuroki/stat... pic.twitter.com/HT5BcTdSt2

タグ：統計

posted at 19:43:17

なのわに @naniwa_osaka9

人口密度順に市町村を塗っていき、国内人口の半分になったところで止めた結果がこちら。 pic.twitter.com/TwkOPHzwev

タグ：

posted at 19:44:15

#統計例えば、以下のリンク先の場合には、母尖度=-1.2なので、Z²とWの共分散は負になり、Z²とX²の共分散も負になります。X²が増えるとZ²は小さくなる傾向になるはず。

実際、以下のリンク先の右上のプロットを見ると、縦軸のX²が大きくなると、Zの横の広がり方は狭くなっています。 twitter.com/genkuroki/stat...

タグ：統計

posted at 19:49:28

なのわに @naniwa_osaka9

あなたのお住まいは？

タグ：

posted at 20:14:14

なのわに @naniwa_osaka9