黒木玄 Gen Kuroki(@genkuroki)/2020年02月29日

#Julia言語 v1.4.0-rc2.0 のスレッド並列を試してみました。v1.3以降のJuliaではrand()が遅くなっていて、以前の速さで乱数を生成するには

rng = Random.default_rng()
for ～
～rand(rng)～
end

のように書く必要があります。詳しくは

nbviewer.jupyter.org/gist/genkuroki...

を参照。

タグ： Julia言語

posted at 00:10:56

#Julia言語

v1.3以降からスレッド並列が非常に便利に使えるようになったと思うのですが、私よりもその手のことに得意な人達が簡単な例を沢山みせてくれると助かります。

nbviewer.jupyter.org/gist/genkuroki...
Distributed and Threads examples (Julia v1.4.0-rc2.0) pic.twitter.com/HjXnuQxCxn

タグ： Julia言語

posted at 00:17:53

#Julia言語 rand(rng) スタイルでないrand()での遅くなり具合が v1.3.1 より v1.4.0-rc2.0 の方が緩和しています。

nbviewer.jupyter.org/gist/genkuroki...
v1.3.1

と

nbviewer.jupyter.org/gist/genkuroki...
v1.4.0-rc2.0

を比較してみて下さい。シンプルなモンテカルロで遊んでいる人はこれだけでもメリットがあるかも。

タグ： Julia言語

posted at 00:22:07

まえはる @maeharu

妻の小言。はいつ読んでも面白いよね。最近また更新多くてうれしい😊
inaw.exblog.jp/5958673/

タグ：

posted at 01:19:16

ニケ @LingkoNIKI

蛇が怖い方は見ないでね。

タオルを誤飲しちゃった蛇さん。獣医師さんがずりーーっと引き抜いてます。 twitter.com/cbsnews/status...

タグ：

posted at 01:19:43

#Julia言語

v1.3以降で遅くなったのは、rand()の繰り返し。
rand(10^4) の類は遅くなっていない。

タグ： Julia言語

posted at 01:21:21

Net Zero Watch @NetZeroWatch

〈 Berger | Dillon 〉 @InertialObservr

Freeman Dyson, World-Renowned Physicist & GWPF Founding Member, Has Died Aged 96 - The Global Warming Policy Forum (GWPF) www.thegwpf.com/freeman-dyson-...

タグ：

posted at 04:16:10

eⁱˣ = cos(x) + i sin(x) pic.twitter.com/4fCuF4LOY3

タグ：

posted at 05:24:35

きなこ @3h4m1

〈 Berger | Dillon 〉 @InertialObservr

前倒し春休みにサンシャイン池崎級の歓声をあげて帰宅したウチの息子も今朝

「学校も塾も算盤教室も休みに..寧ろ外に出るなと言われたのだが」
「これから家でする事..妹の子守？」

何かを悟ったらしく

スンッ..

となっていてオイなんだその騙されてマグロ漁船に乗せられた人みたいな目は。

タグ：

posted at 06:00:04

A natural question is: Why is this animation smooth?

The answer is that you can define the so-called fractional derivative, which generalizes nth derivative to a real number of derivatives

The fractional Taylor series smoothly interpolates between integer derivatives pic.twitter.com/757q5yvu7k

タグ：

posted at 06:22:30

Shuhei Kadowaki @kdwkshh

Today I wrote (detailed) update instructions of Juno:
hope this help people solve the Juno update messes !

docs.junolab.org/stable/man/upd...

#JuliaLang #Juno

タグ： JuliaLang Juno

posted at 06:36:59

sandship @s_and_shi_P

julialangやるぞの機運

タグ：

posted at 06:52:21

岸田直樹,MD,MPH,PhD,CEO【 @kiccy7777

ニャンコ先生‼️ twitter.com/teramoto_toyo/...

タグ：

posted at 07:21:22

【新コロの検査・治療を理解する事は重要】
ただ
その前にインフルの現状理解に挑戦を

似て非なる部分もありますが
この理解なく新コロのそれを理解するのは
簡単ではありません
【間違いだらけ！インフルの正しい｢怖がり方｣｢検査もしないし薬ものまない｣という選択肢】toyokeizai.net/articles/-/147...

タグ：

posted at 07:24:16

母数がparameterの訳語として定着した背景なのだが、どうやら、日本語に訳された時点で｢母集団分布の形を決める｣の意味は込められていたようだ。ガーン。
tarotan.hatenablog.com/entry/2016/08/...

タグ：

posted at 08:27:48

英語でも調べたら、このサイトで紹介される限りではparameterは、初出では母集団の特性を表す言葉として導入されたらしい。ガーン。｢（母集団分布に限らず）分布の形状を決める数｣というニュートラルな意味でも使われるようになったのはいつ頃なんだろう？
www.statisticshowto.datasciencecentral.com/what-is-a-para...

タグ：

posted at 08:27:49

このガーンは｢知らなかったので驚いた｣のガーンです

タグ：

posted at 08:32:42

#統計いつもお世話になっているTarotanさんのブログではないですか！ガーン！私も知らなかった。

以前は、「母数」が「母集団分布の代表値」と「確率モデルのパラメータ」(母数はparameterの訳語)の両方の意味で使われていて、悲惨なことになっているという問題に気付いて無かった。 twitter.com/not_identified...

タグ：統計

posted at 08:40:14

#統計どのように悲惨なことになっているか。

数学的モデルと現実の混同は非科学的な考え方の典型例です。

だから、「確率モデルに入っているパラメータ」と「現実の未知の母集団の特製値」を混同することは非科学的。

ところがそういう非科学的な考え方が統計学教育では普通に観察されるのです。

タグ：統計

posted at 08:54:24

#統計 Tarotanさんのブログより

【「分布函數に含まれる媒介變數又はその函數を[原文ママ]母集團の特性を表す常數という意味で，母集團常數又は母數という．】（統計科學研究會編『新編統計數値表 I』 1952, p.7）

分布函数に媒介変数(確率モデルのパラメータ)と母集団の特性値を混同している。

タグ：統計

posted at 08:57:15

#統計孫引き元のTarotanさんのブログは

tarotan.hatenablog.com/entry/2016/08/...

です。

ただし「母集団」の意味を「現実とは無関係のモデル内部の仮想世界における母集団」と解釈すれば「モデルと現実の混同」の問題は無くなります。

タグ：統計

posted at 08:59:57

#統計しかし、現実の統計学教育で使われている教科書を見る限り、「現実の母集団分布」(例: S市の中学3年生の身長の分布)と「数学的モデル内の仮想世界における母集団分布」(例: 正規分布モデルにおける正規分布)を完璧に行っているようには見えません。

タグ：統計

posted at 09:04:54

#統計実際、統計学入門の教科書には

　　母集団分布は正規分布であると仮定する。

とように書いてあり、その母集団分布は統計分析で使うサンプルを取得する分布であり、例題でも現実の母集団から取ったサンプルが正規分布のi.i.d.であるかのように扱われていたりします。

非科学的な雰囲気が満載！

タグ：統計

posted at 09:08:07

@kankichi57301 @kankichi57301

@Ra_koyama @sekibunnteisuu ２段が連続しないって条件がなければフィボナッチです。

タグ：

posted at 09:10:45

@kankichi57301 @Ra_koyama 漸化式が、an＝a(n-1)+a(n-2)　になるとは思ったけど、確かにフィボナッチだ。気づかなかった^^

タグ：

posted at 09:12:25

#統計例えば、

S市の中学3年生全体の身長の分布は正規分布であるはずはないのだが、正規分布モデルで推定を行ってみる。S市の中学3年生全体の身長の分布と正規分布の違いの大きさ(それも不明)によっては推定誤差が無用に大きくなってしまう可能性に注意せよ！

のように書いてあれば「普通」だと～

タグ：統計

posted at 09:14:47

#統計続き～思うのですが、

　S市の中学3年生全体の身長の分布は正規分布であると仮定する。

と言われたら

　はい、その仮定は間違っています。

と思わざるを得ないし、もっと抽象的に

　母集団は正規分布に従っていると仮定する。

と言われても、

　え？そんなん仮定していいの？

と感じる。

タグ：統計

posted at 09:14:51

#統計以上のような意味で、統計学入門の分野は

　　学生に非科学的な考え方を吹き込んでいるので要注意

とみなさざるを得ないと思います。

タグ：統計

posted at 09:16:56

#統計以上のような当然あるべき批判に対して、

　何らかの仮定がないと統計的推測は不可能なので仕方がない

と答えることはひどい詭弁であり、非難に値すると思います。

現実のある側面が採用したモデルを使った分析で適切につかまると期待するという意味での「仮定」については問題にしていない。

タグ：統計

posted at 09:21:38

@kankichi57301 @kankichi57301

@sekibunnteisuu @Ra_koyama え”？消防に３項間の漸化式ぃ、えらい時代になったもんだ＞
www.chugakujuken.com/koushi_blog/ta...
＞この結果だけ（フィボナッチ数列）は覚えるように教えている塾もあるようですね。
消防から暗記数学やめれ！

タグ：

posted at 09:22:35

#統計最近紹介した島谷健一郎さんの本では、ポアソン過程モデルによる分析例について解説した後に

【そもそもこれらの樹種がポアソン過程に従って分布している可能性は皆無である】

と書いてあって爽快な気分にさせてくれます。

より正確な理解を目指す健全な科学的態度がそこにあると思いました。 twitter.com/genkuroki/stat...

タグ：統計

posted at 09:27:52

@kankichi57301 @Ra_koyama こういうのって、あれこれ考えるのが面白いのに、中寿対策で覚えるべきこと、となった瞬間に

つまんねーもの

となってしまいますね。

この問題もフィボナッチを知っているかどうかの知識問題。そもそもこれどんな数を入れても正しいはずだけどね。
www.suguru.jp/Fibonacci/nyuu...

タグ：

posted at 09:49:45

#統計悪い解説の具体例

東京大学教養学部統計学教室編『統計学入門』より

節のタイトル__のみ__で「母集団分布は正規分布である」という仮定が明示されている。

〈例〉では厚さの測定値と正規分布の関係が何も明示されていない。

何が何で近似されていると期待しているかが全て不明瞭。 pic.twitter.com/oVD4JsjdsP

タグ：統計

posted at 09:51:17

#統計

 twitter.com/genkuroki/stat...

タグ：統計

posted at 09:56:05

#統計関連

twitter.com/genkuroki/stat...

タグ：統計

posted at 09:56:53

@kankichi57301 @Ra_koyama 中寿対策→中受対策

タグ：

posted at 09:59:23

@tsatie とりあえずの計算には何の問題もないでしょう。

摩擦がないとみなすことは理想化の一種であり、現実がそうなっているわけではない、という類のことを強調しない物理の人を見たことがないので問題なし。

統計学入門の解説のひどさと物理のまともさを並べるのはナンセンス。続く

タグ：

posted at 10:01:16

@tsatie 統計学も勉強すると楽しいです。

タグ：

posted at 10:03:42

#数楽

 en.wikipedia.org/wiki/Dyson_con...
Dyson予想のウィキペディア twitter.com/thegwpfcom/sta...

タグ：数楽

posted at 10:21:04

Shinji Kono @shinji_kono

「大学では何も教えてくれなかった。必要なことは全部自分で学んだ」と思って卒業していくのは大学教育的には成功。だめなのは授業で要求されたことだけで卒業してしまうこと。

タグ：

posted at 10:25:15

峰　宗太郎 @minesoh

〈 Berger | Dillon 〉 @InertialObservr

有名な反ワクチン、トンデモですね👶
#マシュマロを投げ合おう
 marshmallow-qa.com/messages/b6190... pic.twitter.com/jCMciwRsNi

タグ：マシュマロを投げ合おう

posted at 10:43:57

There are of course other ways to interpolate between terms in Taylor series

タグ：

posted at 11:09:10

math26 @math26

私から何度もお伝えしたいのはこれです
twitter.com/math26/status/...

タグ：

posted at 11:12:01

#統計 population parametersに関する記事

www.statisticshowto.datasciencecentral.com/what-is-a-para...

で紹介されているFisherさんの論文(1922)は

royalsocietypublishing.org/doi/abs/10.109...

からOCR済みのPDFをダウンロードできますね。

文書内検索すると確かにpopulation parametersと書いてあった。 twitter.com/not_identified... pic.twitter.com/agqAGgBHGj

タグ：統計

posted at 11:13:15

お，モジュール毎に最適化レベルを設定できる変更が。
github.com/JuliaLang/juli...

タグ：

posted at 11:13:41

これは読む価値がある！

【苅谷剛彦】教育改革にひそむ「主体性」「平等主義」という名の落とし穴 globe.asahi.com/article/13166629

タグ：

posted at 13:04:45

ううむ、読む価値はあったが、経済への認識が弱過ぎて、有害な言説もたれながしている点が気になりました。(【財政赤字が巨大になり～】と書いてある。)

データに基いた意見だけを述べてくれれば良かったのになと思いました。

「同調圧力」についても「日本では高い」という説は疑わしいのでは？

タグ：

posted at 13:13:25

現在の日本政府の財政は充分に健全であり、教育への投資のために財政赤字を増やすという選択肢は十分考慮に値します。

どんなに「論」を積み重ねても、「金」を出さないとどうにもならない。

タグ：

posted at 13:15:44

米国の有名大学教授のピンカーさんは米国流の「人物を見る」大学入試に激おこ‼️

ダメなやり方は日本に限らずどこでやってもダメだということなんじゃないですかね。

note.com/kodaikusano/n/...

タグ：

posted at 13:25:52

@OokuboTact @metameta007 @takusansu @sekibunnteisuu 高橋誠さん、おバカさんはそろそろ黙った方がよいと思いました。

大学生向けの線形代数学の教科書でスカラーaとベクトルの積を av と書いて統一しているものは多いのですが、そういう教科書を引用して「著者は va と書くことを否定している」と主張し出したら、バカ扱いして問題ない。それと同じ。

タグ：

posted at 13:31:15

@OokuboTact @metameta007 @takusansu @sekibunnteisuu 最初に「読む価値がある！」と書いてしまったが、正直、かなり落胆してしまった。

タグ：

posted at 13:38:38

#超算数メタメタさんの発見した文献を流用しちゃう。
林鶴一編『算術教科書』上巻東京、金港堂、1899年。
林編(1899)は32頁で #掛算を被乗数と乗数の積と定義dl.ndl.go.jp/info:ndljp/pid...。乗数は必ず不名数と。45頁で積=被乗数×乗数=乗数×被乗数とdl.ndl.go.jp/info:ndljp/pid...。名数に関する留保なし。 twitter.com/metameta007/st... pic.twitter.com/WSbNANNCZa

タグ：掛算超算数

posted at 14:29:49

K.B.砂糖 @KB_satou

allocation をへらすだけで 10 倍早くなるんか…… pic.twitter.com/Yqd3Sp7OtC

タグ：

posted at 15:46:05

#超算数暗算について再説。51年指導要領の暗算説は47年指導要領のそれwww.nier.go.jp/guideline/s22e...を引き継いでいた。曰く【暗算とは，計算過程を全く意識することなく，反射的に結果を得ることをいう。したがって，計算過程を意識するような計算は，外見上暗算のように見えても筆算と考えてよい。】と。

タグ：超算数

posted at 16:35:24

#Julia言語 #julialang

exp(it) (0 ≤ t ≤ 2π) のn次までのTaylor展開をnについて連続的に補間してプロット。 fractional derivative に基く補間になっている。

nbviewer.jupyter.org/gist/genkuroki...

twitter.com/InertialObserv... pic.twitter.com/NyTmSDMLlx

タグ： julialang Julia言語

posted at 16:37:15

#Julia言語 #JuliaLang

補間していないバージョン

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/SFH0iM9jLI

タグ： JuliaLang Julia言語

posted at 16:37:15

#Julia言語 #JuliaLang

数式による説明

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/QD3DUfPouO

タグ： JuliaLang Julia言語

posted at 16:37:16

#数楽関連

twitter.com/abel_prize/sta...

タグ：数楽

posted at 16:39:14

K.B.砂糖 @KB_satou

TODO
これはちゃんとソースをつけるべき twitter.com/KB_satou/statu...

タグ：

posted at 16:43:35

#Julia言語この動画では Plots.jl の pyplot() バックエンドを使って動画を作っています。

線の色を連続的に変化させる方法の例としても貴重かも。

Plots.jl での連続的な色の変化の取り扱いについては

docs.juliaplots.org/latest/colors/

nbviewer.jupyter.org/github/JuliaPl...

を参照。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/5m5hhvnpSB

タグ： Julia言語

posted at 16:43:40

#超算数暗算をこう定義すると、紙と鉛筆を使わない計算でも、計算過程を意識するなら、暗算とは呼べなくなります。また、【反射的に結果を得る】ことを重視すると、教師が幾ら計算を指導しているつもりでも、児童は暗記ととらえる傾向が強まるでしょう。

タグ：超算数

posted at 16:47:09

#数楽三角函数のテイラー展開の収束の様子はよくプロットされていると思うのですが、複素平面上にプロットする方がずっと面白いんですね。

あと、n次までの展開のnを連続的に補完した方が動画にしたときに気持ちよくなる。

しかも、1つの自然な補完法として非整数階微分の考え方を使えると。

タグ：数楽

posted at 16:50:09

万博@盲学校マジック @bampaku

「祖母の家の居間で数分に1回『ピー』と音がする。」をトゥギャりました。 togetter.com/li/1475065

タグ：

posted at 16:52:07

データ圧縮などの「圧縮」の反対は？

タグ：

posted at 16:53:36

#超算数以前引用した内藤美城男(1989)によれば、塩野直道は最初期指導要領は暗算を軽視しているとし、それはGHQの方針だと考えていたようです。それが当たっているのかどうか、47年指導要領については「在米資料」も十分調査されていると思うので、既に研究成果があるかもしれません。

タグ：超算数

posted at 16:53:39

#数楽「三角函数のn次までのベキ級数展開のプロット」のような易しい課題(おそらく大学のコンピュータの授業の類でやっている課題)であっても、遊び心があれば結構楽しめるんですね！😊

タグ：数楽

posted at 16:54:15

日本将棋連盟【公式】 @shogi_jsa

#超算数塩野直道は、アメリカ占領当局が暗算軽視だったのは、彼らが暗算が苦手だからだと論じ、その原因を英語における数え方の不規則性に求めます。

タグ：超算数

posted at 16:59:13

詰将棋（3手詰）です。
「まいにち詰将棋」から過去作を厳選してお届けしています。
わかったら「いいね」ボタンをお願いします。
→ヒント、解答はこちら　buff.ly/2Ot0NtA
#詰将棋 #まいにち詰将棋 #3手詰 pic.twitter.com/amomjO10fC

タグ： 3手詰まいにち詰将棋詰将棋

posted at 17:00:26

#超算数内藤美城男(1989:77,8)に塩野の『日本の暗算』(30頁。書誌不明、絶版)が引用【暗算は、まず、数観念に基づいて計算する…数観念が命数法を根底としたものであるから…命数法は国によってちがう。これが、暗算の方法、したがって、また、その難易の度合いが国によって違う原因となっている。】

タグ：超算数

posted at 17:01:04

#超算数【数観念が命数法を根底としたものである】というのはあまり納得がいきません。内藤の引く塩野がtwenty-fiveにthirteenをたすことは不規則性があって二十五に十三をたすより難しいといっています。しかし、数字やドットを思い浮かべれば大した困難ではないです。4が○○○○になるのと同じ。

タグ：超算数

posted at 17:16:15

確率モデルは人間が定義するものであって、その確率モデルを特徴づける数（＝parameter ）は母集団分布（未知の分布）とはなんら関係がないかもしれない。現代的な統計学はそこから出発すると理解しているけど、歴史的にparameterは母集団分布の形を決める定数として使われてきた経緯があるようだ。

タグ：

posted at 17:50:53

@temmusu_n #超算数

 ameblo.jp/metameta7/entr...
「命数法ありて後数あり」

という複数の意味でトンデモっぽいブログ記事を思い出しました。

「命数法ありて後数あり」という説だけではなく、ピダハンも信じているようだし。

「命数法ありて後数あり」には元ネタがあったんですね。 twitter.com/genkuroki/stat...

タグ：超算数

posted at 17:52:30

#統計訂正

❌～を完璧に行っているようには見えません。
⭕️～の区別を完璧に行っているようには見えません。

「の区別」が抜けていた。

実際には、正規分布モデルに基く区間推定を、正規分布に近似的にもしたがっていない可能性が残っている現実の未知の母集団分布に適用する場面は多いはず。 twitter.com/genkuroki/stat...

タグ：統計

posted at 17:57:40

Katsushi Kagaya @katzkagaya

We have updated the 'animal personality in sponge crab' preprint on bioRxiv with a new author @nhayashi1994 contributing to a mathematical description of the Bayesian model selection by information criteria such as (W)AIC, (W)BIC in the discussion:
www.biorxiv.org/content/10.110...

タグ：

posted at 18:54:27

@genkuroki #超算数スレッドはごちゃごちゃしているのでゆっくり読みますが、ブログ記事は読みました。ある数にどんな言語記号を割り当てるかはあまり重要ではない証拠としてシュメール・アッカド語による数学の成果を挙げます。粘土板上の楔形文字による数学のこと。詳細は『数学史』の室井氏担当箇所を参照。

タグ：超算数

posted at 19:21:57

@genkuroki #超算数シュメール・アッカド語による数学では60進法で数を記しましたen.wikipedia.org/wiki/Sexagesimal。これは位取りがうまくいかないことがあるので、現代の記数法には及ばないが極めて有用であったことが今までの研究で明らかに。一方、シュメール語またはアッカド語で文章されていない数の名前がある。

タグ：超算数

posted at 19:25:09

@genkuroki #超算数 it.stlawu.edu/~dmelvill/meso...によると、シュメール語で100をなんと呼んだか、我々はまだ知らないんです。これから出土する粘土板で分かることもあるかもしれないが。アッカド語でも19をなんと言ったか知られていないです(セム語だし他の1xから外挿は可能)。名前は重要じゃないという感じです。

タグ：超算数

posted at 19:29:49

分数べき微分にあたる積分変換の最初は，ODEの教科書にも書かれているEuler変換であろう
books.google.co.jp/books?id=HKQe3...
Abelは積分方程式の祖とよく言われるが，何を持って「積分方程式」と思うかによるだろう。ノルウェーの雑誌掲載時は1823なので，もうすぐ積分方程式200年である
books.google.co.jp/books?id=XiUAA... twitter.com/genkuroki/stat...

タグ：

posted at 20:02:19

Abel全集は何度か紹介しているが，フランス語に翻訳されている2巻本が手軽かつ今の定番だろうと思う。

積分方程式の論文は1巻の11ページ
1: gallica.bnf.fr/ark:/12148/bpt...
2: gallica.bnf.fr/ark:/12148/bpt...
#HDDの肥やし

タグ： HDDの肥やし

posted at 20:09:33

島谷さんは修士は京大理で戸田さんがいちおう指導教授だったはず。博士でRIMSに移ったことは知っている人も多いが、指導教授が佐藤幹夫だったことを覚えている人は少ないかもしれない。その後，統数研にポストを得るまでの人生は，，140字に書くには狭すぎる。 twitter.com/yoshitakeh/sta...

タグ：

posted at 20:15:21

高橋誠さんのこの「数量感」ってなんなんだろうか？
数量感が同じとか違うとか言う以前に、「何それ？」という感じ。
ameblo.jp/metameta7/entr...

タグ：

posted at 20:55:10

ameblo.jp/metameta7/entr...
＞その後，遠山啓や森毅が，ペアノの公理系は過大評価されている，たとえば位取り法を無視していると言っているのを知ったとき，少し安心した。

これ、あやしいと思って直接聞いたことがある。

タグ：

posted at 20:59:53

twitter.com/sekibunnteisuu...

タグ：

posted at 21:01:47

OokuboTact　大久保中二病中年 @OokuboTact

twitter.com/metameta007/st...

この遠山啓の言っているのも「何つまんないこと仰々しく言ってんの？くだらね～」と思ってしまうけど、

それでもメタメタさんが要約しているような話じゃなくて、「ペアノ公理を算数教育に持ち込むことには反対」というだと思う。そりゃそうだろう。

タグ：

posted at 21:08:56

@sekibunnteisuu ＞　それでもメタメタさんが要約しているような話じゃなくて、

本気で誤読している人からすれば、「自分は嘘は言っていない、これが証拠です」となるわけだから、
自分の誤読に気づかなれば平行線状態のままになるわけで・・・ pic.twitter.com/HRJmdUxSEB

タグ：

posted at 21:14:40

@OokuboTact これのどこをどうすれば、【遠山啓や森毅が，ペアノの公理系は過大評価されている，たとえば位取り法を無視していると言っている】とう解釈になるのやら。

タグ：

posted at 21:16:40

#Julia言語

 discourse.julialang.org/t/mcmc-landsca...

タグ： Julia言語

posted at 21:22:18

@OokuboTact しかしこれ、遠山啓もつまんね～こと言っていると思う。
順序数だの数え主義だのを仮想敵としたドンキホーテに思えてしまう。

タグ：

posted at 21:23:32

@OokuboTact pic.twitter.com/aebAOaFvvM

　私も多少は「遠山啓はこういう粉と言っているんだな」というのある程度わかるが、普通の人が読んだら「なんだこれ？」という感じだと思う。

　数教協会員もよく分かってないんじゃないかな？

タグ：

posted at 21:26:45

@OokuboTact ameblo.jp/metameta7/entr...
＞ペアノの公理系は順序数からの記述であり，私の数観が（おそらく多くの日本人と共通する）基数的なものであり，またペアノの公理系の発想の底には一神教的なものがあるから，一神教的なものに馴染めない私が（ほとんど憎悪に近い）違和感を感じたのは無理がなかったのだと思

タグ：

posted at 21:32:49

@OokuboTact 突っ込みどころ満載。メタメタさんは順序数vs基数、という対立構造を勝手に想定して、ペアノ、一神教は順序数、自分や日本人は基数、というようなワケノワカラナイ話に持って行っている。

　ペアノの公理はそんな話じゃないのに。

タグ：

posted at 21:34:53

舞田敏彦 @tmaita77

突然の休校、給食業者に衝撃　「涙止まらない」国に補償要求へ　発注済み2万食どうなる　headlines.yahoo.co.jp/hl?a=20200229-...

タグ：

posted at 22:28:53

OpenBLASと互角(というかちょっと速い?)なマルチスレッドでpure Juliaの行列積!?🤔
discourse.julialang.org/t/ann-loopvect...

タグ：

posted at 22:55:44

にこ( ´ω` ) @nikoyky

我らが上昌広先生が出演されたそうなので、よく知らない方のために

上先生は、PCR検査のビジネス化に非常に熱心（意味深）なお方です。

#ニュースキャスター pic.twitter.com/qRFQjujpPC

タグ：ニュースキャスター

posted at 22:55:51

コード見ると下付き小文字使ってる。
github.com/MasonProtter/G... pic.twitter.com/obJPOkCIxz

タグ：

posted at 22:58:11

Juliaのコードはやっぱりキュートやな。うまく説明できんが。

タグ：

posted at 23:00:27

gzipくらいJuliaで実装してみるか(誰得)

タグ：

posted at 23:06:32

「積分方程式」も日本語の良い本が少ないなあと書こうと思ったが、先週に上村豊さんにお会いしてたことを思い出したのでセーフ：
「積分方程式―逆問題の視点から―」
www.kyoritsu-pub.co.jp/bookdetail/978...
上村さんのこの本は専門書であるが，第1章は「Abel積分方程式とその遺産」である。

タグ：

posted at 23:12:20

#統計階層 WAIC を検索して発見しました。

仮にMCMCで近似的に解いているなら、階層ベイズでWAICをどのように計算しているのだろうか？既存のパッケージでWAICやLOOCVを計算すると、正しい解釈でそれらを正しく計算してくれない可能性が高い。

解説に続く twitter.com/wnk_datasci/st...

タグ：統計

posted at 23:12:49

タグ：統計

posted at 23:22:11

残り12件のツイートを見る

#統計続き～、サンプルY_kを生成した未知の分布の推測において、

(B) Y_k ～ p(y_k|z)
　　z ～ φ(z)

と同値なモデルになります。内部パラメータ w_k が積分によって消去されている。

以上の(A)と(B)の違いがWAICやLOOCVを求めるときに問題になります。続く

タグ：統計

posted at 23:22:14

しみずとしやす @watermanWork

今日、10歳娘が粘土をこねて一所懸命作っていたものは、なんと、「耳」です。
耳…粘土で作ろうとか思ったこと1回もないなぁ…。 pic.twitter.com/zIObuGq2Dd

タグ：

posted at 23:29:36

#統計 (A)のモデルをMCMCで解くと、内部パラメータも含むパラメータ (w_1,…,w_n,z) の事後分布のサンプル(chainとよく呼ばれる)が得られます。

その事後分布のサンプルには少なくとも2通りの解釈がありえます。続く

タグ：統計

posted at 23:32:10

#統計 ① モデル(A)中の(Y_1,…,Y_n)を、確率モデル

p(y_1,…,y_n|w_1,…,w_n,z)=Π_k(p₁(y_k|w_k)p₂(w_k|z))

のサイズ1のサンプルとみなす解釈。この解釈はモデル(A)とモデル(B)と同値だとみなさない立場に立つことになり、サンプルサイズも1になります。続く

タグ：統計

posted at 23:32:16

タグ：統計

posted at 23:36:49

Brian Roemmele @BrianRoemmele

It was 1867 and James Clerk Maxwell suggested how the second law of thermodynamics might hypothetically be violated.

In the thought experiment he imagined a demon decreasing entropy and violating the second law of thermodynamics.

Meet a Maxwell demon:

pic.twitter.com/msKR3lVDmx

タグ：

posted at 23:38:23

#統計階層ベイズでWAICを求めて利用している人は検索しても(松浦さんや清水さん以外には)なかなか見つかりません。

その理由は②の解釈で数値積分が必要になるからだと思われます。

計算のノウハウを公開すれば助かる人は多いと思ったので、このスレッドを書くことにしました。

タグ：統計

posted at 23:39:59

#統計階層ベイズであっても、事後分布の形が綺麗に単峰型になっているなら、自由エネルギーの方であればブリッジサンプリング(複数の変種がある)が使えるはず。

自由エネルギーの方であれば(パラメータを全部積分して対数を取って-1倍したものなので)2通りの解釈のどちらでも同じ値になります。

タグ：統計

posted at 23:43:35

#統計 WBICの方は逆温度 1/log n のMCMCを実行して得た z に関するサンプルと p(y|z) = ∫p₁(y|w)p₂(w|z)dw を使えば近似計算できます。しかし、WAICと同様に一般には p(y|z) の計算に数値積分が必要になります。

その代わり、事後分布がシンプルな単峰型にならないケースでも使える。

タグ：統計

posted at 23:47:05

K.B.砂糖 @KB_satou

Julia は反復動作で速度があがるのできめつ

タグ：

posted at 23:50:10

Arturo Erdély @ArturoErdely

The #Distributions package in #julialang provides a large collection of probabilistic distributions and related functions. juliastats.org/Distributions....

タグ： Distributions julialang

posted at 23:52:22