参考区间分析目的参考区间referenceintervalRI由参考样本群体的测量结果界定用于解释个体结果的正常范围对应 CLSI EP28-A3c。ivdtools提供三种方法百分位数法nonparametric、参数法parametric正态或对数正态、稳健法robust。三种方法对样本量、分布形态和离群值的敏感性不同建议同时对同一数据应用多种方法进行比较并在建立参考区间前先做异常值检测与正态性评估。本文演示用outliers_test()与normal_test()做预先分析再调用reference_interval(method all)同时输出三种方法的参考区间。示例数据为确定性的教学数据n120轻度右偏并含离群值。::: celllibrary(ivdtools) library(readr):::函数概述函数 主要用途 关键输入或输出outliers_test()异常值检测 Grubbs / ESD / Dixon / IQRnormal_test()正态性检验与图形methodauto按样本量选检验reference_interval()三种方法参考区间methodall、interval、ci、idplot.reference_interval()分方法绘制区间 jitter 上下限示例三法并用的参考区间读取和核验数据::: cellri_data - read_csv(./data/reference-interval.csv, show_col_types FALSE) ri_data - as.data.frame(ri_data) str(ri_data):::预先分析一异常值检测EP28 建议在建立参考区间前审查离群值。样本量 n120不适合 DixonDixon仅适用于 n≤30 使用 Grubbs 与广义 ESD::: celloutliers_test(ri_data, col value, method grubbs):::::: celloutliers_test(ri_data, col value, method esd, r 3):::参数说明outliers_test()的method可选grubbs、esd、dixon、iqr ESD 通过r指定最多检出数Dixon 通过type指定尾部both/min/maxIQR 通过coef指定倍数默认1.5。检测到离群值不代表应删除------应回查原始记录、实验过程与方案规定在有充分理由时进行包含/不包含的敏感性分析。预先分析二正态性检验::: cellnt - normal_test(ri_data, col value, method auto) nt:::::: cellplot(nt, type c(qq, his)):::methodauto对 n120 自动选择 Anderson-Darling 检验p值很小提示数据不服从正态分布 该数据右偏并含离群值。QQ图与直方图可直观确认偏态。三法并用的参考区间::: cellri - reference_interval( ri_data, col value, interval 0.95, ci 0.95, method all, id id ) print(ri):::::: cellplot(ri):::参数说明reference_interval()中interval是参考区间覆盖度默认 0.95对应 2.5%–97.5%ci是区间端点的置信水平两者含义不同。method可取percentile、parametric、robust、all或其中任意组合。id用于重复样本检测。本示例三种方法结果比较百分位数法73.4–157.6非参数、不作分布假设但上限被两个大离群值抬高参数法73.3–145.9自动检测到数据非正态后做了 log变换transformedTRUE再回代稳健法72.2–134.7基于 Huber M 估计与 MAD对离群值最稳健。百分位法上限明显高于稳健法正是离群值影响的表现。参考区间应在样本来源、样本量、分布特征与研究方案指导下选择报告方法不宜只挑结果好看的一种。结果判读要点方案先定明确参考人群、纳入排除标准、样本量百分位法 EP28 建议≥120与分区策略。先检后算先做异常值检测与正态性评估记录离群值处理决定是否排除、依据。方法可比三法同出便于发现离群值或分布假设对区间的影响报告时说明所选方法与理由。参数区分区分interval覆盖度与ci端点置信度不要混用。转换透明参数法自动 log 变换时要报告变换类型与回代方式。非唯一结论参考区间不等同于临床决定限也不自动构成可接受性判定。