本文介绍了如何在R语言中利用sample函数生成基于概率的二元混合正态分布样本，并通过nlminb函数进行参数估计。在估计过程中，详细讨论了初始值的选择和参数约束，以确保计算的准确性。
- comp <- sample(c(0, 1), size = 100, prob = c(0.7, 0.3), replace = T)
- rnorm(100, mean = ifelse(comp == 0, 0, 1), sd = ifelse(comp == 0, 1, 2))
猜测分布是两个正态分布的混合，需要估计出函数中的5个参数：p、μ1、μ2、σ1、σ2。
在R中编写对数似然函数时，5个参数都存放在向量para中，由于nlminb()是计算极小值的，因此函数function中最后返回的是对数似然函数的相反数。
|
1
2
3
4
5
6
7
8
|
> l1=
function
(para)
|
最低0.47元/天 解锁文章
扫一扫
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
