こんにちは、えんせきです。
今日は技術的な要素はありません。いつも読んでいただきありがとうございます。
つまりなにしたの？
半年以上Blogを2日に1回投稿をほそぼそと続けた。
完璧に出来ないくせに完璧主義者で完璧にできないとすぐに飽きてしまう自分が
何でBlogを30週間も続けられたのか覚えておけば別の習慣を始めなきゃいけない時に参考になるかと思って書いた。
基本スタンスは心理的な障壁（完璧のレベル）を下げること。
その試行錯誤をこのBlogを例に紹介する。
Python+NumpyでFIML（完全情報最尤推定法）試してたら色々勘違いをしていたっぽい話
こんにちは、えんせきです。
タリーズブラジルコレクション ブラジル ファゼンダ プラナウトがうまくてハマりました。
いただきものの電動ミルがあるので、豆を買ってきて家で挽いて飲んでます。
*1
つまりなにしたの？
欠損値の補完方法を最近色々試していて、FIML（完全情報最尤推定法）という手法に行き着いた。
FIMLは欠損値を"補完"するわけじゃなくて、欠損値が含んでいても統計量を推定できるってことだったっぽい。
【追記: 2018/03/12】
正規化したら文献中MPlusで推定された結果と同じ値になりました。
この辺、今でもちょっと疑っているので詳しい文献があったら読もうと思います。また訂正いれるかもしれません。
前回はこちら
ensekitt.hatenablog.com
