为什么“效应检验”下的平方和不等于整体模型的平方和？大家好，
我使用标准最小二乘法拟合了一个包含一个连续变量和一个名义变量的简单模型。现在我想知道为什么方差分析和效应检验下的平方和似乎不匹配。我使用的是 JMP 18.0.1。
3.16 和 0.01 并未给出整体模型的 SS 3.37。好像缺少了 0.2。有人能解释一下这是为什么吗？
这让我特别困惑，因为我正在查看同事提供的 minitab 中相同数据的分析，其中因子“月份”确实获得了 0.209 的增加 SS。所有其他东西，如 DF 和 SS 都匹配。最后，这对 p 值判断“月份”对拟合的重要性有影响，这就是它一开始困扰我的原因。
问候，
康斯坦丁
基础数据分析与建模回复：为什么“效应检验”下的平方和不等于整体模型的平方和？如果您附上实际数据表而不是仅附上输出，这将非常有助于查看您的结果。有点像在学校里“展示您的工作”。
在您与 Minitab 的比较中，我可以立即得出的一个结论是 Minitab 报告的是顺序 SS（类型 1），而不是部分 SS（类型 3）。顺序 SS 依赖于顺序，而部分 SS 对模型顺序更稳健（并且是首选，通常是默认的）。JMP 将使用类型 3 作为默认设置，并使用 REML（我认为 Minitab 仍然在 EMS 估计平台上）。
我不确定你的因子是什么，哪些是连续的，哪些是名义的。批次和月份似乎都不连续？你有 9 个批次和 2 个月。但我不知道你的 MSE DF 来自哪里？
如果您查看参数估计表并添加 VIF 列，这些值是什么？
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？非常感谢，现在清楚多了！
我已将数据附加到此回复。VIF 相当小。
您的观察对我很有启发。昨天我实际上想知道模型的顺序是否起作用，并切换了月份和批次，但 JMP 并不关心。您对类型 3 的解释完全合理。我现在还明白，JMP 中的效应会根据其他效应的存在进行调整，而 minitab 不会这样做，这解释了差异。从模型中删除“批次”后，JMP 给出“月份”的 SS 为 0.2，这是 minitab 从包含两个因子的模型中输出的结果。
我相信这个问题已经得到解决。
但是，我也和你一样对 DF 感到疑惑。月份是一个连续因子。那么为什么最终 DF 为 1？也许你看完数据后就会明白。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？你好@Bashburz1 ：“但是，我也和你一样对 DF 感到疑惑。月份是一个连续因子。那么为什么它的 DF 最终为 1？”
这是因为作为一个连续因子，它只向模型引入了一个参数：“斜率”，即与月份相乘的项。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？知道了，谢谢！
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？顺便说一句，如果你将月份的数据类型更改为名义型，你将获得 9 个月份 DF。我不确定你为什么将月份视为连续变量？
此外，如果您单击“拟合模型”输出中响应旁边的红色三角形，选择“估计”>“顺序检验”，您将获得类型 1 SS。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？嗨@statman我的猜测是“月”确实是连续的；它实际上是时间……但它的单位是月（可能是秒、分钟、天、周等）。我猜这是制药行业的稳定性分析（保质期）。我说得对吗@Bashburz1 ？如果是这样，“月”就是连续的。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？嗯，很有趣，如果你关心的是稳定性，我不会用线性模型来建模月份。我会使用抽样思路和控制图来评估稳定性。如果月份很重要，那意味着什么？是否有一个不寻常的月份（用戴明的术语来说很特别）还是每个月都有很大的变化？
您当然不能完成研究并说月份很重要，1 月份更好，让我们就此打住。现在，如果您在烘焙过程中问时间是否有影响，您肯定会认为“时间”是一个连续因素。
但不可否认的是，我并不是一家“保质期”的中小企业。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？@statman适当的类比是您的烘焙过程；“稳定性”的线性模型（这是行业中通常所说的保质期……不要与 QC 术语中的“稳定过程”混淆）是以时间为协变量的 ANCOVA 模型。
可能比您要求的还多，但是您可以继续！
https://database.ich.org/sites/default/files/Q1E_Guideline.pdf#page13
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？感谢您今天早上的轻松阅读。哇，这显然是委员会编写的指导方针。
“稳定性研究的目的是建立基于
在对至少三批药物物质或产品进行测试后，重新测试
适用于所有未来批次的保质期或保质期以及标签储存说明
在类似情况下制造和包装”。
3 从何而来？3 个批次代表所有批次……在哪个世界？
正如您所说，稳定性一词的使用让我感到困惑，文件中也没有定义。它似乎是指药物疗效的下降？
我喜欢这个部分：
“数据不适合进行统计分析”
无论如何，对于 OP，如果您遵循本指导，月份应该是协变量，您应该对协变量使用顺序测试，然后对研究中的其余因素使用部分 SS（类型 3）。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？是的，“稳定性”是指监测药物任何属性随时间的任何下降或变化。感谢您提供有关使用 JMP 中的顺序测试的提示。在这种情况下，您通常有三个批次（可能是出于实际原因，历史上是这样），并且想知道它们是否表现相同。因此，您从批次、时间和它们的交叉项开始。您首先测试是否可以去掉交叉项，这意味着所有斜率都相等。然后，您想去掉“批次”，即显示截距相等。您将所有数据点放在一个图中并拟合一条线。有时您可以有更多因素，理想情况下您希望它们没有影响。因此，月份/时间通常是这里的最后一位。
我在 JMP 网站上看到了一门关于混合模型及其如何应用于此主题的精彩课程，但它附带一个免责声明，即这种方法根本不符合所提到的指导方针，所以我们开始吧:)
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？没错，月只是时间单位。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？你好@Bashburz1 ：不确定您是否知道，JMP 中有一个用于稳定性分析的模块。
https://www.jmp.com/support/help/en/18.1/index.shtml#page/jmp/stability-analysis-in-the-degradation-platform.shtml
它的功能与拟合模型平台相同，但它有一些特定于 ICH Q1E 稳定性分析的选项。
回复：为什么“效应检验”下的平方和不等于整体模型的平方和？我不知道，谢谢！理解这些步骤并手动执行仍然很好，但该模块实际上是自己完成这项工作的。
