Биномиальное Распределение в в R

Question

Я применил этот код, и я получаю эту ошибку снова и снова.

library(mixtools)

simulate <- function(lambda=0.3, mu=c(0, 4), sd=c(1, 1), n.obs=10^5) {
    x1 <- rnorm(n.obs, mu[1], sd[1])
    x2 <- rnorm(n.obs, mu[2], sd[2])    
    return(ifelse(runif(n.obs) < lambda, x1, x2)) 
}

x <- simulate()

model <- normalmixEM(x=x, k=2)

Ошибка: объект 'C_normpost' не найден

Может ли кто-нибудь помочь мне с этим?

Если, с другой стороны, кто-то может придумать способ разделения двух биномиальных распределений, было бы здорово.

Спасибо!

Answer 1

Если у вас правильно установлен пакет mixtools, ошибки не должно быть. У меня установлена ​​версия 1.1.0, которая, по словам CRAN, является самой последней выпущенной версией. .C() вызов C_normpost осуществляется в трех разных местах в теле normalmixEM. После этого кода я вижу:

str(model)   # don't type just `model` ... appears to have no print method for it.
List of 9
 $ x         : num [1:100000] 3.33 3.17 -2.04 3.66 5.11 ...
 $ lambda    : num [1:2] 0.299 0.701
 $ mu        : num [1:2] -0.00253 3.9986
 $ sigma     : num [1:2] 0.994 1.001
 $ loglik    : num -197426
 $ posterior : num [1:100000, 1:2] 1.94e-03 3.66e-03 1.00 5.17e-04 1.42e-06 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : NULL
  .. ..$ : chr [1:2] "comp.1" "comp.2"
 $ all.loglik: num [1:39] -427524 -216959 -208812 -204523 -200531 ...
 $ restarts  : num 0
 $ ft        : chr "normalmixEM"
 - attr(*, "class")= chr "mixEM"

summary(model)
summary of normalmixEM object:
           comp 1   comp 2
lambda  0.2992979 0.700702
mu     -0.0025302 3.998599
sigma   0.9935503 1.001362
loglik at estimate:  -197426.2 

Таким образом, он проделал довольно хорошую работу по оценке средних (0 и 4), а также очень хорошую работу по оценке дисперсий этих компонентов (1 и 1). Возможно, вам следует попытаться переустановить mixtools в новом сеансе и повторно запустить код.