python - Python 中的 Bland-Altman 图

Question

是否可以制作一个Bland-Altman plot在 Python 中？我似乎找不到任何相关信息。

此类图的另一个名称是 Tukey 均值差图。

示例:

Answer 1

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

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 和data2 中的相应元素用于计算绘制点的坐标。

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

from numpy.random import random

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