Bland Altman分析的Origin實現(xiàn)（附帶MATLAB實現(xiàn)）
Bland Altman分析的Origin實現(xiàn)（附帶MATLAB實現(xiàn)）
參考教程：
https://cloud.tencent.com/developer/article/1556951
教程是用Python做的，我比著做了MATLAB和origin兩個版本，MATLAB用于理解Python，origin實現(xiàn)點點點傻瓜式操作。
重點講解origin操作，MATLAB代碼附在最后
step1：導(dǎo)入數(shù)據(jù)，true和pred兩列，都設(shè)置為y軸數(shù)據(jù)
? ?
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?
?
step2：計算兩列數(shù)據(jù)差，在F(x)窗格輸入col(a)-col(b)即可，同理，計算pred和true數(shù)據(jù)的均值，如圖
?
?
step3：使用統(tǒng)計功能對差值進行統(tǒng)計，默認選項即可，將獲得的結(jié)果，新建立兩列進行存儲mean和standard deviation
?
?
?
?
?
step4：計算95的置信區(qū)間，up=col(e)+1.96 * col(f),down=col(e)-1.96 * col(f)
?
?
?
?
step 5:對數(shù)據(jù)進行排序，方便繪圖，選中一列數(shù)據(jù)，然后右鍵，按圖中選擇即可
?
step6: 繪散點圖，選中mean和diff數(shù)據(jù)，plot-》scatter
?
?
?
?
step7：增加置信區(qū)間線和差值平均值，繪制好的散點圖激活，然后insert，plot to layer ,line
?
?
彈出窗口，點ok
?
?
彈出的窗口，添加需要添加的線條，先up然后按住ctrl選中down，點向右箭頭，
?
同理，添加diffmean數(shù)據(jù)，如圖
?
?
點擊ok，得到初步的圖像了
?
?
step8：調(diào)整顏色，字體，線條類型，邊框類型等元素
粗調(diào)如下：
另附MATLAB代碼如下：
%繪制blandaltman圖形
clear all;clc;close all;
ytrue=[123.25, 126.83, 79.70, 129.23, 110.23, 116.53, 88.35, 90.04, 129.32, 211.44, 142.93, 135.41, 135.63, 92.58, 75.82, 93.90]';
ypred=[86.32, 136.65,73.33, 133.19, 119.34, 101.29, 88.40,113.25, 131.44, 210.50, 124.82, 112.81, 139.93, 92.13,77.04,90.14]';
x=(ytrue+ypred)/2;
ydiff=-ytrue+ypred;
ydiffmean=mean(ydiff);
ydiffstd=std(ydiff);
up95=ydiffmean+1.96*ydiffstd;
down95=ydiffmean-1.96*ydiffstd;
figure
scatter(x,ydiff,'filled','r')
hold on?
yline(up95,'-r','95up','LineWidth',2)
hold on?
yline(down95,'-r','down95','LineWidth',2)
hold on?
yline(ydiffmean,'--b','diffmean','LineWidth',2)
?
