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 » 2020年12月24日
並び順 : 新→古 | 古→新

2020年12月24日(木)

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 Distributions.jl に fit_mle(Beta, X) を追加する方法。

Optim.jl を何の工夫も無しに使う場合。

工夫すればもっと高速化できるはずです。

gist.github.com/genkuroki/5ecd... pic.twitter.com/OPWZK6NoM9

タグ: Julia言語

posted at 00:08:48

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 の使用時の気楽さは本当に独特な感じ。

タグ: Julia言語

posted at 00:15:17

Katsushi Kagaya @katzkagaya

20年12月24日

シャコパンチの研究でPatek研にフィルさんという古生物学者がいたんだけど「物理法則は変わってないけど生物は基本的に変わってきた」ということを言っていて、なるほどなぁと。サカナのアゴやシャコの捕脚を梃子とみると、形がわかればその比から機械効率が分かってそれがどう進化したか分かる、と。

タグ:

posted at 00:26:24

ファインマンbot @feynmannnn

20年12月24日

何かを見付けようというときには、深遠な哲学的議論をたたかわすよりも、注意深く実験をする方がよい。

タグ:

posted at 00:34:04

積分定数 @sekibunnteisuu

20年12月24日

これ、死刑判決書いた元裁判官が、「間違っていた」と言っているんだよね。 twitter.com/nhk_news/statu...

タグ:

posted at 00:59:19

積分定数 @sekibunnteisuu

20年12月24日

ja.wikipedia.org/wiki/%E7%B4%85...
>担当した事件で数多くの冤罪を作ったことで知られる。

タグ:

posted at 01:01:28

積分定数 @sekibunnteisuu

20年12月24日

しかしこれ、現役教員の発言だけどドン引きするよね。 twitter.com/MegJoytoy/stat... pic.twitter.com/VW6ziBiwI1

タグ:

posted at 01:22:16

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy >おかしな行為だと思う理由ですか?文科省が定めたテストに文句を言っているからです。

と言っていますよ。

文科省に文句を言うのはおかしな行為、という認識なんですよね?

twitter.com/MegJoytoy/stat...

タグ:

posted at 01:28:33

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy 学校教育法 

>第五十一条 高等学校における教育は、前条に規定する目的を実現するため、次に掲げる目標を達成するよう行われるものとする。

>三 個性の確立に努めるとともに、社会について、広く深い理解と健全な批判力を養い、社会の発展に寄与する態度を養うこと。

タグ:

posted at 01:29:23

黒木玄 Gen Kuroki @genkuroki

20年12月24日

@Salieri91381586 #超算数 文科省が教師に掛算順序指導を強制しているという事実はありません。

基本的に小学校の先生は自由意志で子供を傷つけるようなひどい教え方をしているとみなされます。そのような教師の各個人が社会的に厳しく批判されるのは仕方がないことです。

続く

タグ: 超算数

posted at 01:29:38

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy >健全な批判力を養い

タグ:

posted at 01:29:41

黒木玄 Gen Kuroki @genkuroki

20年12月24日

@Salieri91381586 #超算数 ただし、小学校内その他の周囲の圧力に負けて、掛算順序でバツをつけてしまう先生もいます。

そのようなケースでは外部から見えない小学校内部での教師に対する不当な指示を表沙汰にして、指示を出した側に厳しい処分が下されることが妥当だと思いますが。

タグ: 超算数

posted at 01:31:50

黒木玄 Gen Kuroki @genkuroki

20年12月24日

@Salieri91381586 #超算数 例えば、学年主任の先生が、若い先生が「掛算順序でバツをつけるのは間違っている」と言っているのに、無理矢理強権的にバツを付けさせたというような場合には、その学年主任の先生に社会的な非難が集中するべきだと思います。

タグ: 超算数

posted at 01:34:02

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy つまり、文科省が「カラスは白い」と言ったら、カラスは白いと信じると言うことですね。

で、文科省が「やっぱりカラスは黒でした」というと、あなたも「やっぱりカラスは黒でした」と言うのですね。

タグ:

posted at 01:34:04

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy >わたしは組織を大切にするタイプ

国家も会社もやくざ組織も、トップの周りがイエスマンばかりで組織がつぶれる、というのはよくあると思うのですがね。

タグ:

posted at 01:35:12

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy >おかしいというか謎です。

つまりあなたは、文科省を批判する人がいると言うことが謎ということでしょうか?

そういう認識をどこで身に付けたかの方が謎なんですが。

タグ:

posted at 01:36:03

黒木玄 Gen Kuroki @genkuroki

20年12月24日

@Salieri91381586 #超算数 私がリアルで知っている範囲にも、保護者側に「掛算順序でバツをつけるにはおかしい」とはっきり言ってくれる良い先生がいます。

そういう先生達が周囲の圧力に負けずに普通の教育ができる環境が整うようになれば良いと思います。

タグ: 超算数

posted at 01:36:13

黒木玄 Gen Kuroki @genkuroki

20年12月24日

そういう経験は皆無。 twitter.com/phasetrbot/sta...

タグ:

posted at 01:41:53

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語

Juliaは begin-origin.

例えば0始まりだったり、-100始まりのインデックスであることもあり得る。

a[begin] で配列aの最初の要素にアクセスできる。その次の要素は a[begin+1].

size(a)ではなく、axes(a)で配列のインデックスの範囲のタプルが得られる。 twitter.com/cocomoff/statu...

タグ: Julia言語

posted at 01:45:26

積分定数 @sekibunnteisuu

20年12月24日

@MegJoytoy あなたみたいな人ばかりが警察組織にいたら、上司の暴走を止める人がいなくて、冤罪が多発しそうですね。

タグ:

posted at 01:47:26

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語

別言語で0始まり配列にa[i]とアクセスしている部分は、Juliaではa[begin+i]に書き換えればよい、という感じだと思う。

OffsetArrays.jl を使っても良いが、線形代数で困る。
その辺もまだ発展の余地が残っている感じ。

タグ: Julia言語

posted at 01:51:00

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 2次元配列aについて

m, n = size(a)
for j in 1:n
for i in 1:m
a[i,j] = f(i, j)
end
end

のように書く代わりに

is, js = axes(a)
for j in js
for i in is
a[i,j] = f(i,j)
end
end

と書けば、1始まりでない2次元配列にも対応できる。

タグ: Julia言語

posted at 01:53:54

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 おまけ

Pythonからの引っ越し組で、速さが欲しい人は

a = []
for i in is
x = f(i)
push!(a, x)
end

は避けた方がよい。

aがAnyの配列になってaを使った数値計算が遅くなってしまう。

a = typeof(f(first(is)))[]
以下同じ

とか

a = f.(is)

とか。

タグ: Julia言語

posted at 01:59:33

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 私はこうしている。

⓪ a = f.(x)

①最終的にできる配列のサイズがわかっているなら、配列を一挙に確保して各成分に順次代入する。

②速度無視の場合または作られる配列にサイズが計算終了まで分からない場合には、a = (型指定)[] で空の型付配列を作ってから、push!を繰り返す。

タグ: Julia言語

posted at 02:03:57

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 配列 a の型が Vector{Real} になっているが、すべての成分が Float64 になっているならば、 a = Vector{Float64}(a) で型を変換してからaを使用した方が計算が速くなることがある。

例えばaがMCMCで得たサンプルの類で、そのサンプルを使って期待値を計算する場合とか。

タグ: Julia言語

posted at 02:09:16

ファインマンbot @R_P_Feynman_bot

20年12月24日

権威は尊敬に値しない。

タグ:

posted at 03:49:36

Noboru Hagino @Noboru_Hagino

20年12月24日

例によって例のごとく、本筋と関係ないところで用語の独自解釈があって、そちらにイラッとさせられる。
yomidr.yomiuri.co.jp/article/202012...

タグ:

posted at 04:49:36

Dean Markwick @DeanMarkwick

20年12月24日

Latest blog post, converting a machine learning model from #Python to #JuliaLang dm13450.github.io/2020/12/19/Edw...

タグ: JuliaLang Python

posted at 05:17:06

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#Julia言語 とても良い!

Julia 1.6の個人的に気になる新機能 qiita.com/ho-oto/items/e... #Qiita

タグ: Julia言語 Qiita

posted at 05:38:58

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#数楽 反応拡散方程式の大幅な単純化であるYoungモデルの解説

#Julia言語 で遊ぶネタが欲しい人向け情報。

17-7. Youngのモデルでお手軽にチューリング・パターン #ビネクラ vigne-cla.com/17-7/

タグ: Julia言語 ビネクラ 数楽

posted at 05:53:50

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計

不完全ベータ函数も統計学で必要な基本特殊函数なのですが、現在よく使われているコードはこれ↓

github.com/JuliaStats/Rma...

2000行以上ある!

こういうのを「もっと整理しろ」と言われると泣きそう。

不完全ベータ函数がコンピュータで計算できることは結構重要。

タグ: 統計

posted at 06:00:42

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 表の出る確率がθのとき

(n回中表の回数がk以下の確率) = 1 - I_θ(k+1, n-k)

(k回表が出るまでの試行回数がn回以上の確率) = 1 - I_θ(k, n-k)

となることは、仮説「表の出る確率はθ以上」の片側検定のP値がベイズ統計でのimproper事前分布の事後分布で仮説が成立する確率に等しいことを意味。

タグ: 統計

posted at 06:19:30

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 ここで、I_z(a,b)は正則化された不完全ベータ函数:

I_z(a,b) = ∫_0^z x^{a-1}(1-x)^{b-1} dx/B(a,b).

仮説検定でのP値がベイズ統計での事後分布内での「仮説が成立する確率」にぴったり一致するので、

P値の代わりにベイズ統計での「仮説が成立する確率」を使え

という主張は

無意味!

タグ: 統計

posted at 06:19:31

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「頻度論のP値は分かりにくいが、ベイズ統計における事後分布で計算した仮説が成立する確率は分かり易い」とか、「頻度論の95%信頼区間の95%は確率ではないので分かりにくいが、ベイズ統計における95%信用区間の95%は確率なので分かり易い」の類の主張をしている人達は、理解していない人達。

タグ: 統計

posted at 06:30:06

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 一般に、困った人達が「頻度論」「頻度主義」と呼んでいる事柄(P値や信頼区間)の統計学入門の教科書における説明はひどくずさんで分かりにくいことが多い。正確な理解も一般的な性質の理解も困難な説明の仕方になっている。

だから、P値や信頼区間は分かり難いと誤解されているのだと思う。

タグ: 統計

posted at 06:30:07

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計

①統計学入門の教科書ではP値や信頼区間の説明は数学的前提がクリアに説明されておらず、P値と信頼区間の関係さえ説明されいない。

②読者は理解不可能になる。

③FisherとNeyman-Pearsonで検定の意味が違うというような無用なことを言い出す人が偉そうにしていてさらに混乱が広がる。

タグ: 統計

posted at 06:36:04

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 P値や信頼区間について正確な説明をしてもらったことがない人達がそれらを誤用するのは当たり前。

FisherやNeyman-Pearsonはそういう初歩レベルの話題には無関係。

そして、P値や信頼区間について正確に理解しているのに、それらを誤用している場合は論理的能力と科学的常識の欠如を疑うべき。

タグ: 統計

posted at 06:42:03

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#数楽 とてもいい話だ。

教員志望の学生にこういう話を伝えて行きたい。 twitter.com/ototo_/status/...

タグ: 数楽

posted at 06:46:53

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 このスレッドでは信頼区間について扱わなかかったが、パラメータθとデータDをP値に対応させる函数pval(θ,D)が与えられているとき、信頼度1-αの信頼区間は「pval(θ,D)≧αとなるθ全体の集合」とシンプルに定義できます。

要するに有意水準αで棄却されないθの範囲が信頼区間になる。

タグ: 統計

posted at 06:57:54

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 これを知っていると、pval(θ, D)がコンピュータ上ですでに実装されている場合には、本質的に、函数 θ ↦ pval(θ, D) - α が0になるθ達を求めてくれるライブラリを使うことで、信頼区間も計算できるようになります。

自前で信頼区間を計算する函数のコードを書きたい人には有用な知識。

タグ: 統計

posted at 07:07:37

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「ベルヌイ分布モデルの場合に、何種類も pval(θ, D) を書いて、それらのグラフをプロットし、信頼区間も求めてみる」のは手頃な練習問題になり得ます。

自分でコードを書いた方が理解が進むタイプの人はこういう方針で統計学に入門すると楽をできると思います。

タグ: 統計

posted at 07:10:42

Dave @SpotTheLoon2010

20年12月24日

We all need a little help sometimes pic.twitter.com/FBhrydMlGt

タグ:

posted at 08:25:09

Yuki Nagai @cometscome_phys

20年12月24日

LGのテレビ、標準でAirPlayがついているので無線でMacのディスプレイにできる。55インチにスライドを表示させたら研究室のプロジェクターにうつるスライドに近いサイズ感

タグ:

posted at 08:49:05

中村 良幸 (Nakamura Yosh @nakayoshix

20年12月24日

@shigerufujita 導入教育をRustにするかそれともJuliaにするかは、情報工学向けの情報教育なら前者、それ以外の(一般の)理工学向けの情報教育なら後者という感じですかね…

(文系向けの導入教育にはやっぱりPython?w)

タグ:

posted at 09:54:23

いほやきゅ @K0TSU_K0TSU

20年12月24日

内部構造描いてみました

カニのハサミの中には腱が2本あります、動画では開く時も閉じる時も赤の腱を動かしてますが、実際は閉じる時だけ赤い腱が動きます、逆に開く時は緑の腱が動きます、赤の腱は面積が大きく大量の筋肉が付着しています、つまり閉じる力が大きいことが分かります

※大体の形です pic.twitter.com/8KQnoJCP0z

タグ:

posted at 10:23:42

古庄严 @iwaofurusho

20年12月24日

周→春秋戦国勢力変遷を30秒で勉強

春秋戦国時代は、中国史において、紀元前770年に周が都を洛邑(成周)へ移してから、紀元前221年に秦が中国を統一するまでの時代。 pic.twitter.com/V47jo8YCNz

タグ:

posted at 10:42:11

非公開

タグ:

posted at xx:xx:xx

統計たん @stattan

20年12月24日

死の商人よろしく、主義の対立を煽ることで論文を増やすことができる人たちがいる。

タグ:

posted at 11:16:58

あみあみ @amiami114114

20年12月24日

小学校は分からない子はそのままで、分かっている子をさくらんぼ計算やくもわ・はじきなどで「分からない子」に仕立てているように見えます。分からなさに向き合っているようには見えない。 twitter.com/kiya__na/statu...

タグ:

posted at 11:19:33

家庭教師curren @curren49590380

20年12月24日

まずは運動部で教えたいから社会や国語の中学校教諭になりたい

という生徒を減らしていかないといけないですよ。(彼らは決して体育とは言わない)

タグ:

posted at 12:03:29

砂___の___女 @vecchio_ciao

20年12月24日

次男が小5当時の「振り子」案件。

この時の教員の「実験した結果とは違うけど、教科書に書いてある通りに書き直せ」と言った指導は、多くの小学校教員に象徴的な学力不足が原因だったと思う。 twitter.com/vecchio_ciao/s...

タグ:

posted at 12:11:04

非公開

タグ:

posted at xx:xx:xx

砂___の___女 @vecchio_ciao

20年12月24日

Twitterで皆さんに教えていただいて、実験の条件・道具に重大な不備があったのが大きな原因だとわかって、次男も納得した。
担任は朗らかでいかにも性格良さそうな人だったが、残念ながら教科指導力には決定的に欠けていた。

タグ:

posted at 12:14:00

砂___の___女 @vecchio_ciao

20年12月24日

いかにも「性格良さそう」に見えたのは、次男が担任のお気に入りの生徒であったので、親の私に対してもいつも笑顔で対応してくれたから、というバイアスは否めない。 twitter.com/vecchio_ciao/s...

タグ:

posted at 12:17:05

砂___の___女 @vecchio_ciao

20年12月24日

「センター7割」に“噴き上がってた”教員アカウントは、「全員に要求するのは無理」などと謙虚なことを言ってなくて、「学力よりも大切なものがある」「学力高い人よりもむしろ初等教育者として能力が高い」と言ってたのを知らずに入ってきた人かな? twitter.com/sgym50600406/s...

タグ:

posted at 13:48:33

非公開

タグ:

posted at xx:xx:xx

Tadashi Nagasawa@Arm @Sonopapa

20年12月24日

バイデン次期大統領が、年収12万5千ドル(≒1300万円)以下の家庭には、公立大学の学費を免除する、という政策を発表したようだ。 twitter.com/JoeBiden/statu...

タグ:

posted at 14:33:56

砂___の___女 @vecchio_ciao

20年12月24日

「センター7割」の件を見て、教員志望で7割ボーダーくらいの(少しできる、そして謙虚な)学生さんが「自分に教員は向いてない…」と断念して居なくなり、それ以下しか残らない、さらに低レベル化が加速する未来を懸念している。

タグ:

posted at 14:35:48

Tomoya TENDO @end_tt

20年12月24日

盛田みずすまし @nosiika さんのアイデアで、ねじれ正多面体
ja.wikipedia.org/wiki/%E3%81%AD...
のような物を菱形十二面体でもできないか試してみました。
菱形の鋭角の頂点では6面、鈍角の頂点では4面が集まるようにうまく繋げると、空間を同体積に2分割できるようです。
全体が連動するので動かしても楽しい。 pic.twitter.com/l4tH4XMbOn

タグ:

posted at 14:47:11

カイヤン @389jan

20年12月24日

@ceptree @TomiyaAkio 世界のSumio Watanabe「両方の視点で量子統計力学reviewを書きました」
www.jstage.jst.go.jp/article/soken/...

タグ:

posted at 15:00:36

増田の准教授 @ProfMasuda

20年12月24日

底抜けのナントカ。単なる破壊に終わると思う。
「痛みを伴う改革を断行する大学長や自治体首長のリーダーシップを重視。」
学部解体の断行も!地方国立大学は破壊的改革へ(ニュースイッチ) - Yahoo!ニュース news.yahoo.co.jp/articles/aa3c5...

タグ:

posted at 15:01:46

Steven Pinker @sapinker

20年12月24日

The Best Illusions of the Year are 2020's Last Attempt to Break Your Brain gizmodo.com/the-best-illus... via @gizmodo

タグ:

posted at 15:09:51

Shuhei Kadowaki @kdwkshh

20年12月24日

影響受けまくってるので、JET.jlのREADMEでruby/typeprofをreferenceした:

github.com/aviatesk/JET.j...

タグ:

posted at 15:14:55

matheca @paulerdosh

20年12月24日

私小学生のとき父兄参観日の授業中教諭からこういうの食らって、とても傷つきました。図書館で確認しても私の方が正しかったのに。 twitter.com/esumii/status/...

タグ:

posted at 15:27:28

abap34 @abap34

20年12月24日

Juliaでカラフルに点滅するクリスマスツリーを書いてみました
#Julia言語 pic.twitter.com/eZZr8fCgAe

タグ: Julia言語

posted at 15:30:02

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

20年12月24日

初等教育での「振り子の等時性」問題は深刻な気がするな。等時性は小振幅での近似的性質だということを教員が理解していないと、誤った教育が行われる

タグ:

posted at 15:30:17

abap34 @abap34

20年12月24日

出力を戻すところはごまふさんのコードを拝借した

タグ:

posted at 15:31:28

積分定数 @sekibunnteisuu

20年12月24日

そうそう、ハジキって、能力のない教員が、教育委員会の目を盗んでこっそり教えているとかじゃないんだよね。三島市教委と話をした時も、「ハジキを教えることがあるけど、意味も教えている」という言い方で、「ハジキで教えないようにと教員に言っているのだが徹底できていない」ということではない。 twitter.com/jarajaraKei/st...

タグ:

posted at 15:38:13

@kuri_kurita

20年12月24日

しかも「正しい」観察結果を得た子供が、「自分の予想に合うようにデータを処理した」と教師に責められる可能性すら。 (予想に合うようにデータをいじっているのはどっちだ、という超絶皮肉な倒錯) twitter.com/kikumaco/statu...

タグ:

posted at 15:44:35

Akinori Ito @akinori_ito

20年12月24日

地方国立大学を閉学した最後の学長はたくさん退職金がもらえるとか、そういうやつ?

タグ:

posted at 15:47:12

Akinori Ito @akinori_ito

20年12月24日

教授を一人リストラすると学長の給料が100万円ずつ上がっていくとか

タグ:

posted at 15:47:56

Akinori Ito @akinori_ito

20年12月24日

苦しい会社に金融畑の社長が送り込まれるとまずコストカットするよね(類推)

タグ:

posted at 15:49:56

あやう@10y5y @ayau15

20年12月24日

かけ算の順序問題、担任から連絡があった。
かける数とかけられる数の順序が入れ替わると意味が違ってしまうとのこと。
しかもひとつあたりの数×いくつ分の順序で統一して教えてるから、逆になると本人が苦労するとか、逆を正解とすると他の保護者からクレームになるとか、そんな言い訳ばかり。→

タグ:

posted at 16:00:31

l_ppp @ppp3141592ppp

20年12月24日

くそっ、タプル使えば

a,b=b,a

で値の交換できるじゃないか
俺のtmpは何だったんだ

タグ:

posted at 16:10:25

akky @tall_tree_desu

20年12月24日

@yuky_az RustやJuliaならニーズはあるかもしれませんね

タグ:

posted at 16:12:25

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「尤度」=「もっともらしさ」と説明している人は何も理解せずに説明しているので要注意。

尤度の定義は「モデル内で観測されたデータと同じ数値が生成される確率または確率密度」です。モデルのデータへの適合の良さの指標と言って良いのですが、「もっともらしさ」と言えるとは限らない。 twitter.com/rkmsfc/status/...

タグ: 統計

posted at 16:45:25

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 例:「正規分布のi.i.d.でデータが生成されている」というモデルの確率密度函数は

p(x_1,…,x_n|μ,σ²) = p(x_1|μ,σ²)…p(x_n|μ,σ²)
(ここでp(x_i|μ,σ²) = exp(-(x_i-μ)/(2σ²))/√(2πσ²)).

観測データの数値X_1,…,X_nが得られたとき、このモデルの尤度は

p(X_1,…,X_n|μ,σ²)

になる。続く

タグ: 統計

posted at 16:45:26

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 尤度が大きいほどモデルはデータによくフィットしていると思ってよいです。

観測データに関するモデルの尤度をモデルのパラメータの函数とみなしたもの、例えば上の例では

L(μ,σ²) = p(X_1,…,X_n|μ,σ²)

を観測データX_1,…,X_nに関するモデルの尤度函数と呼びます。

続く

タグ: 統計

posted at 16:45:26

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 尤度の大きさは、モデルのデータへの適合の良さを意味すると考えて良いので、尤度函数を最大化するモデルのパラメータを求めることは、モデルをデータに最もよく適合させるパラメータを求めているとみなせます。これを最尤法と呼びます。

続く

タグ: 統計

posted at 16:45:27

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 ただし、データに最も良くフィットするようにモデルのパラメータを調節しても、もっともらしい結果が得られるとは限らないことに注意が必要です。

モデルをデータにフィットさせ過ぎたせいで、モデルによる予測が外れまくるようになる場合があります。これをオーバーフィッティングと言う。

タグ: 統計

posted at 16:45:27

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 オーバーフィッティングが引き起こされていないという条件のもとで、最尤法はモデルの範囲内でベストに近い結果を近似的に得る方法として便利に使えます。

例えば、上で例に挙げたパラメータが2個しかない正規分布モデルの最尤法はそのような最尤法が有効な場合になっています。続く

タグ: 統計

posted at 16:45:28

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 仮に観測データが未知の真の分布のi.i.d.としてランダムに生成されているならば、nが十分大きなとき、上の正規分布モデルによって、未知の真の分布の平均と分散を正しく推定できます。

未知の真の分布に最もフィットする正規分布は同じ平均と分散を持つ正規分布になります。

タグ: 統計

posted at 16:45:29

黒木玄 Gen Kuroki @genkuroki

20年12月24日

しかし、未知の真の分布が正規分布から程遠い場合には、正規分布を正規分布と違う分布にフィッティングさせていることになり、予測の誤差は大きくなります。

モデルの適切さは非常に重要です。

タグ:

posted at 16:45:29

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 最尤法はモデルのパラメータ数が増えたり、構造を持ったモデル(特異モデル)では実用的でなくなります。

そういう場合には事前分布を使う方法がより実用的な推測法になります。

(主観確率でベイズ統計について解説している人達は社会的に負の貢献をしているので要注意!)

タグ: 統計

posted at 16:45:30

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 観測データに関するモデルの尤度函数を最大化するようなパラメータを求めることは、観測データにもっともモデルをフィットさせるようなパラメータを求めることになっている、と理解しておき、「もっともらしさ」のような曖昧でミスリーディングな言い方で理解したつもりにならないことが重要。

タグ: 統計

posted at 16:48:18

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「尤度」はlikelihoodの訳語で、その単語は「もっともらしさ」という意味を持つので、英語では「尤度=もっともらしさ」だとひどく誤解してしまう可能性が高くなっているものと思われます。

科学的合理性ではなく、歴史的な理由で不適切な単語が専門用語として固定されてしまうことがあります。

タグ: 統計

posted at 16:53:21

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 以上のように、「尤度=もっともらしさ」ではないとはっきり説明してくれればクリアなのに、統計学入門の教科書がそうなっていないせいで、多くの人が誤解したり、「尤度、理解できん」となっていると思う。

そこに「尤度原理」などの「主義」に関わるデタラメの解説が投入されると悲惨!

タグ: 統計

posted at 16:56:14

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 以下のリンク先動画は、正規分布とは違うが正規分布っぽくも見えるガンマ分布のi.i.d.で生成されたデータによる正規分布モデルの最尤法の様子です。

正規分布とは違う分布が正規分布で最良近似されて行く様子を見れます。 twitter.com/genkuroki/stat...

タグ: 統計

posted at 17:04:33

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計

データはドットで表示。データは y = sin(πx) + ノイズ で生成。

添付画像1: 3次式によるパラメータ数5個のモデルの最尤法によるフィッティング。最大対数尤度=4.2

添付画像2: 21次式によるパラメータ数23個のモデルの最尤法によるフィッティング。最大対数尤度=27.8 (上より大きい) twitter.com/genkuroki/stat... pic.twitter.com/lteSOILLlj

タグ: 統計

posted at 17:16:08

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 尤度は、添付画像1より2の方が高いです。

だから、「尤度=尤もらしさ」だと信じている人は、「添付画像2のフィッティングの方がもっともらしい」と言わなければいけません。

理屈を知らなくても、普通の感覚でそれはおかしいと思うのが普通だと思います。

これがオーバーフィッティング。 pic.twitter.com/GAgY1xPmDM

タグ: 統計

posted at 17:16:10

TaKu @takusansu

20年12月24日

@OokuboTact こういう資料を見ると、教科書改善の為に見る先は「生徒」ではなく「教員」になっていますね。

www.mext.go.jp/a_menu/shotou/...
【およそ3/4の教員(72.0パーセント)が、算数指導では「基礎・基本の定着」が最も重要と考えていることがわかる。】
この手の発想も、 #超算数 がのさばる要因かも。

タグ: 超算数

posted at 17:19:16

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 大事なことなので繰り返す。

尤度は観測データへのモデルの適合の良さ(フィッティングの良さ)の指標とみなせるが、「もっともらしさ」の指標ではない。

尤度を大きくすることは単に観測データにフィットするモデルを探しているに過ぎず、データにぴったりフィットすれば良いわけではない。

タグ: 統計

posted at 17:20:08

OokuboTact 大久保中二病中年 @OokuboTact

20年12月24日

@takusansu 一般的な小学校教員の想定する「算数の基礎・基本」はチョー算数ですから

タグ:

posted at 17:22:59

なべきち @nabekichi32

20年12月24日

「算数と数学は違う」という人の,数学をわかってない率は異常。

タグ:

posted at 17:34:40

Kouji Ohnishi 大西浩次 @koujiohnishi

20年12月24日

授業でも言っているように「#振り子の等時性」は #振幅#小さい ときに #近似的#単振動 に一致する事による。特に、短い糸(数十㎝)での振り子の実験をすると、振幅が変わることで大きく周期が変わる様子が簡単に見える。実験的には正しいのに、小学校では「間違い」とされる。 twitter.com/kikumaco/statu...

タグ: 単振動 小さい 振幅 近似的

posted at 17:40:15

Kouji Ohnishi 大西浩次 @koujiohnishi

20年12月24日

小学校学習指導要領(平成29年)解説を見てみると、P65 の下7行目より「振り子の長さや振れ幅を一定にして・、変える条件と変えない条件を制御しながら実験・。その際、振れ幅を極端に大きくならないように適切な振れ幅で実験をおこなようにする。」と書いてあるが、教科書の指導者には書いていない?

タグ:

posted at 17:40:15

黒木玄 Gen Kuroki @genkuroki

20年12月24日

その問題には私も興味があるので、ツイッターで「周囲の人たち」をよく観察しているのですが、一般に数学的能力が高い人はアホなタイポが結構目立ちます。

研究者であるかどうかは無関係だと思います。

自戒するより、細かなノイズに強くなる方法を積極的に教えて行った方が生産的だと思う。 twitter.com/anairetta/stat...

タグ:

posted at 17:48:20

黒木玄 Gen Kuroki @genkuroki

20年12月24日

教える側が細かなミスを無くすことに膨大なリソースを割くのは非生産的で、みんなでノイズに強い思考法を身につけた方がお得。

細かなミスを無くすことに膨大なリソースが必要であることを理解していれば、細かな誤りの存在を親切に教えてくれる人のありがたみも深まる。

タグ:

posted at 17:48:20

黒木玄 Gen Kuroki @genkuroki

20年12月24日

そして、ノイズに強い思考法の存在が常識化していれば、細かなミスを見つけたときにイライラしてしまう原因がノイズに弱い自分自身の思考法にあるのではないかと考えて成長できる余地が増えるかもしれない。

タグ:

posted at 17:48:21

黒木玄 Gen Kuroki @genkuroki

20年12月24日

数学でよくありがちなのは、具体的な計算で理解できる例を1つも用意せずに本を読み進めようとしているせいで、ちょっとしたノイズですぐに分からなくなってしまうというパターン。

私は学生時代に「例を3つ作りましょう」のように教わりました。

タグ:

posted at 17:51:26

Kouji Ohnishi 大西浩次 @koujiohnishi

20年12月24日

ミスタイプの訂正と追加。小学校学習指導要領(平成29年)解説以降、教科書の指導書などには、振幅を大きくすると周期が伸びると書いていないかな?どうでしょう。なお、以前の文科省の小学校理科の観察,実験の手引き詳細には、めちゃめちゃひどい誤りが書いてある。

www.mext.go.jp/component/a_me...

タグ:

posted at 18:01:03

Kouji Ohnishi 大西浩次 @koujiohnishi

20年12月24日

文部科学省
小学校理科の観察,実験の手引き詳細
「おもりの重さや振れ幅が,おもりの1往復する時間に関係すると予想や仮説をもった児童の中には,実験の誤差をおもりの1往復する時間が変化したととらえる児童がいる。これは,自分の予想や仮説に合うようにデータを処理しているからである。」

誤り!

タグ:

posted at 18:01:03

よ @133cmJS

20年12月24日

PythonよりJuliaを学習したい今日この頃

タグ:

posted at 18:07:16

黒木玄 Gen Kuroki @genkuroki

20年12月24日

私のツイッターでの主要テーマの1つは、初歩的な統計学を理解するために必要な計算例とそのプロット(場合によっては動画)およびコンピュータで実際にそれを行う方法の拡散。

統計学は初歩的な段階ですでに結構混み入った複雑な話になっているので、例が頭に入っていないと理解できるはずがない。

タグ:

posted at 18:12:25

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「尤度の定義を見ても、どうしてこれが『もっともらしさ』なのか理解できない」と思った人は非常に正常に頭が働いていると思います。

「尤度」を「もっともらしさ」と言い換えて理解した気分になった人はすでにあっちの世界に行きかけていると思う。少なくとも論理的思考は諦めている。

タグ: 統計

posted at 18:17:14

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 統計学入門の教科書的解説はそういうパターンが実に多い。

論理的な理解を諦めさせることに繋がるような書き方がほとんど至る所でされている。

例えば、論理的な人なら「正規分布を仮定して有意差を出して本当にいいの?」と当然思うはずで、実際に多くの人がそう思っている。

タグ: 統計

posted at 18:20:21

数学とか語学とか楽しいよね @sasaburo

20年12月24日

1次元アレン-カーン方程式を解きました。-1と1に相が分離します。フェーズフィールド法などで使います。
ja.wikipedia.org/wiki/%E3%82%A2... pic.twitter.com/TdUiAZOmj6

タグ:

posted at 19:23:28

非公開

タグ:

posted at xx:xx:xx

@Ryu_chaaa

20年12月24日

C勢からJulia勢になるか

タグ:

posted at 20:09:24

せんりちゃん @EqfLV06jtQ6zEEz

20年12月24日

将棋の渡辺くん懸賞当選💕
ありがとうございました。
イナ先生、渡辺先生
最高のクリスマスプレゼント🎄🎅🎁✨
嬉しすぎて
今夜は眠れそうにない😆 pic.twitter.com/p20V2AvvVg

タグ:

posted at 20:16:31

Hiroyasu Kamo @kamo_hiroyasu

20年12月24日

ええと、日本の官僚は大学受験生がすべて高校卒業見込み者だと思い込むほど馬鹿なのでしょうか? twitter.com/yuru6654/statu...

タグ:

posted at 20:47:51

積分定数 @sekibunnteisuu

20年12月24日

内部告発で自分の所属している組織の不正を公にするなて、とんでもないこと

という認識なんだろうな。

政府のやることに水を差す学術会議なんてケシカランとい認識なんだろうな。 twitter.com/MegJoytoy/stat...

タグ:

posted at 21:08:00

数学とか語学とか楽しいよね @sasaburo

20年12月24日

2次元アレン・カーン方程式を解きました。だんだんくっきりわかれていきます。きれいです。 pic.twitter.com/okSPoVE57k

タグ:

posted at 21:08:02

数学とか語学とか楽しいよね @sasaburo

20年12月24日

3D版です。にょきにょき生えます。 pic.twitter.com/B1SwjvOK6v

タグ:

posted at 21:08:29

天むす名古屋 Temmus @temmusu_n

20年12月24日

【文科省が定めたテスト】! こんなことをいう人は私が学テや新学テのお為ごかしを批判していると知ったら卒倒するかも。というのは冗談ですが、現場教員がテスト会社に甘く、言いなりになっているようにさえ見える理由の一端がわかる気がします。#超算数 twitter.com/sekibunnteisuu...

タグ: 超算数

posted at 21:15:54

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 ディガンマ函数ψ(s)=(d/ds)log Γ(s)やトリガンマ函数ψ'(s)などはコンピュータで気軽に計算できる函数になっている、という情報を知っていれば、後で統計学の勉強のための数値実験で活かせる可能性が高い。

基本特殊函数一発で計算できるのと、MCMC法などで計算するのでは、手間が全然違う。 twitter.com/genkuroki/stat...

タグ: 統計

posted at 21:34:15

いほやきゅ @K0TSU_K0TSU

20年12月24日

ノコギリガザミのハサミの突起

哺乳類の顎のようになっている、先端の突起は犬歯のように鋭く、根本では臼歯のようになっている、突起だけでなく左右のハサミの形状も違う、ただ挟むだけではなく、使い分けていることが分かる、飼育してどのように使うのかを観察してみたい pic.twitter.com/EMP45qNYZ1

タグ:

posted at 21:44:46

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「尤度原理はBirnbaumによって証明されている」とか、「証明されていると〇〇さんが書いていた」の類のことを言っていた人について、あきれるのは、自分でその「証明」とやらを説明しようとせず(おそらく説明できない)、誰か偉い人に名前を出すだけですましていることです。続く

タグ: 統計

posted at 21:48:02

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 そのBirnbaumの「証明」とやらは、前提がアドホックである疑いや、正確に定式化すると穴があるという疑いを払い切れるものでなかったことは、歴史が証明しています。

すでに私はその証拠として、赤池弘次さんの論文とMayoさんのウェブサイトを紹介しています。

タグ: 統計

posted at 21:48:03

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 「尤度原理はBirnbaumによって証明されていると〇〇さんが書いていた」のように述べてしまった人は、その〇〇さんが証明の詳細をきちんと説明していたか、強い疑いがかけられて続けていることを正直に説明していたかを確認し、自分自身が騙されている可能性を疑うべきです。

タグ: 統計

posted at 21:48:03

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 いずれにせよ「尤度原理はBirnbaumによって証明されている」と口走ってしまった人には、その証明とやらの詳細を自分自身が理解していることを示して、さらにそれが正しいことをギャップのない議論で示す立証責任が生じています。

タグ: 統計

posted at 21:48:04

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 私が、非常に残念だと思いかつ、あきれてしまうのは、自分自身の言葉で論証部分を説明せずに、「偉い人の名前」を挙げてすます恥知らずが堂々と発言していることです。

自分で論証し、自分でコードを書き、自分でプロットし、…としていれば、十分に尊敬可能になるのですが。

タグ: 統計

posted at 21:48:06

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 赤池さんによる尤度原理潰しのための論文情報↓

twitter.com/genkuroki/stat...

タグ: 統計

posted at 21:51:58

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 統計学の哲学が専門のMayoさんのウェブサイトでは、「尤度原理潰し」が結構ホットな話題の1つになっています。

例えば↓

errorstatistics.com/2019/12/30/a-p...

タグ: 統計

posted at 21:54:15

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 赤池さんの論文もMayoさんの議論も、ほとんど「詭弁論証潰し」的な様相を示しています。

普通の科学畑の人は「尤度原理」なんて知らないだろうし、知らなくても困らないし、「尤度原理」の正確な主張を知れば極端な考え方過ぎるという理由でノータイムで相手をするのをやめると思います。

タグ: 統計

posted at 21:57:39

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 〇〇〇氏の本に本に書いてあったとか、〇〇〇〇氏が言っていたの類は、「尤度原理」に真面目に相手をする価値があることの証拠には一切なっていません。

一般に、どんなに高名な人であっても馬鹿げた考え方にはまり込んだり、馬鹿げた考え方をまじめに受け取ったりするものです。

タグ: 統計

posted at 22:01:06

shikihuiku @shikihuiku

20年12月24日

あ。これは確かに酷い。
こちらに詳しく取り上げられてる。
oku.edu.mie-u.ac.jp/~okumura/stat/... twitter.com/genkuroki/stat...

タグ:

posted at 22:10:38

Hal Tasaki @Hal_Tasaki

20年12月24日

@Historyoflife むしろ GOTO など(本来は疫病が収束してからやるべきこと)をやって増やしていたわけですから本当に辛いです。政府がここまで平然と医療関係者に破滅的な負担を強いることは想像できなかった。

タグ:

posted at 22:35:26

黒木玄 Gen Kuroki @genkuroki

20年12月24日

#統計 個人的には、「本や論文に書いてあると言っても正しいことを言っていると認めてもらえず、議論を再構成してギャップを埋めることを要求される」という教育を受けた人達には、「統計学においてもそうすること」を結構期待しています。

偉い人の名前を出しても無力であることを知っているはず。

タグ: 統計

posted at 22:59:35

Tootsie Roll @TootsieRoll2581

20年12月24日

きたよ!きたよ!!
高校生が数理モデルを用いた計算をするようになる時代が!!新時代の幕開けだぁぁぁ!!!

news.yahoo.co.jp/articles/c2f9a...

タグ:

posted at 23:01:18

ごまふあざらし(GomahuAzaras @MathSorcerer

20年12月24日

Linux環境でビルドできました.

タグ:

posted at 23:10:12

非公開

タグ:

posted at xx:xx:xx

Tootsie Roll @TootsieRoll2581

20年12月24日

ってか、割とマジでこの子が羨ましい。私も医学部時代にこの分野と出会っていれば、もっと長い期間この研究ができたのに、とつくづく思う。

タグ:

posted at 23:22:23

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

ごまふあざらし(GomahuAzaras @MathSorcerer

20年12月24日

Julia の master ブランチ,make 通るけれど make install すると適切な so がないからセグフォする(笑)

タグ:

posted at 23:35:09

ごまふあざらし(GomahuAzaras @MathSorcerer

20年12月24日

なるほど.これは面白い.

Pythonで数式・グラフのアニメを描ける「Manim」 qiita.com/Tonoyama/items... #Qiita

タグ: Qiita

posted at 23:37:40

OokuboTact 大久保中二病中年 @OokuboTact

20年12月24日

> 「政府がそう言っている」とはぐらかされたりした例があった。

???

#謎の校則

www.nishinippon.co.jp/item/n/676010/

タグ: 謎の校則

posted at 23:48:59

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

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

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

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

送信中

送信に失敗しました

タグを編集しました