我正在尝试使用加权样本在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))