混合ガンマ分布モデルをデータにフィットさせる
Rのmixtoolsパッケージでは，混合ガンマ分布モデルもデータへフィットさせることができる。 #λ= .25,.75，α - 4, 10，β= 1, .2が正解 set.seed(0) dat<-c(rgamma(100,4,1),rgamma(300,10,.2)) hist(dat,breaks=20, main="",col="lightblue") model<-gammamixEM(dat) model$gamma.pars model$lambda まあまあ，て感じ。これ…
