Question

由于对我前一个职位作出了令人 wonder的回答,我使用了下文链接中规定的程序,使我的数据与三个正常分布相适应:

https://stats.stackchange.com/questions/10062/ ,-r- Pack-to-use-to-calculate-component-parailes-for-a-mixture-model>https://stats.stack Exchange.com/questions/10062/ -r- Package-to-use-to-calculate-component-parailes-for-a-mixture-model

在符合我的数据之后,三个正常分发的参数如下:

      pi      mu sigma
1 0.5552 -0.4868 2.044
2 0.2739  8.3846 1.399
3 0.1709 12.5317 1.036

为了检查我的数据(x)和模型分发(ee)之间的协议,我采取了以下步骤:

e1 <- rnorm(5552, mean=-0.4868, sd=2.044)
e2 <- rnorm(2739, mean=8.3846, sd=1.399)
e3 <- rnorm(1709, mean=12.5317, sd=1.036)
ee <- c(e1,e2,e3)
qqplot(x, ee)

I got the qqplot as follows: (https://i.stack.imgur.com/3favy.png

看来不是坏的,因此,我想计算这一模式人口获得价值等于或低于2.0的数值。你们能否想教我如何利用R来计算这一价值?

示范人口的密度图见本文(。

Answer 1

我认为,你只是想用类似东西计算尾矿的概率。

> sum(ee <= 2) / length(ee)
> [1] 0.4936

Answer 2

您的参数载于数据框架,param(而且可能的话,这将与单列的矩阵相配合)。

个人对“2线”的贡献:

> probs <- with(param, pi*pnorm(2, mu, sigma) )
> probs
[1] 4.930888e-01 6.883473e-07 2.409615e-25

总数:

> prob <- with(param, sum(pi*pnorm(2, mu, sigma)) )
> prob
[1] 0.4930895

仅仅看,你会猜测,几乎所有捐款都将来自第一部分,因为另两个捐款的冲动手段远不错。第一笔捐款的价值主要以(批量)估计数为主,因为其大部分数额是“2英尺”。

友情链接