大家好,
我已经使用模型驱动的多元控制图有一段时间了,当我开始更深入地研究该平台的一些统计细节时,我只是刷新了我对这个主题的记忆/技能。在此过程中,我想尝试根据 JMP 在线帮助中的定义方式手动计算 T^2 控制上限( 此处) ;请注意,UCL 有两种定义:历史数据和当前数据。


两者相似,除了 Beta 或 F 分布外,仅取决于行数 n、主成分数 k 和 1-alpha。默认情况下,alpha 为 0.05,因此 1-alpha = 0.95。
作为示例,我使用的是 Steam Turbine Historical.jmp 文件 - 事实上,您可以在此处的在线帮助中找到相同的示例。当您运行 MDMVCC 平台并选择所有列(在本例中没有历史行数据,您使用所有行)时,T^2 图表明 UCL 约为 7。当我尝试使用以下命令手动计算 UCL 时列公式,我得到非常不同的数字。对于 Beta 发行版本,我得到 26.04,对于 F 发行版本,我得到 1.91。
我猜测原因是 JMP 帮助中定义/引用 UCL 中的分布的方式与在列公式或 JSL 代码中调用/使用函数的方式有所不同。以下是我定义列公式的方式。与上述定义进行比较。
贝塔
F
根据脚本索引,这些函数的参数听起来与在线帮助中给出的定义略有不同。我什至尝试交换分布函数中的最后两个参数,但 UCL 仍然不匹配。


毫不奇怪,运行具有 3 个分量的 PCA 分析会产生与 MDMVCC 平台所示相同的中值和 UCL。与 DModX 相同。希望有人可以帮助解释我遇到的差异及其原因。我希望能够协调手动获取的 UCL 值与 MDMVCC 和 PCA 平台生成的 UCL 值之间的差异。
谢谢!,
DS
你好@SDF1 ,
看来您需要使用 Beta Quantile 函数而不是 Beta Distribution 函数。
F 分布函数也是如此。您需要使用 F 分位数函数。
希望能帮助到你。
雄一
嗨@yuichi_katsumur ,
感谢您的清晰解释。从在线帮助中的统计细节来看,并不明显需要使用相关分布的分位数函数。但是,现在您已经解决了这个问题,我可以验证它在这两种情况下都可以手动计算,无论是仅历史数据还是历史加新数据。再次感谢您澄清我的错误!
谢谢!,
DS