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 を見て下さい。
Favolog ホーム » @genkuroki » 2017年09月30日
並び順 : 新→古 | 古→新

2017年09月30日(土)

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang 定義域を外れたときに無限大の値を返すようにしおかないと、optimize函数が定義域を外れたところを調べようとしたときにエラーが出て止まってしまう。この問題のせいで、変更前には初期条件をうまく選ばないと最尤法が実行不可能なコードになっていた。😭

タグ: JuliaLang 数楽

posted at 22:50:36

@kuri_kurita

17年9月30日

下がるよりいいでしょ。😐
まだ足りない、もっと上げろ、と言えばいいのであって(実際、今の傾向が続けば確実にそうなる)、「これだからアベノミクスは失敗、やめるべき」とか言いたがる人がいるのはナゼなのか。 twitter.com/mugu_ri/status...

タグ:

posted at 22:48:29

積分定数 @sekibunnteisuu

17年9月30日

@kankichi573 @AS_Insects #超算数
ちなみに、中学数学では単項演算としての「-」を扱わないようです。
-(ー3)=3 というような式が見た範囲では教科書の殿にも書いてありません。

ーx にx=-7を代入する場合は

ーx
=(ー1)×x
=(ー1)×(ー7)
=7

としています。学校図書1年教科書より

タグ: 超算数

posted at 22:36:10

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang

nbviewer.jupyter.org/gist/genkuroki...
測定スキルの異なる7人の科学者たち

を更新。optimize函数で最小化する函数が定義域を外れたとき値無限大を返すように変更した。最尤法をoptimize函数を使ってやっている。

タグ: JuliaLang 数楽

posted at 21:25:45

Haruhiko Okumura @h_okumura

17年9月30日

美文書第7版のTeX(TeX Live 2016)を入れたMacをHigh Sierraにしてしまってヒラギノが使えなくなって困っているかたのためのダーティハックをとりあえず書きました okumuralab.org/bibun7/ (あわてて書いたのでやや心配)

タグ:

posted at 20:54:45

非公開

タグ:

posted at xx:xx:xx

ウェザーニュース @wni_jp

17年9月30日

新潟県上越市で発生した、竜巻のたまご「ろうと雲」を捉えた動画がこちら。段々と消えていく様子が、間近の視点から撮影されています。
weathernews.jp/s/topics/20170... pic.twitter.com/x2JYJLTzBe

タグ:

posted at 20:24:05

黒木玄 Gen Kuroki @genkuroki

17年9月30日

もしかしたら、

gist.github.com/genkuroki/bc35...
測定スキルの異なる7人の科学者たち

#JuliaLang によるグラフのプロットの仕方の良いexampleになっているかも。PyPlotを使っているので #Python 初心者にも参考になるかもしれません。

タグ: JuliaLang Python

posted at 20:22:54

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang 登場して来る逆ガンマ分布は、0<α<1 の場合なので期待値は無限大になっています。分散の期待値が無限大(笑)。

en.wikipedia.org/wiki/Inverse-g...

タグ: JuliaLang 数楽

posted at 20:16:36

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang 添付画像はそれぞれのモデルにおける標準偏差の予測分布の累積分布函数のグラフです。よく見ると各モデルが科学者達をどう思っているかの「気持ち」が分かって来る(笑)。 pic.twitter.com/97SFt2eiNM

タグ: JuliaLang 数楽

posted at 20:11:55

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang TDist2さんは、Cauchy分布モデルさんよりも、標準偏差が0.25以下の科学者がずっとたくさんいると思っていますが、それと同時に標準偏差が30以上の科学者もずっとたくさんいると思っています。

タグ: JuliaLang 数楽

posted at 20:09:37

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang TDist1さんは、標準偏差1以下の測定誤差しか出さない科学者は3割もいないと思っていますが、準偏差30以上の測定誤差を出す科学者はTDist2さんより少し多いと思っています。

タグ: JuliaLang 数楽

posted at 20:07:57

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang 更新で追加されたのは1.6.2節です。

正規分布モデルさんは、すべての科学者は標準偏差が13もの誤差を出すと推定してしまっています。正規分布モデルさんは外れ値の存在をそのように解釈している。

タグ: JuliaLang 数楽

posted at 20:05:28

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #JuliaLang

nbviewer.jupyter.org/gist/genkuroki...
測定スキルの異なる7人の科学者たち

を再更新。t分布を使ったモデルを正規分布と逆ガンマ分布の階層モデルとみなしたときの、正規分布の分散の予測分布をプロットしました。各モデルの「気持ち」がわかります。

タグ: JuliaLang 数楽

posted at 20:03:31

ano_ano @ano_ano_ano

17年9月30日

辻元清美「小池氏は寛容な保守の立場、私は現実的なリベラルの立場で、安倍晋三政権を右と左から挟み撃ちに…」

あ~あ辻元、言っちゃった。右と左に挟まれた安倍政権が中道だって認めちゃった。

『民進・辻元清美氏、希望の党に公認申請せず』
www.sankei.com/politics/news/...

タグ:

posted at 19:15:03

でえもん @GreatDemon1701

17年9月30日

#算数 #超算数 たまたま見かけたフランス語での数のまとめ。自然数は{1,2,3,…}または{0,1,2,3,…}のどちらかと説明しています。 pic.twitter.com/RqXvVwGhiC

タグ: 算数 超算数

posted at 18:39:22

内田良:新刊『だれが校則を決めるのか』  @RyoUchida_RIRIS

17年9月30日

学校内での生徒の暴力。
対応できない場合は,さっさと警察等に相談すべきという旨の啓発資料。国がそう言ってるんだ。学校が丸抱えする必要なんてない。

▼文部科学省国立教育政策研究所:生徒指導リーフ「学校と警察等との連携」
www.nier.go.jp/shido/leaf/lea... pic.twitter.com/15Bp2ATTej

タグ:

posted at 18:20:41

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

「間違った説の流布に手を貸すのはやめて欲しい」という話をしているのに「文科省にきいてください」という反応が来るのは毎度のこと。 #掛算

タグ: 掛算

posted at 17:56:21

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan 【文科省にでも問い合わせたら】
それは #掛算 の順序強制反対と言っている人の一部がやりました。
そして「文科省は順序アリとは言ってない」という結論が出ているのです。私がAngel_takenyanさんに申しあげたいのは「間違った説に加担するのはやめてください」という事です。

タグ: 掛算

posted at 17:55:06

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan 【順序はずっとどうでもいいことなんてないと思いますけど】
いいえ。中学以降の勉強では小学校での #掛算 の順序のルールはむしろ通用しません。そして、小学校の範囲内でもけっこう違反しているのは、先に示した通りです。

タグ: 掛算

posted at 17:51:59

A級3班国民 @kankichi573

17年9月30日

@sekibunnteisuu @AS_Insects #超算数
なんか発火点の低そうな可燃物をばらまいてるような。そういうのは言及せずにあいまいなままほっとくほうがいいのでは。-3の-は符号なのか演算子なのか-aならどやねん、とか。

タグ: 超算数

posted at 17:51:07

☆大天使たけにゃん☆ @Angel_takenyan

17年9月30日

@golgo_sardine 順序はずっとどうでもいいことなんてないと思いますけど文科省にでも問い合わせたらどうですか? 自分はいつから気にしなくなったのかわかんないですし。

タグ:

posted at 17:49:07

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan 【式が増えてくれば当然掛け算の順番なんて答えに直結してくるところではないのでそこが採点の基準にならないから順番は気にしなくてもいいんじゃないかと】
そこでですが、何年生から順序どうでも良いのですか?明文化されているのではないですか? とお聞きしているのですが。#掛算

タグ: 掛算

posted at 17:41:20

ryugo hayano @hayano

17年9月30日

「一家に1枚元素周期表」stw.mext.go.jp/series.html は,写真入り.どんな元素がどこで使われているか解説しています.無料ダウンロードできます. twitter.com/ogugeo/status/... pic.twitter.com/otxwb67gpL

タグ:

posted at 17:23:19

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 「外れ値を出した科学者」というストーリーがつぼにはまりまくり。😆

タグ: 数楽

posted at 17:00:34

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 「パラメーターの推定」というパラメーター思考にとらわれると交差検証のような基礎的なアイデアは決して出て来ないと思う。「未知の確率分布自体の推定」という予測分布自体の精度を問題にする発想にシフトして行くべきだと思います。

タグ: 数楽 統計

posted at 16:58:40

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 測定スキルに個人差があるという事実をモデルに組み込みたければ、自然に階層モデルを考えることになります。測定スキルの分布の仕方にある仮定をおくと、Cauchy分布やt分布が自然に出て来ます。その辺のことも公開したファイルに書いてあります。

タグ: JuliaLang 数楽 統計

posted at 16:52:49

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 正規分布をCauchy分布やt分布に置き換えることの「意味」はCauchy分布やt分布がある種の階層モデルに同値であることを見ればわかります。

ある科学者は特別に実験がへたくそでその人の測定結果は全く信用できないかもしれない。続く

タグ: JuliaLang 数楽 統計

posted at 16:51:22

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang

詳細は次のリンク先のソースコードおよび数式入りの解説を見ればわかると思います。ちょこまか無断で更新しているので注意して下さい。

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

タグ: JuliaLang 数楽 統計

posted at 16:43:52

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang

t分布モデル2(パラメーターの個数は3)
これの成績が一番よいのですが、やはり、外れ値での挙動が安定しています。 pic.twitter.com/uuztfeRbpA

タグ: JuliaLang 数楽 統計

posted at 16:41:36

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang

t分布モデル1(Cauchy分布モデルとは異なる方法で2個のパラメーターを入れた場合)
やはり、外れ値での挙動も安定しています。 pic.twitter.com/eRhgibBijk

タグ: JuliaLang 数楽 統計

posted at 16:40:36

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang それでは他のモデルではどうなっているのでしょうか?

Cauchy分布モデルの場合

上段左から二番目の外れ値の場合も含めて一個抜き予測分布(の確率密度函数の対数の-2倍)の形はほぼ同じになっています。Cauchy分布モデルは外れ値に強い。 pic.twitter.com/p7WYp8rBqM

タグ: JuliaLang 数楽 統計

posted at 16:39:11

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 上段左から2番目のグラフは、外れ値を除いて作った予測分布で外れ値がどれだけ予測できるかを表わしています。その場合には損失は200を軽く超えており、これがLOOCVの数値が跳ね上がる原因になっています。正規分布モデルは外れ値にものすごく弱い。

タグ: JuliaLang 数楽 統計

posted at 16:37:15

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 続き。そのことは、正規分布モデルの場合には、外れ値を外して作った予測分布は当然外れ値以外を除いた場合および何も除かなかった場合の予測分布(左上端)とは全然違ったものになることを表わしています。続く

タグ: JuliaLang 数楽 統計

posted at 16:35:25

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 正規分布モデルの場合

左上の黒線はサンプル全体から作った予測分布の対数の-2倍です。丸印はサンプルです。それ以外は一個抜き予測分布の対数の-2倍です。

上段左から二番目だけ、縦軸の縮尺のオーダーが違っており、お椀の底の位置もずれています。 pic.twitter.com/5ar4rheTEs

タグ: JuliaLang 数楽 統計

posted at 16:33:29

☆大天使たけにゃん☆ @Angel_takenyan

17年9月30日

@golgo_sardine 単純に求めるからこそ順番が大事になるんじゃないですか? 式が増えてくれば当然掛け算の順番なんて答えに直結してくるところではないのでそこが採点の基準にならないから順番は気にしなくてもいいんじゃないかと

タグ:

posted at 16:27:34

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 正規分布モデルの場合だけ、交差検証の結果がとても悪いものになってしまった理由と、他の裾野が重い分布を使ったモデルではそうなっていない理由を知るために、サンプルサイズが7しかないことを利用して一個抜き予測分布をすべてプロットしてみました。続く

タグ: JuliaLang 数楽 統計

posted at 16:25:29

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang

loocv_モデル名

の数値が以上で説明した計算の結果です。正規分布モデルの LOOCV の数値だけが跳ね上がっている!(数値が大きいほど悪いと考えられる。)

aic_モデル名

は所謂AICです。 pic.twitter.com/6lVLUbOAmz

タグ: JuliaLang 数楽 統計

posted at 16:23:29

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#統計 予測分布 p_i(x) によるi番目の数値の予測は p_i(i番目の数値) が大きいほど当たっていることになります。その対数の-2倍の和を計算しているので、計算結果が小さいほど予測精度が高いことになります。

タグ: 統計

posted at 16:18:24

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan 【5年生用と言うことで別の話かと】
何年生で解除になるのでしょうか?
明文化されていないのではないですか?
そのことが明文化されていれば、#掛算 の順序強制ハンターイ! と言っている人々の騒ぎはもうすこし沈静化するのですが。

タグ: 掛算

posted at 16:17:53

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 実際の計算では、大体において、サンプル

-27, 4, 8, 9, 10, 10, 10

からi番目を取り除いて作ったデータから予測分布 p_i(x) を作り、-2 log p_i(i番目の数値) の和を取ったものを計算してみました。

タグ: 数楽 統計

posted at 16:16:52

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 一個抜き交差検証(LOOCV)とは、もとのサンプル全体から数値を1つを取り除いて作ったデータから得られる予測分布で取り除いた数値をどれだけ予測できているかを確認することです。取り除く数値はサンプルに含まれる数値全体を動かす。

タグ: JuliaLang 数楽 統計

posted at 16:13:19

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang

nbviewer.jupyter.org/gist/genkuroki...
測定スキルの異なる7人の科学者たち

を更新しました。(一個抜き)交差検証(LOOCV)の節を付け加えました。

交差検証をしてみると、分析結果をどれだけ外れ値が支配しているかがわかります。

タグ: JuliaLang 数楽 統計

posted at 16:10:25

☆大天使たけにゃん☆ @Angel_takenyan

17年9月30日

@golgo_sardine 2年生までは順序固定にしてるとあるのでこの回答を作った人も昨日の発言のようなクレーム対策として書いてるんだと思います。
実際昨日の発言は小学2年生のテストの回答で順番が違うから×になったという記事に対するもので今回は5年生用と言うことで別の話かと

タグ:

posted at 13:17:54

森下屋 菊助 @morishitaya

17年9月30日

@4rikodama @genkuroki @RyoUchida_RIRIS ちさな事だからと握りつぶしていたら、いずれ学校に警備員や警察官を常駐させなくてはならなくなるかもしれません。

学校にも割れ窓理論を活用してこのような事を減らす努力をしないとなりませんね。🤔

タグ:

posted at 13:07:37

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan どうして、餅の数を求める計算が
 8臼×55個×6学年
なんでしょうね? #掛算
e-school.e-tokushima.or.jp/awa/es/donari/... pic.twitter.com/3BcxGbOteU

タグ: 掛算

posted at 12:42:07

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan 「10分間の2倍 は 10×2」
「10分間の3倍 は 10×3」が正しいのでしょうか?
どうして「自宅での学習時間は 学年×10分」と言われるのでしょうか?
www.pref.chiba.lg.jp/kyouiku/shidou...
#掛算 pic.twitter.com/T7UDtlRhiW

タグ: 掛算

posted at 12:38:56

ゴルゴ・サーディーン @golgo_sardine

17年9月30日

@Angel_takenyan いきなりで失礼します。
【タイヤの数が前に来なきゃいけないから2×5なんよね】
そう言われるかたにはこれをお聞きしたいと思います。
この受験研究社の本は間違っているのでしょうか? #掛算 pic.twitter.com/CZHel8lZq3

タグ: 掛算

posted at 12:34:12

黒木玄 Gen Kuroki @genkuroki

17年9月30日

まあ、いずれにせよ、「国も借金で大変なので予算削減も仕方がない」「自分達だけ予算を要求できない」とか言っている大学の先生達はどちらかと言えば(大学の行く末と無関係に)社会的に負の貢献をしてしまっていることは間違いないと思う。財政緊縮を当然視する発言を平気でできることが大問題。

タグ:

posted at 12:26:38

黒木玄 Gen Kuroki @genkuroki

17年9月30日

twitter.com/yamauchitaiji/...
マクロ経済政策との関係を理解せずにこういう反応をする人達は、まさに「知的レベルに問題がある人達」(カギ括弧付き)だと思うのですが、どうなんでしょうかね?

自覚が無いのか、それとも自分にとって都合が悪いからとぼけているのか?

タグ:

posted at 12:21:55

七誌 @7shi

17年9月30日

回線が繋がらないときにテザリングしまくったので、今月は結構容量を使った。いつもは1GBすら行かないので契約をミスったかと思っていたけど、お陰で助かった… pic.twitter.com/DkuSYQtnG0

タグ:

posted at 12:02:43

@kuri_kurita

17年9月30日

たとえば、(アベノミクスが始まる前の)不況期に実質賃金が「上がった」のは、低賃金の非正規雇用者からまず解雇されていくため、平均の賃金が上がったからだとか。

twitter.com/kuri_kurita/st...

タグ:

posted at 11:24:30

@kuri_kurita

17年9月30日

雇用状況が改善し始めた時には実質賃金というのはまずいったん下がるものなのですが、「実質賃金ガー」の人も「ドル建てGDPガー」の人と同じく、まず何を言われても理解を拒みますよね。なぜなんだか。

しかしどうやら実質賃金も上がり始めてきたようですけど。 twitter.com/ranrando/statu...

タグ:

posted at 11:18:32

Haruhiko Okumura @h_okumura

17年9月30日

High SierraにしてからRのas.POSIXct()の動作が変だ。Sys.timezone()はNAを返すし

タグ:

posted at 11:10:58

黒木玄 Gen Kuroki @genkuroki

17年9月30日

twitter.com/corollary2525/...
twitter.com/genkuroki/stat...
ラピュタもぼくがお気に入りのプロットもでっぱりは7個だじぇ!

タグ:

posted at 10:55:43

山内太地『偏差値45からの大学の選び方』 @yamauchitaiji

17年9月30日

シビアな話ですが、自分でお金を稼ぐ能力のある人が生き残るのです。 twitter.com/genkuroki/stat...

タグ:

posted at 10:51:23

黒木玄 Gen Kuroki @genkuroki

17年9月30日

で、この手の発言の方が、私の好きな数楽話よりも注目をあびる怖さ。

タグ:

posted at 10:41:11

黒木玄 Gen Kuroki @genkuroki

17年9月30日

自分の方が知的に優っていると思っている人達が実は知的なレベルの低さのせいで政治的に負けている場合の典型例は、「国も借金で大変なので予算削減も仕方がない」「自分達だけ予算を要求できない」とか言っている大学の先生達だな。そろそろ自分自身の知的レベルの低さに気付くべきだ。

タグ:

posted at 10:39:59

黒木玄 Gen Kuroki @genkuroki

17年9月30日

「安倍政権をぶっ潰せ」と運動するのはよい。大賛成だ(笑)。しかし、雇用を増やすことに成功した政策=金融政策を重視せずに「ぶっ潰せ」と言う人達は極めて有害。緊縮財政でもこれだけの成果を出せたのだから、財政赤字を思いっ切り増やす政策をやればさらに素晴らしい成果を出せるだろう。

タグ:

posted at 10:36:36

グレッグ @glegory

17年9月30日

金子洋一先生は素晴らしいよ。沈没寸前の民進党に落選の身でいながら、今回の騒動の中でも微動だにしなかった。腹の据わり方が他の議員とは違う。そう、彼はモノが違う。うろうろしている奴ばっかりだから、尚更彼の凄さが際立つ。自民党は是非、彼をリクルートすべき。

タグ:

posted at 10:01:59

Nyoho @NeXTSTEP2OSX

17年9月30日

@taketo1024 なるほどGL(n,k)の境界でしたか〜! O(n)も楽しみにしています!

タグ:

posted at 09:03:16

さのたけと @taketo1024

17年9月30日

@NeXTSTEP2OSX これは GL の二つの連結成分の境界を見たかったので =0 にしました。SL の形を見たければ =1 ですね👍

タグ:

posted at 09:00:19

Nyoho @NeXTSTEP2OSX

17年9月30日

@taketo1024 これってイコール0なんですか? イコール1じゃなくて

タグ:

posted at 07:38:07

Masahiro Hotta @hottaqu

17年9月30日

昨日の #QTTF2017 のバンケットでは、東大の上田さんに冷却原子系実験の最近の凄まじい進歩の話をお聞きしているうちに、思わぬ展開があって嬉しい驚きあった。4次元ブラックホール蒸発の熱的性質を再現する、僕らが作った崩壊N体量子ビットモデルの実験ができる可能性を指摘して頂いた。

タグ: QTTF2017

posted at 05:14:31

あ〜る菊池誠(反緊縮)公式 @kikumaco

17年9月30日

ベイズと統計力学が数学的に等価であること自体は渡辺さんの発見ではないと思ったのだけどな twitter.com/hayashiyus/sta...

タグ:

posted at 02:33:56

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 階層モデルを内部変数 u について積分してしまえば、普通に最尤法も使えます。内部変数 u で積分した結果がよく知られている分布になる階層モデルの分類表があると便利だと思いました。すでにありそうな気がしないでもない。あるならばありかを教えて下さい。

タグ: 数楽 統計

posted at 02:19:40

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 たぶん以上のような作業は #JuliaLang を使うよりも、Rにおける既存のパッケージを使った方がずっと楽に行えると思います。しかし、統計用のパッケージにあまり頼らずに様々な量を計算する函数を自前で書くと勉強になると思う。

タグ: JuliaLang 数楽 統計

posted at 02:04:34

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 裾野が重い t 分布を使ったモデルは、「各々の科学者は測定スキル σ² を持っており、測定値の分布は真の値を中心とする正規分布に従っている。そして、科学者ごとの σ² は逆ガンマ分布に従ってばらついている」というモデルと同値。

タグ: JuliaLang 数楽 統計

posted at 02:01:47

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliLang 外れ値を出した実験が苦手な科学者がめっちゃすき(笑)

タグ: JuliLang 数楽 統計

posted at 01:52:00

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliLang 正規分布モデルの予測分布のプロットにサンプルの位置も付け加えました。正規分布モデルだとこういう感じになります。「外れ値を出した科学者は特別に実験が苦手である」と考えるなら、他の裾野が重いモデルを採用する必要があります。 pic.twitter.com/fwrvfuitBh

タグ: JuliLang 数楽 統計

posted at 01:43:25

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliLang 正規分布モデルの予測分布のプロット。他の3つと同時にプロットすると平べった過ぎるので単独でプロットしてみました。外れ値のせいで、μ=3.5、σ=13になってしまっています。裾野が重い他の3つのモデルでのμはどれも10弱になっています。 pic.twitter.com/ubfdwQFUQK

タグ: JuliLang 数楽 統計

posted at 01:38:21

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang ベイズ推定のMCMCの方法ならば内部変数で積分する前の情報をコンピューターに教えてあげれば事後分布が得られます。しかし、AICのベイズ版のWAICを計算しようとすると内部変数での積分が必要になってしまうのだ。その辺は面白い問題だと思います。

タグ: JuliaLang 数楽 統計

posted at 01:27:55

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 階層モデルとの関係の解説も含めて、詳しいことについては、

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

を見て下さい。「最尤法を使う場合であっても尤度函数と予測分布の形をきちんと他人に見せる」という習慣になるとよいと思う。

タグ: JuliaLang 数楽 統計

posted at 01:24:43

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 「科学者ごとに測定したときの分散σ²が違っていて、σ²は逆ガンマ分布に従ってばらついている」という仮定に基づいたモデルは、内部変数のσ²について積分すれば、我々が扱ったt分布を使ったモデルと同値になります。

タグ: JuliaLang 数楽 統計

posted at 01:23:15

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 以上で扱った裾野が重いt分布を用いた確率モデルは次の形の階層モデルと同値です。

y ~ Normal(μ,σ²)
σ² ~ invGamma(θ,k)
θとkはパラメーター (ベイズ推定では事前分布を与える)

続く

タグ: JuliaLang 数楽 統計

posted at 01:21:21

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
4つの確率モデルを最尤法で解いた結果を紹介したのですが、それらのモデルは単純な正則モデルなので、ベイズ推定で解いても(すなわちベイズ推定の方法で予測分布を求めても)ほとんど同じ結果しか得られないと思います。

タグ: JuliaLang 数楽 統計

posted at 01:17:56

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
予測分布のプロット
外れ値はこの範囲のプロットに含まれていません。

正規分布モデルによる予測分布は平べったく薄く広がってしまっていています。 pic.twitter.com/PfYRhokuyy

タグ: JuliaLang 数楽 統計

posted at 01:15:36

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
「自由度νのt分布に従う確率変数のσ倍足すμ」の尤度函数
これは以上の2つのモデルの同時拡張になっています
パラメーターが3つあるので切断面をプロットしました。
さらにμ方向の幅が狭くなりました。 pic.twitter.com/RICISjaBIE

タグ: JuliaLang 数楽 統計

posted at 01:13:36

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
「自由度νのt分布に従う確率変数足すμ」のモデルの尤度函数
Cauchy分布によるモデルほどではないのですが、μ方向の幅が狭くなっています。 pic.twitter.com/yN3V40ONoG

タグ: JuliaLang 数楽 統計

posted at 01:11:27

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
「Cauchy分布に従う確率変数のσ倍足すμ」のモデルの尤度函数
外れ値の影響が減って、μの範囲が狭くなっています。 pic.twitter.com/qdUbMbEKva

タグ: JuliaLang 数楽 統計

posted at 01:09:34

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
正規分布モデルの尤度函数
平均値μの範囲がものすごくひろい。
外れ値の影響です。 pic.twitter.com/vLQPvpYoxu

タグ: JuliaLang 数楽 統計

posted at 01:07:52

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang サンプルサイズがn=7のケースなので、十分信頼できることは何も言えないと思うのですが、ひとまずの判断をしたい場合にはこんな感じで分析をするしかないと思う。

タグ: JuliaLang 数楽 統計

posted at 01:03:07

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 正規分布モデルのAIC(もしくはもっと単純に対数尤度の-1倍)は裾野が重い確率分布を用いた残りの3つと比較すると、かなり大きな値になっています。AICは小さい方がよいので、正規分布モデルは予測精度が他の3つより低いと推定されます。

タグ: JuliaLang 数楽 統計

posted at 01:02:02

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 我々が扱っているケースにおける7人の科学者の測定結果から作った「予測分布」は、さらにもう1人の科学者が同一の物理量を測定したときの測定結果の推定確率分布になります。その予測精度はAICの小ささで推定できる。

タグ: JuliaLang 数楽 統計

posted at 01:00:06

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 尤度函数はベイズ統計における一様事前分布のもとでの事後分布と(定数倍の違いを除いて)一致するので、「尤度函数も見せろや!」と言うことは、「ベイズ統計の結果も見せろや!」と言うことにかなり近いです。「予測分布」についても同じ雰囲気があります。

タグ: JuliaLang 数楽 統計

posted at 00:58:07

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 私は「尤度函数も見せろや!」「予測分布はどうなったんだ?」と言いたくなるタイプなので、最尤法の結果だけではなく、尤度函数そのものもプロットしているし、予測分布もプロットしています。みんな、そうして欲しいと思います。

タグ: JuliaLang 数楽 統計

posted at 00:55:26

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 「岡田先生資料」はベイズ統計推しで、私もベイズ統計は好きなのですが、今回は4種類の確率モデルを使った最尤法で推定してみた。

正規分布モデル
Cauchy分布モデル
t 分布を使った2種類のモデル

外れ値には裾野が重い分布が便利です。

タグ: JuliaLang 数楽 統計

posted at 00:53:22

☆大天使たけにゃん☆ @Angel_takenyan

17年9月30日

掛け算の件についてはこの書き方だと5×2でもよく見えるけど単位がタイヤの個数だからタイヤの数が前に来なきゃいけないから2×5なんよね

タグ:

posted at 00:51:52

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang 7人の科学者が同一の物理量を測定したら、それぞれ次の数値を得た:

-27.020, 3.570, 8.191, 9.808, 9.603, 9.945, 10.056

-27.020というひどい外れ値が含まれています。続く

タグ: JuliaLang 数楽 統計

posted at 00:49:51

黒木玄 Gen Kuroki @genkuroki

17年9月30日

#数楽 #統計 #JuliaLang
drive.google.com/drive/folders/...
の「岡田先生資料」の例2に4つの確率モデルを適用し、最尤法で解いてみました。
「7人の科学者の測定の中にひどい外れ値がぁー!」の問題です。
nbviewer.jupyter.org/gist/genkuroki...

タグ: JuliaLang 数楽 統計

posted at 00:48:10

さのたけと @taketo1024

17年9月30日

まぁそれを見るためには Grapher という便利なソフトが最初から入ってて…😇 pic.twitter.com/5eILFk7GlN

タグ:

posted at 00:44:47

さのたけと @taketo1024

17年9月30日

先に GL(2; R) の [-1, 1]^4 の部分を見てみた。w = 0 のとこで境界は yz = 0 ってなってる感じは見えるけど、 w を動かしたときに境界がどう動くか見えないと連結成分が 2 個なのは見えづらい。 pic.twitter.com/yj40oi3xxp

タグ:

posted at 00:29:05

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

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

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

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

送信中

送信に失敗しました

タグを編集しました