Twitter APIの仕様変更のため、「いいね」の新規取得を終了いたしました

黒木玄 Gen Kuroki

@genkuroki

  • いいね数 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 を見て下さい。
並び順 : 新→古 | 古→新

2019年09月07日(土)

三村智保 @igomimu

19年9月7日

挑戦!創作詰碁【大西智之34☆☆☆】黒先結果如何 bit.ly/2KBmvdK
#創作詰碁

タグ: 創作詰碁

posted at 00:08:29

小口貴宏 / Takahiro Kogu @TKoguchi787

19年9月7日

Galaxy Fold、控えめに言って感動しますね。HUAWEIの外折りより感動したかも pic.twitter.com/H2C0Olw3Wp

タグ:

posted at 00:09:41

(「・ω・)「ガオー @bicycle1885

19年9月7日

Juliaするか,超Juliaするかの2択で考えてください。

タグ:

posted at 03:10:26

Space Explorer Mike @MichaelGalanin

19年9月7日

Jupiter viewed from its South Pole.
Credit: NASA/JPL pic.twitter.com/9ZVCplssa2

タグ:

posted at 04:26:58

SAKAUE Akkuş Tatsuya @sakaue

19年9月7日

[stats] ベイズ統計学の概論的紹介 buff.ly/2HyEbo1

タグ:

posted at 05:30:08

Simon Byrne @SimonPJByrne

19年9月7日

Most of the time #JuliaLang helps solve the two language problem. Sometimes, however, it makes it worse: github.com/JuliaInterop/M...

タグ: JuliaLang

posted at 06:12:30

K.B.砂糖 @KB_satou

19年9月7日

N重振り子3D.jl (N=10)
ルンルン♪ ルンゲ=クッタ pic.twitter.com/9LSFGoHZZb

タグ:

posted at 06:43:57

山田太郎 ⋈(参議院議員・全国比例) @yamadataro43

19年9月7日

35人が亡くなり多くの負傷者も出た京都アニメーション放火事件。25億円を超える寄付金への税制優遇を決めました。当初から粘り強く働き掛けてきた結果です。これが京アニ事件だけでなく放火事件で苦しむ他の事案についても税制優遇の道を切り開く前例になればと考えています。www3.nhk.or.jp/news/html/2019...

タグ:

posted at 07:36:30

こなみひでお @konamih

19年9月7日

今日も日本システム企画の熊野社長がカンカンになるようなネタを投下しますので,みなさんお楽しみに。火に油を注ぎ返してやるのだ。

タグ:

posted at 07:58:36

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 指数型分布族にのみ適用できる方法を見たら、「特別に簡単な場合にのみ適用できる方法であり、より複雑なモデルには適用できない」と判断するべき。

指数型分布族はかなり一般的な確率分布族であるというのはひどい誤解。

twitter.com/genkuroki/stat...

タグ: 統計

posted at 09:47:09

たまさか @TamasakaTomozo

19年9月7日

娘さんから、自分に孫ができることをサプライズで知らされた男性のリアクションが、反則級に可愛い。
pic.twitter.com/etoTXxWXlZ

タグ:

posted at 09:56:14

ルーン@( 。・ω・)⊃━☆゜*・。 @rune_worth

19年9月7日

いずれ来るとは思っていたけどやはり来たか。ウチにはまだ広告来てないけど。結局最後の稼ぎ場所は医療系に落ち着くんだよなぁ…。
NMRパイプテクターの原理を応用した謎の「血液還元装置」 www.gohongi-beauty.jp/blog/28817/ @kuwamitsuosamuさんから

タグ:

posted at 10:15:23

ヤマンザ子 @ktykgucci

19年9月7日

悲しみの実戦詰将棋です(44金打てば悪いなりにもう少し長引くと思った) pic.twitter.com/yGq5dwMQwr

タグ:

posted at 10:42:56

かず@猫の下僕 @ns_pooh_

19年9月7日

コバ先生にhpvワクチン絶対に娘にも自分にも打たないとリツイートしていた女性。

まだお子さんは5歳のようですので、まだ猶予があります。あと7〜10年くらい、、、
それまでに厚労省が推奨に舵きりをしていただけたらいいなと思います。

タグ:

posted at 11:12:51

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#Julia言語 から #Wolfram言語 をMathLinkライブラリ経由で使う方法

まだ表示が数式に整形されないのですが、JuliaからWolfram言語を使う道が開かれたので、今後の発展が期待されると思います。誰かが表示の整形をなんとかしてくれるだけで、相当に実用的になると思う。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/4meE1jlXvh

タグ: Julia言語 Wolfram言語

posted at 11:15:19

飯間浩明 @IIMA_Hiroaki

19年9月7日

「『十回』の『十』は『ジュッ』でなく『ジッ』と読む、本来の字音『ジフ』の変化だからだ」という主張があります。歴史的にはそうですが、現在の常用漢字表の「十」には「ジッ」の備考欄に「『ジュッ』とも」と記されています。現在一般的な「ジュッ」がお墨付きを得た格好。テストでは両方マルです。

タグ:

posted at 11:48:39

非公開

タグ:

posted at xx:xx:xx

Andrew Whitby @EconAndrew

19年9月7日

This is top shelf trolling, because thanks to Excel “1 in 5” genetics papers contain errors in gene names. www.sciencemag.org/news/2016/08/o... twitter.com/msexcel/status...

タグ:

posted at 12:35:00

K.B.砂糖 @KB_satou

19年9月7日

N重振り子3D.jl (調子に乗ってN=100)
ちゃんと発散せずに動く。ルンゲ=クッタならね。 pic.twitter.com/MUPmO4hf2E

タグ:

posted at 13:09:38

K.B.砂糖 @KB_satou

19年9月7日

N=500, 紐っぽい pic.twitter.com/tMX8h5UBR4

タグ:

posted at 13:24:12

Lukas Schlögl @LukasSchlogl

19年9月7日

Just priceless. An estimated 20% of genetic research papers contain errors because Excel converted some gene names into calendar dates. twitter.com/EconAndrew/sta... pic.twitter.com/ts6sZbFa1C

タグ:

posted at 13:37:45

ごまふあざらし(GomahuAzaras @MathSorcerer

19年9月7日

漢字で書いてある日本語の本をひながなにして
大人に渡すとナメてんのかって思われるのに,
数学書を数式を使わずに書いた本を書店に置くと
バカ売れする現象に名前をつけたい.

タグ:

posted at 13:54:59

K.B.砂糖 @KB_satou

19年9月7日

重力も入ってます、github にソースを置いたのでご確認いただければ。
twitter.com/ktonegaw/statu...

タグ:

posted at 13:59:51

黒木玄 Gen Kuroki @genkuroki

19年9月7日

1つ前のツイートで紹介したボールチェーンもしくはニュートンビーズを数値計算で再現してくれないものかしら。

twitter.com/kb_satou/statu...

タグ:

posted at 14:00:17

K.B.砂糖 @KB_satou

19年9月7日

計算に対して動画出力が長すぎる……
だが record 以外の動画出力方法がわからん。 pic.twitter.com/caBnyPovfh

タグ:

posted at 14:03:22

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#Julia言語 #Makie.jl

ソースコードの yossha7.mp4 にうけた!

twitter.com/kb_satou/statu...

タグ: Julia言語 Makie

posted at 14:11:20

黒木玄 Gen Kuroki @genkuroki

19年9月7日

ボールチェーンがカップから飛び出る様子を再現するためには、カップの底面に接している質点には重力が働かないようにして、初期条件を適切に(←これが問題)すればいいのかな?

手を動かしていないのでまだ感覚がつかめていない。

タグ:

posted at 14:14:16

広島ベイズ塾 @bayesjuku

19年9月7日

【塾生の活動記録】Stanでポン・デ・リングを作ってみた bayesmax.sblo.jp/article/181987...

タグ:

posted at 14:31:15

naoki10007 @naoki10007

19年9月7日

@konamih なんでしょう 楽しみですね

タグ:

posted at 14:45:51

K.B.砂糖 @KB_satou

19年9月7日

はい、みなさんが静かになるために空気抵抗を設定しました.jl pic.twitter.com/oOol5TdRVB

タグ:

posted at 14:52:13

natsuo @natsuokawai

19年9月7日

駆け足でしたがベイズ学習に入門した話をしました。
発表に使った notebook はこちらです。#StudyAIchi

github.com/natsuokw/machi...

タグ: StudyAIchi

posted at 15:20:06

栗原裕一郎 @y_kurihara

19年9月7日

文芸誌が「お金」特集するっていって呼んで来たのが大澤真幸とか進歩しないねえ。

タグ:

posted at 15:44:36

S (ツイートはスレッド全体をご確認く @esumii

19年9月7日

@shiba_masa 国大協は国立大学の「総意」として明確に反対意見を述べましたが、その発言を座長が遮ってまで強引に推し進めた、というのが事実です。審議会の議事録をご確認ください。web.archive.org/web/2015022415...

タグ:

posted at 15:56:09

S (ツイートはスレッド全体をご確認く @esumii

19年9月7日

@shiba_masa なお、現行の入試が「暗記中心」というのも、少なくとも多くの国立大学については(重大な)事実誤認です。

タグ:

posted at 15:57:30

kagshun/EMANON @kagshuntravel

19年9月7日

カバ。ライオンの5倍、サメの50倍人間を殺している獰猛な生き物ですが

このスピードを見てください。
そして何より驚くのは、泳いでるのではなく、水中の底を走っているという事実。 twitter.com/thescarynature...

タグ:

posted at 16:27:07

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#Julia言語 のDynamicHMC.jlに混合正規分布

p(x|a,b)=(1-a)N(x)+aN(x-b), N(x)は標準正規分布

で作った尤度と事前分布φ(a,b)=exp(-b^2/20)/√(20π)の積の対数を与えてMCMCさせてみたのだが、尤度部分を

log((1-a)+a*exp(b*x-b^2/2))

に変えたら、サンプルサイズn=2^10で3割程度速くなった。 #統計

タグ: Julia言語 統計

posted at 17:16:48

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#Julia言語 #統計 a,bに関する定数函数を除いて、よく呼ばれる対数尤度函数を簡素化した。

DynamicHMC.jl は対数尤度函数+対数事前分布を与えればNUTSでMCMCしてくれるだけの素朴な作りになっているので、こういうことを試すには都合がよいです。

github.com/tpapp/DynamicH...

タグ: Julia言語 統計

posted at 17:16:51

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#Julia言語 #統計 さらに試しに、MCMCで事後分布のサンプルを作らずに、数値積分パッケージ

github.com/stevengj/HCuba...

を使って、(a,b)に関する積分も実行してみましたが、この場合には2変数の積分で済んでいるので結構うまく行きます。

この情報は自分で実験したい人には参考になるかも。

タグ: Julia言語 統計

posted at 17:16:52

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #Julia言語 要するに、多変数函数の数値積分を計算したいだけなのですが、それが結構大変だという話。

タグ: Julia言語 統計

posted at 17:18:34

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #Julia言語 ベイズ統計の場合には被積分函数(本質的に尤度函数)をコンピューター上で適切に作ること自体に工夫が必要です。

logsumexpという函数の実装で注意するべきことを知っている人にはお馴染みの問題が生じる。

en.wikipedia.org/wiki/LogSumExp

タグ: Julia言語 統計

posted at 17:21:51

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 確率論がらみの数値計算では「沢山の正の実数を掛け合わせたもの」を計算したくなることがよくあります。尤度函数

(a,b)↦p(X₁|a,b)…p(X_n|a,b)

はまさにそういう例になっている。しかし、これはコンピューター上で標準的に使われている浮動小数点数とは極めて相性が悪い。

タグ: 統計

posted at 17:26:18

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #Julia言語 Float64では 10^{-383}~10^{384} 程度の正の浮動小数点数を扱えます。log_{10} を取ると -384~384 と非常に狭い!

沢山の正の実数をかけるとこの範囲を容易に飛び出して、浮動小数点数側では0や∞になってしまいます。

これを防ぐには、正の実数の対数を足し合わせればよい。

タグ: Julia言語 統計

posted at 17:32:07

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 正の実数値函数f_i(w)の積

f_1(w)…f_n(w)

ではなく、対数和

log f_1(w) + … + log f_n(w)

なら浮動小数点数の範囲で安全に計算できる。

しかし、f_1(w)…f_n(w)が被積分函数の場合には対数和の exp を取って元の積に戻さなければいけない。そしてexpを取った瞬間に0や∞になったりする。

タグ: 統計

posted at 17:36:43

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 そういう場合には、積分領域K内での対数和の最大値をMとして

g(w) = exp(log f_1(w)+…+log f_n(w) - M)

とおけば、g(w)は真っ当に使える函数になりやすい。そしてf_i(w)の積の積分を

Z = ∫_K g(w) dw + M

で数値計算。これの対数log Zが重要な量になっているというパターンがよくある。

タグ: 統計

posted at 17:41:24

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 欲しい計算が

expを取る → 和もしくは積分 → logを取る

のパターンになっていることはよくある。logsumexpのパターン。和ではなく積分版も結構出て来る。

タグ: 統計

posted at 17:44:12

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #Julia言語 ただし、logsumexpのくそ真面目な実装では最後のlogを取るときのunderflowの問題も気にする必要がある。詳しくは以下のリンク先のスレッドを参照。

Juliaではlog(1+a)の代わりにlog1p(a)を使えばよい。

log(1+a)のaの函数としての直接的実装を使う。

twitter.com/genkuroki/stat...

タグ: Julia言語 統計

posted at 17:50:24

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 以上に書いたようなことは、

p(x|a,b)=(1-a)N(x)+aN(x-b), N(x)=exp(-x²/2)/√(2π)

の尤度函数をサイズ2000程度のサンプルで作って、プロットしてみれば実感できると思います。うまく計算できれば、正則ではない確率モデルの尤度函数の様子を観察できる。

複数の意味で非常に良い演習問題。

タグ: 統計

posted at 17:56:26

おばけ @triwave33

19年9月7日

コメントいただいてPython (Numba) の方が3倍速いよって言われてマジか!ってなってる。なるべく早く検証する

qiita.com/triwave33/item...

タグ:

posted at 18:24:14

S (ツイートはスレッド全体をご確認く @esumii

19年9月7日

中身がないのにプレゼンテーション能力だけ高い、は最悪のパターン。

タグ:

posted at 18:53:03

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#数楽 Laplace近似の弱バージョンもlogsumexp関係の話。

log Σ_i exp(n x_i) = n max{x_i} + o(n) as n→∞.

Σ は ∫ でもよい。

大偏差原理もLaplace近似の弱バージョンなので、大偏差原理もlogsumexp関係の話。

logsumexpパターンは普遍的。

タグ: 数楽

posted at 19:15:52

Norio Nakatsuji @norionakatsuji

19年9月7日

グーグルが、医学的根拠がない「幹細胞治療」の宣伝を禁止するらしい。日本の新聞や雑誌やTVやネットメディアも見習って欲しい。「認可なき/根拠なき治療」は、病気で苦しみ藁をもすがる気持ちの患者を騙して金儲けをするだけでなく、患者の命を危険にさらす、社会的だけでなく人間として最低の罪悪。 twitter.com/washingtonpost...

タグ:

posted at 21:16:31

狸穴猫/松村りか @mamiananeko

19年9月7日

@AS_Insects @genkuroki リンク先のスライドの61~65も興味深いですよ~。
mutsumino.info/ppt/Dr.N_world...

タグ:

posted at 22:22:38

こなみひでお @konamih

19年9月7日

記事を書いたのですがあまりにも先方にとって刺激的なものになってしまっているので,サーバーを止めろという仮処分を出されると理科探みたいなことになりかねません。公開の仕方を考慮中です。著作権侵害やデマとは無縁なのですけど,彼らは見境がないですからね。どうするかなあ。 twitter.com/konamih/status...

タグ:

posted at 22:45:41

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #信頼区間 最近の事例

www.dr-kid.net/pand95ci
疫学からみたP値と95%信頼区間について
2019-08-31

* p値の解説は非常に良い。

* 95%信頼区間の説明も良いが、「使用した統計モデルが正しい」という仮定との関係で、現実には95%という数値の信頼性の問題があることに触れていない。続く

タグ: 信頼区間 統計

posted at 23:03:53

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #信頼区間

www.dr-kid.net/pand95ci
疫学からみたP値と95%信頼区間について
2019-08-31

の非常によろしくない極めて不快な点は内容ではなく、マウスでコピーして引用できないようにしてあること。

批判するための引用に開かれていることはこの手のサイトでは絶対に必要なことだと思う。

タグ: 信頼区間 統計

posted at 23:03:58

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #信頼区間 続き

p値の解説では、現実の母集団に(近似的にも)適用できない統計モデルを使った場合には(具体的には例えば正規母集団の仮定が不適切な場合には)、p値は信頼できない数値になることを強調してある点が非常に良いです。

タグ: 信頼区間 統計

posted at 23:03:58

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 #信頼区間 続き

信頼区間についても、95%信頼区間を求めるために使った統計モデルの現実の母集団への適用が不適切な場合には、95%という数値が信頼できないものになってしまうことも強調した方が良かったと思う。

タグ: 信頼区間 統計

posted at 23:03:59

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 個人的な意見では、仮説検定での帰無仮説の中には用いた統計モデルの仮定も入れてしまった方が誤解が少ないと思う。例えば「母集団分布の平均はμである」ではなく、「母集団分布は正規分布でかつその平均はμである」を帰無仮説にしてしまう。続く

タグ: 統計

posted at 23:13:43

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 「母集団分布は正規分布でかつその平均はμである」という帰無仮説が棄却される場合には

(1) 母集団分布に正規分布モデルを適用することが妥当。
(2) 母集団分布の平均はμである。

の少なくともどちらかが棄却されることになる。続く

タグ: 統計

posted at 23:16:56

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 続き。もしも本当に棄却したい仮説が

(2) 母集団分布の平均はμである。

であるならば、別の議論で

(1) 母集団分布に正規分布モデルを適用することが妥当。

が正しいことを示さなければいけない。(母集団分布がほぼ正規分布になっていることまで示す必要はない。それは普通は無理)

タグ: 統計

posted at 23:21:05

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 正規分布モデルを使って求めた平均の95%信頼区間の定義は「正規分布モデルM(μ)が有意水準5%で棄却されないμの範囲」なので、仮説検定の場合と同じ注意が必要です。

タグ: 統計

posted at 23:27:35

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#統計 別の議論で正規分布モデルの使用が妥当であることを示すことができない場合には

【100回研究を繰り返した場合、95回以上は真の値~を含んでいる区間が95%信頼区間といえます】

という主張は誤りになります。

おそらく、論文査読者の多くはこの事実を無視しているのではないか?大丈夫? pic.twitter.com/9ADBSzwKgc

タグ: 統計

posted at 23:27:40

James Tanton @jamestanton

19年9月7日

How many essentially distinct ways are there to draw five circles on a page? pic.twitter.com/b6jmH5cIpW

タグ:

posted at 23:30:55

黒木玄 Gen Kuroki @genkuroki

19年9月7日

#数楽 一般に「権威的に響く」もしくは「魅力的に感じられる」可能性のある用語に引っ張られる「理解」を目指すとろくなことがないと思う。

「Kullback-Leibler」「情報量」「大偏差原理」などの用語自体に理解のための重要性はない。

どれも簡単で素朴で一貫した普遍的な考え方で理解できる。

タグ: 数楽

posted at 23:40:11

鈴木みそ @MisoSuzuki

19年9月7日

安いエアを乗り継いで行くので荷物重量制限が厳しい。と聞いたので嫁が急遽作った内ポケットカーディガン。パンパンにポケットに詰めて3キロ太って飛行機に乗り込んだらしい。万引き犯か。 pic.twitter.com/3agfinsXI8

タグ:

posted at 23:40:25

@genkurokiホーム
スポンサーリンク
▲ページの先頭に戻る
ツイート  タグ  ユーザー

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

掛算 統計 超算数 Julia言語 数楽 JuliaLang 十分 と教 モルグリコ 掛け算

※タグはスペースで区切ってください

送信中

送信に失敗しました

タグを編集しました