python三个变量互换值_连续变量测量方法一致性评价BlandAltman：不同测量值的互换公式...

临床研究常见问题：不同实验室/仪器/测量员检测的某项指标，测量的差异可能对结果有影响，如何处理呢？

• 我自己不说，审稿人怎么知道？心存侥幸，不可取。审稿人通常是同行，你遇到的问题，他很可能早就遇到过。

• 埋头苦看、勇往直前。我把所有患者用不同仪器全部都测一遍。通常时间、财力、患者配合等问题很难开展。

学学高手是怎么做的？定量性的测量指标一致性评价有专门的统计学方法：Bland-Altman。可以评价不同测量方法的系统性差异、测量误差大小、构建不同测量方法的可互换函数、评估转换结果准确性等。

Bland-Altman 方法的临床研究应用：得出不同测量方法的互换公式，用于不同方法测量值的互换。例如将A方法的测量值转化为B方法的测量值。原理很简单，10元人民币和1美元直接比较不了，用汇率公式转换后1美元等于人民币6元，这样就可以比较了。

数据结构

每个研究对象一行，都测量了两种方法。例如：脂肪含量测量方法有两种：甘油三酯酶解法(trig)和Gerber法(gerber)。用两种方法检测45例研究对象，得出不同测量方法的互换公式。下载练习数据methcomp.xls：链接: https://pan.baidu.com/s/1KzWzoWt3nsoAgsp4zqygSg 提取码: tfxm

软件操作

数据分析 - 诊断试验与预测方程 - 定量测量方法比较。“方法1”设置为trig法(TRIG)。“方法2”设置为Gerber法(GERBER)。此处选择重复方法为“1：一次性测量”，查看结果。

 

得到一系列的图表结果：

 

结果解读

本例测量方法互换公式为：

TRIG = 0.07794 + 0.97212 * GERBER

GERBER = -0.08017 + 1.02868 * TRIG

Bland-Altman的图横坐标是两种测量方法的均数，纵坐标是两种测量方法的差值。如果差值=0，表明测量值完全一致。结合临床意义和差值的95%CI可以判断测量值的一致性。

实例

北京大学第一医院的范博士2016年发表在Hypertension 的血压与肾功能关系的研究。Y是肾功能，关键指标是血肌酐。对3153例研究对象在基线和随访时测量了血肌酐。是在不同单位用不同方法检测的：

• 基线的肌酐是在解放军总医院应用血清酶法检测的。

• 随访的肌酐是在北京医科大学第一医院应用血清苦味酸法检测的。

随访和基线的肌酐差异，可能是不同单位测量方法导致的。为确保测量的一致性。研究者选了123例患者在第三方单位(南方医科大学南方医院)进行测量，得到不同测量方法的转换公式，将不同方法的测量值都标化到一个标准。

具体做法：挑选了100个患者，挑选原则是肌酐测量值十等分后，每个分组里面有10个患者。另挑选了23个极端值患者，其中8位的肌酐 <40 μmol/L，15位肌酐 >120 μmol/L。这样做的目的是尽量使研究对象的测量值范围接近总人群。

数据类型

(1) 测的是连续变量。

(2) 两种方法都测了。具体包括三种情况：

① 一次性测量：每个受试者两种方法各测量一次。

② 可互换的重复性测量：每个受试者两种方法各测量2次或以上。

③ 配对的重复性测量：如每个受试者由不同的测量员或在不同的时间分别用两种方法测量一次，每个测量员或每个时点，都有两种方法的测量结果。这样测量员或时点就是配对因素。

(3) 如果从两种方法测量值散点图及方差分析得出D的方差不固定(test for constant variance)，则需要尝试把原始数据进行转换后分析。易侕软件本模块有“数据转换”设置菜单，可以选择对数转换等各种方法。

参考资料

http://www.empowerstats.com/analysis/menuHelp.php?ch=1&module=538

Martin Bland, J. and D. Altman, STATISTICAL METHODS FOR ASSESSING AGREEMENT BETWEEN TWO METHODS OF CLINICAL MEASUREMENT. The Lancet, 1986.

Fan, F., et al., Noninvasive Central Systolic Blood Pressure Is More Strongly Related to Kidney Function Decline Than Peripheral Systolic Blood Pressure in a Chinese Community-Based Population. Hypertension, 2016.