Я пытаюсь подогнать модель GMM в R, используя взвешенную выборку, где у меня есть наблюдения, для которых взвешена соответствующая ошибка измерения.
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))