0

我应用了这段代码，我一遍又一遍地得到这个错误。

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”

谁能帮我这个？

另一方面，如果有人能想出一种方法来分离两个二项分布，那就太好了。

4

1 回答 1

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在新会话中重新安装并重新运行代码。

于 2018-05-10T22:15:52.937 回答