Bland-Altman一致性评估学习

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

最近一段时间一直在研究Bland-Altman作图的SAS实现方法，阅读了一些文献和SAS help，拿出来和大家分享一下。
Bland-Altman方法多用于评价测试方法与金标准方法的一致性情况，主要用于连续性变量的测定，而对于分类变量，则多采用敏感度、特异度、准确度、Kappa等参数进行评估。
其基本思路为：利用原始数据的均值和两种检测的差值（或比值），分别以均值为横轴，差值为纵轴做散点图，计算差值的均数和差值的95%一致性界限（均数±2SD，也有文献使用均数±1.96SD）。
在SAS中，没有相应的程序直接绘制Bland-Altman图，有大神曾编制了Macro宏程序来实现，相应的PDF说明我贴在了附件里。该方法非传统的Bland-Altman图，有需要的道友可自行尝试。 BlandAltman.zip (544.39 KB, 需要: 2 个论坛币) 本附件包括：

• ba.sas7bdat
• BlandA_Fernandez 2009 WUSS.pdf
• BlandAltman.Mac
• BlandAltman_maccall.sas

2023-5-6 22:38:27 上传

宏程序
需要: 2 个论坛币  [购买]

更一般的方法为使用PROC SGPLOT程序来绘制该图像，具体代码见楼下。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Bland-Altman Altman Bland blan Land

1. data diffs;
   
2. diff = wn12-pn12; #此处计算两连续变量差值
   
3. mean = (wn12+pn12)/2; #此处计算两连续变量均值
   
4. run;
   
5. proc sql noprint; #获得95%置信区间数值
   
6. select mean(diff) - 2*std(diff), mean(diff)+2*std(diff)
   
7. into :lower, :upper
   
8. from diffs;
   
9. quit;
   
10. proc sgplot data = diffs; #绘制B-A图，包括散点、上下95%区间线、参考线y=0及loess拟合曲线
    
11. scatter x = mean y = diff / MARKERATTRS = (color = black symbol = CircleFilled);
    
12. refline &upper &lower / lineattrs = (pattern = 2);
    
13. refline 0 / lineattrs = (thickness = 2);
    
14. loess x = mean y = diff / CLM CLMTRANSPARENCY = 0.7;
    
15. title 'Bland-Altman Plot';
    
16. xaxis label = "Mean of Measurements";
    
17. yaxis label = "Difference between Devices";
    
18. run;

复制代码