1

私は本HMMのテーブルを複製しようとしています:

ここに画像の説明を入力してください

具体的には、k=2モデルです。

私はコードを使用しています:

dat<-as.data.frame(c(13,  8,  23,   22, 18, 15, 14, 11, 24, 18, 14, 16, 8,  14, 27, 15, 10, 13, 10, 23, 41, 20, 15, 15, 16, 18, 31, 15, 8,  16, 26, 17, 27, 22, 15, 11, 32, 19, 35, 19, 6,  11, 27, 20, 26, 16, 11, 18, 22, 28, 30, 8,  32, 19, 36, 27, 7,  36, 13, 39, 29, 18, 24, 26, 21, 23, 16, 22, 13, 17, 20, 13, 23, 14, 22, 16, 12, 22, 22, 17, 21, 13, 18, 24, 19, 21, 20, 25, 21, 15, 25, 15, 21, 22, 34, 16, 16, 21, 26, 10, 18, 12, 14, 21, 15, 15, 18))
colnames(dat)<-"x"

x2 <- flexmix(x~1, data=dat, k=2,model=FLXMRglm(formula=x~1,family="poisson"))
summary(x2)

事前確率とログの好みを調べると、構文は私が望むことを実行しているように見えます:)。私の質問は、コンポーネントがラムダ(上記の15.777および26.840)を意味するように抽出する方法ですか?それらが単に各クラスターのデータの平均値であるとは思いません。

4

1 に答える 1

4
exp(parameters(x2))

を使用parametersして係数を抽出できますが、モデルはデフォルトでログリンク関数を使用するため、を使用expして元のスケールに戻す必要があります。

于 2013-03-06T23:31:57.143 に答える