R中的Mclust软件包-加权数据样本不起作用?

时间:2018-08-31 14:40:08

标签: r weighted mixture-model gmm mclust

我正在尝试使用加权样本在R中拟合GMM模型,在这里我有观测值被加权了相应的测量误差。

fit.mm <- Mclust(dens,modelNames="V") 

fitnew <- do.call("me.weighted",c(list(weights=w1),fit.mm))

但是当我这样做时,出现以下错误:

Error in eval(expr, envir, enclos) : could not find function "mstepX"

有人知道为什么会这样吗?

这是我正在关注的示例的链接。

https://www.researchgate.net/publication/266350419_Using_Weights_in_mclust

感谢您的帮助。

编辑:

library(mclust)

dens<- c(2.12, 2.71, 3.44, 2.76, 2.72, 0.96, 2, 3.26, 2.5, 1.2, 1.62, 
1.96, 2.6, 1.3, 2.67, 4.4, 1.8, 4.9, 2.39, 1.62, 1.47, 0.89, 
2.52, 1.21, 0.9, 0.8)

err<- c(0.04, 0.11, 0.12, 1.2, 0.12, 0.3, 0.6, 0.6, 0.3, 0.4, 0.3, 
0.34, 0.5, 0.2, 0.03, 2.1, 0.8, 3.9, 0.9, 1.05, 0.95, 0.13, 0.3, 
0.25, 0.1, 0.15)

w <- 1/(err^2)
w1< - w/sum(w) 

fit.mm <- Mclust(dens,modelNames="V") 

fitnew <- do.call("me.weighted",c(list(weights=w),fit.mm))

0 个答案:

没有答案