Bland - Альтман сюжет в Python

можно сделать Bland-Altman plot в Python? Кажется, я ничего не могу найти.

другое название для этого типа сюжета -Tukey mean-разница сюжет.

пример:

enter image description here

9
matplotlib plot python

2 ответов

если я правильно понял теорию сюжета, этот код должен обеспечить базовое построение, в то время как вы можете настроить его под свои собственные потребности.

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()

Bland-Altman Plot

20

может быть, я что-то упускаю, но это кажется довольно просто:

from numpy.random import random
import matplotlib.pyplot as plt

x = random(25)
y = random(25)

plt.title("FooBar")
plt.scatter(x,y)
plt.axhline(y=0.5,linestyle='--')
plt.show()

здесь я просто создаю случайные данные между 0 и 1, и я случайно помещаю горизонтальную линию в y=0.5, но вы можете поместить столько, сколько хотите, где хотите.

0