Python使用robustats计算medcouple时Jupyter与IntelliJ IDE结果不一致

• 第一步：确认包版本一致性
  即使你确认了使用同一个虚拟环境，仍建议分别在Jupyter Notebook和IntelliJ IDE中执行如下代码，输出numpy和robustats的版本号做比对：

import numpy
import robustats
print(numpy.__version__)
print(robustats.__version__)

若版本不一致，说明IDE的解释器配置仍存在偏差，需要重新在IDE的项目设置中指定对应虚拟环境的Python解释器绝对路径。

• 第二步：先对输入数组排序后再计算
  medcouple的部分开源实现版本要求输入数组为升序排列，未排序的输入可能触发边界计算异常返回nan。你可以调整代码如下再测试：

from robustats import medcouple
import numpy as np

x = np.array([9325.06, 6206.00, 10000.00, 9569.78])
x_sorted = np.sort(x)
print(medcouple(x_sorted))

• 第三步：重装依赖包排除编译异常
  如果版本一致、排序后仍返回nan，可以在对应虚拟环境中执行如下命令强制重装robustats，排除包文件损坏或平台编译异常的问题：

pip uninstall -y robustats
pip install --force-reinstall robustats

• 第四步：检查IDE浮点运算配置
  部分IDE的运行配置默认开启了浮点优化选项，可能影响数值计算结果。你可以在IntelliJ IDE的当前项目运行配置中，关闭Enable fast math类的浮点优化选项，再重新运行代码测试。
• 备选方案：使用官方稳定实现
  如果以上步骤都无法解决问题，可以替换为SciPy官方提供的medcouple实现，兼容性和稳定性更高：

from scipy.stats import medcouple
import numpy as np

x = np.array([9325.06, 6206.00, 10000.00, 9569.78])
print(medcouple(x))

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