如何根据 USP <1033> 生物测定验证计算变异成分的 REML 估计值?你好,
我在根据 USP <1033> 中提供的数据和结果计算 REML 分析结果时遇到了问题。我无法计算表 4 中的数据并获得表 9 中包含的相同结果。有没有人可以处理这些数据并解释如何使用此示例进行 REML 分析?

谢谢!
高级统计建模实验设计质量与工艺工程回复:如何根据 USP <1033> 生物测定验证计算变异成分的 REML 估计值?欲了解更多信息,请参阅以下来源: 1033 生物测定验证
回复:如何根据 USP <1033> 生物测定验证计算变异成分的 REML 估计值?嗨@Luk_s ,欢迎。我可以重现表 5 和表 6。但不能重现表 9。表 9 中的结果似乎是使用不同数据的示例。这可以通过 Var(error) 的差异看出。在表 5 中,它是 0.000766。在表 9 中,它是 0.0022。如果两者都使用相同的数据,它们应该是相同的。因此,表 9 不是基于表 4 中显示的数据。它是来自其他内容的示例,正如“从此类分析中获得的方差分量估计值在表 9 中呈现”(斜体是我的)这句话所暗示的那样;请注意,它并没有引用表 4 中的数据。并且,正如他们在表 5 下正确指出的那样,“由于每个因素的级别数量较少(2 个级别),因此未包括分析师和媒体批次。”,因此执行表 9 中显示的分析类型是不明智的。
值得一提的是,尽管这是不明智的,但我得到以下结果(表 4,级别 0.5):

回复:如何根据 USP <1033> 生物测定验证计算变异成分的 REML 估计值?你好@MRB3855,非常感谢你对我提出的问题的关注。我确信我在开发模型数据时做错了什么。我也做了检查,我们得到了相同的结果,对此我感到非常高兴。我也处理过以前的数据表,但是,我对这个例子中得到的结果感到惊讶。

谢谢!