Python 中的 Bland-Altman 图

是否有可能制作一个布兰德-奥特曼情节 http://en.wikipedia.org/wiki/Bland%E2%80%93Altman_plot在Python中？我似乎找不到任何相关信息。

此类情节的另一个名称是Tukey 均差图.

Example:

如果我正确理解了绘图背后的理论，则此代码应该提供基本绘图，而您可以根据自己的特定需求对其进行配置。

import matplotlib.pyplot as plt
import numpy as np

def bland_altman_plot(data1, data2, *args, **kwargs):
    data1     = np.asarray(data1)
    data2     = np.asarray(data2)
    mean      = np.mean([data1, data2], axis=0)
    diff      = data1 - data2                   # Difference between data1 and data2
    md        = np.mean(diff)                   # Mean of the difference
    sd        = np.std(diff, axis=0)            # Standard deviation of the difference

    plt.scatter(mean, diff, *args, **kwargs)
    plt.axhline(md,           color='gray', linestyle='--')
    plt.axhline(md + 1.96*sd, color='gray', linestyle='--')
    plt.axhline(md - 1.96*sd, color='gray', linestyle='--')

中对应的元素data1 and data2用于计算绘制点的坐标。

然后你可以通过运行创建一个图，例如

from numpy.random import random

bland_altman_plot(random(10), random(10))
plt.title('Bland-Altman Plot')
plt.show()