我是高斯混合模型的新手,我正在尝试学习如何在R中使用normalmixEM。我已经阅读了文档,以便了解它的工作原理我试图将混合模型拟合到假数据设置我创建,但我得到以下错误,我不知道为什么:
“错误:外国函数调用中的NA / NaN / Inf(arg 5)”
代码:
S=(t(t(runif(5,0,1))))
Initialprops=S/sum(S)
values=matrix(0L, nrow = 10000, ncol =1)
means=matrix(0L, nrow = 10000, ncol =1)
sigma=matrix(0L, nrow = 10000, ncol =1)
for (i in (1:10000))
{values[i,1]=runif(1,0,1)
means [i,1]=i
sigma[i,1]=(2^(i-1))
}
normalmixEM(values,Initialprops,means,sigma,5)
对于假数据集,我刚使用了一组来自(0,1)均匀分布的10000个随机数。我使用了一组5个归一化随机数作为初始混合比例。对于我使用序列1,2,3,...的方法,对于sigma,我使用了序列1,2,4,8 ......
谢谢!