注册后推荐绑定QQ，之后方才可以使用下方的“用QQ帐号登录”。

您需要 登录 才可以下载或查看，没有帐号？立即注册

x

帮临床一个同学做数据分析，才接触到Bland-Altman图的有关知识。在网上看到很多帖子问及到该图制作问题，想把自己收集的一些知识信息和大家分享下。也希望借此抛砖引玉，请有经验的童鞋踊跃跟帖~

/ ?# J+ q! X. Q% v- E2 L2 z+ d: D0 L) A  `2 O8 N% d

一、理论背景

; _5 `& b: p+ FBland-Altman法：此种方法为定量与定性方法的结合，其原理是对于两种评定间的差异进行随机效应分析，来解释说明一致性问题。Bland与Altman认为，测量误差不会影响变量间的相关，但会影响一致性；Bland-Altman的方法实际上是对两种评定之间差异的一种观察，横轴为每个被观察对象评定得到的均值，纵轴为两种评定间的差异值，考察评定者间均数及差异的关系。Bland-Altman图法主要是观察两种测量之间差异的分布。在纵轴上以差异的均值和理论0值(以比率为纵轴的以理论1值)为均值参考线，另外再添加md±1.96SD即差异均值的95%的置信区间的参考线，也称为95%的可接受的一致性界限。报告结果需结合图形做合理解释，一般情况下，报告在一致性界限外的点的百分比，与理论值相差的最大值，结合临床，进行合理解释，决定对新方法的取舍。

4 D2 P8 t; U* d$ O+ M4 y& k5 T3 u3 v3 r* L; q6 S) i

理论知识可参考文献：

2 f' L9 ^5 S$ P% f, h4 S" @5 X, ]3 u+ PBland, J.M., & Altman, D.G. (1986). Statistical methods for assessing agreement between two methods of clinical measurement. Lancet, 327 (8476), 307-310.

4 e: b% J! {1 t陈卉. Bland-Altman分析在临床测量方法一致性评价中的应用. 中国卫生统计,2007,24(3):308-315." o, O" M% S+ k3 \

……

/ D3 o7 z/ n3 G+ H5 {

9 J  j  G( i) e) e- Y- u7 @) B二、B-A图制作

6 I# X- V' S% q4 W, Q$ g9 a, p3 [    因为对方看到本专业发表的很多文献都是用SPSS制作B-A图，所以强烈要求希望也用SPSS实现。我花了点时间做功课，结论是SPSS果然没有Bland-Altman图的专门模块，但稍费点周折还是可以用SPSS做出来的。因为网上SPSS的制图步骤比较少看到，所以这里详细阐述下自己的制图过程。2 X! @0 V4 T, q- z9 ?

(一)如何用SPSS制作Bland-Altman图？

0 s) m  |2 t/ l) H# V方法一：7 X2 }6 a) K& Y

可直接复制下面语法在SPSS中运行(注意更改变量名称，以及根据diff的均数和标准差手动调整Y轴上的分界线)：

" j; m' i; z; P9 n

- V# o) _" R+ a*步骤1：计算两种方法测量值的差值(X1-X2)和两种方法测量值的均数(X1+X2)/2.# o$ l4 z* O$ K9 H7 @6 |( x) ^

COMPUTE diff = IOLOD-OROD .

- w" D; F% U1 m. ], r2 l4 lCOMPUTE mmean = (IOLOD + OROD)/2 .

! i8 d) |) i+ G$ c# [: ^*步骤2：计算差值diff的均数和标准差.$ g) y- \/ ?( ^' \0 d

DESCRIPTIVES VARIABLES=diff5 v4 }, t& E* h. A0 `1 r

/STATISTICS=MEAN STDDEV MIN MAX .2 ?' G5 Z4 W" R

*步骤3：手动计算“均数+1.96*标准差、均数、均数-1.96*标准差”，在作图时需手动更改Y轴对应数值.:: |2 J% g/ v2 W8 S8 C8 F4 G( H

*步骤4：以两种方法测量值的均数(X1+X2)/2为X轴，以两种方法测量值的差值(X1-X2)为Y轴做散点图，同时显示四条分界线(Y=均数+1.96*标准差、Y=均数、Y=均数-1.96*标准差、Y=0).8 J6 e3 L) V: @6 h$ ]# E. A

*Chart Builder." z/ Y7 o2 U  ^( a

GGRAPH% I: c3 {% L8 ^& B/ O$ Y

/GRAPHDATASET NAME="graphdataset" VARIABLES=mmean diff MISSING=LISTWISE REPORTMISSING=NO3 q8 @  y4 Z2 D0 R) |

/GRAPHSPEC SOURCE=INLINE.

+ s6 E5 T/ y+ M" `BEGIN GPL$ i4 m3 Q; s% q* c8 I, |* l

SOURCE: s=userSource(id("graphdataset"))

. j2 [% @/ W6 |9 L4 S6 fDATA: mmean=col(source(s), name("mmean"))

) P3 h& Z! j% |1 ^6 D. aDATA: diff=col(source(s), name("diff"))8 z+ e" i) N/ V8 Q% O

GUIDE: axis(dim(1), label("mmean"))0 O6 j% S, F* L/ z- J' B4 |* J2 X

GUIDE: axis(dim(2), label("diff"))

- R: C4 ~* a! W& Z6 E$ IGUIDE: form.line(position(*,1.036946))

. N: |$ g* G) A, ]7 `3 z3 jGUIDE: form.line(position(*, 0.4659))" _8 L" b  Y, R- o' B. e

GUIDE: form.line(position(*, -0.10515))& `) @: V  {9 a; s# j4 O" x

GUIDE: form.line(position(*, 0))* Z' q- U9 y+ W7 P1 a  X

ELEMENT: point(position(mmean*diff))

& q! a9 F0 u" n- ZEND GPL." W9 L5 P; B. t

*步骤5：做出图后，图中线条均为实线。若想把其中几条改为虚线或其他线条，可双击输出的图，调出图表编辑器，8 B$ Y4 l% n2 ~5 c# ~5 {: r+ P

双击要改变的那条线，可弹出“属性”对话框，从“样式”下拉菜单中选择你想要的线条." `" O7 W8 s, J& j. W: w; G

$ N  ~$ p) e: ?, k6 K3 H4 m

方法二：

# x. l9 A$ p- e% u       需要指出的是，这里我用语法直接显示出了Y轴上的4条分界线，当然也可以在运行完步骤1、步骤2的语法后，直接调用SPSS菜单命令graphs中interactive的scatterplot以均数为横轴，以差值为纵轴做散点图。然后手动添加Y轴4条分界线。具体为：在spss输出窗口中再对散点图进行编辑，双击图形进入编辑状态，点击右键选择axis中的纵轴，点击reference lines界面，添加四条参考线，一条是理论值线，差值应为0(纵轴为比率值时则为1)，另外加上差值的均值线与95%置信区间两条线，关闭编辑窗口即完成Bland-Altman plot的制作。" m2 Q# F/ k% t4 q) V

7 }  _( X0 \8 s$ z

SPSS作图语法参考材料：6 O, c- _3 X& g" b

http://www-01.ibm.com/support/docview.wss?uid=swg21476730* r5 [4 x( e% {& a7 I

http://www.dxy.cn/bbs/topic/18429986?tpg=5&age=0

0 D+ S( l* I' v; E# K- x1 T9 k. f/ l9 J

(二)其他制图软件推荐：

% c; U' k5 ]/ |! E网上很多人写过用SAS、R制作B-A图的帖子，SAS制图相对更简单容易些，而R做出的图效果更美观。, e: f" Y, S/ R" B! A2 |$ L

此外也有推荐含B-A专门制作模块的Medcalc、Graphpad prism v4.03软件。感兴趣的不妨一试。

6 i. o9 O1 z  W/ I* ~3 |

- @& Q; k# W, d* w*sas Bland-Altman plot procedures.

* D5 A8 Z  z  _7 n7 |  [data temp;6 ]# p. [% i5 F/ t. J

input x y@@;" y% |/ t' `, [& @

cards;. u7 f7 a8 y4 _8 h2 x

775.5 51

& [' ^  H8 ~' Y1550.5 61, C* N& V( D! p' m; }

865 50! ~* b1 `8 v, n1 h" K/ B6 `

617.5 35

0 Y+ o* r  M) Y* i1 Q7 B1 @& W$ m, g9 J1045 -10

$ t, A: k3 C/ o. m" X;% Y/ t. L8 k9 V# V

proc gplot;, d% H' k  E+ }4 `" Y' C4 s9 n

plot y*x='0'/vref=0 a b;/*a和b用置信区间上下界代替*/

7 y; O1 z. b" @, Z% rrun;

, g3 M; V2 r) w8 N; \# j4 P  K3 k具体事例及用法可参考以下资料：

; g; k) s- d# t( z$ A7 ^3 rBland JM, Altman DG. Measuring agreement in method comparison studies.Statistical Methods in Medical Research 1999; 8: 135-160