在本例中,您拟合两个混合分布,然后标识属于第二个混合分布的其中一个聚类的观测。
1. 选择帮助 > 样本数据库,然后打开 Reliability/Mixture Demo.jmp。
2. 选择分析 > 可靠性和生存 > 寿命分布。
3. 选择 Y1 并点击 Y,事件时间。
4. 点击确定。
5. 点击“寿命分布”红色小三角并选择拟合混合模型。
6. 在 Weibull 旁边的数量框中键入 2。
7. 在“起始值方法”面板中选择可分隔的聚类。
8. 点击执行。
图 3.20 拟合“Weibull(2)”的混合模型
JMP 拟合包含两个 Weibull 成分的混合模型。“对应部分 1”的估计值为 0.231688,指示大约有 23% 的观测服从 Weibull 分布,其 alpha = 9.483153,beta = 3.001963。剩余的 77% 估计来自第二个 Weibull 分布。
为了将该模型与另一个模型进行比较,您可以更改“成分”选择和成分“数量”。
9. 在对数正态旁边键入 1,并且在 Weibull 旁边键入 1。
10. 点击执行。
图 3.21 拟合“对数正态(1),Weibull(1)”的混合模型
“叠加”图随之更新以显示两个混合模型。“模型列表”中的图和统计量指出“对数正态(1),Weibull(1)”混合模型看起来提供了与“Weibull(2)”混合模型极为相似的拟合。
1. 点击“对数正态(1),Weibull(1)”旁边的红色小三角,然后选择保存预测。
将向数据表添加两列:
‒ 对数正态(1),Weibull(1) - 以下成分的预测概率: 对数正态
‒ 对数正态(1),Weibull(1) - 以下成分的预测概率: Weibull
2. 选择分析 > 分布。
3. 从“选择列”列表中选择这两个新列,然后点击 Y,列。
4. 选中仅显示直方图。
5. 点击确定。
6. 在对数正态(1),Weibull(1) - 以下成分的预测概率: Weibull 的直方图中,点击与接近 1 的值对应的直条。
图 3.22 混合概率的直方图
在数据表中,选定了 297 个对应的行。这些是可能来自参数 alpha = 29.90 且 beta = 10.41 的 Weibull 分布的观测。