Extrait de paramètres optimaux FLEXMIX: univariée données

Question

Je suis en train de reproduire une table dans le livre HMM:

Plus précisément, le modèle k = 2.

J'utilise le code:

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) 

Examiner les prieurs et le journal lik il semble que la syntaxe est en train de faire ce que je veux :). Ma question est de savoir comment extraire le composant lambda (15.777 et 26.840 ci-dessus)? Je ne crois pas qu'ils soient simplement la valeur moyenne des données dans chaque grappe.

Answer 1

exp(parameters(x2)) 

Vous pouvez utiliser parameters pour extraire les coefficients mais les paramètres par défaut du modèle à l'aide d'un lien fonction de journal de sorte que vous devez utiliser exp pour reconvertir à l'échelle d'origine.