廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

2018-10-20から1日間の記事一覧

EMアルゴリズム:ランダムな欠測があるデータから多変量正規分布のパラメータ推定

R

http://ebsa.ism.ac.jp/ebooks/sites/default/files/ebook/1881/pdf/vol3_ch9.pdf ↑アルゴリズムはこれに出ています。R のコードだけ貼ります。 library(mvtnorm) EM_mGauss <- function(y,muini,sigmaini) { for(j in 1:1000){ z <- y m <-is.na(y) N<-nrow…