Python实现Scheirer-Ray-Hare检验与rcompanion结果不一致原因排查

基础信息 #

• 所用Python实现：GitHub公开仓库提供的Scheirer-Ray-Hare检验Jupyter Notebook代码
• 分析数据集：包含2个自变量（单个自变量均设置2个水平）、1个因变量测量值，数据格式已严格匹配仓库提供的示例数据格式要求
• 对照实现：采用R语言rcompanion包运行完全相同的检验流程，对应R端分析代码已公开可访问，变量设置、数据处理逻辑与Python端完全对齐

已观测现象 #

两套实现输出的最终检验结果存在明显差异，经逐行打印Python代码每一步中间计算值交叉核对，确认差异最早出现在平方和计算环节：rcompanion包输出的平方和结果，与Python代码第8个代码单元输出的平方和数值不相等。

已完成排查操作 #

• 逐行运行Python端实现代码，逐环节输出中间计算结果核对，暂未定位差异产生的具体原因
• 遍历调整rcompanion包对应检验函数的各类可选参数，始终无法得到与Python实现一致的收敛结果

待解决诉求 #

• 定位Python端Scheirer-Ray-Hare检验公开实现中的具体代码逻辑问题
• 确认结果差异是否由对rcompanion包函数参数的理解偏差、设置错误导致

内容的提问来源于stack exchange，提问作者lmnopqk