对青少年脊柱侧凸神经影像研究中多重比较校正与纤维束成像指标质疑的回复与澄清

《Global Spine Journal》：Response to Letter to the Editor Regarding “Outcome Selection before Multiplicity Correction May Overstate Neuroimaging Evidence in Adolescent Scoliosis”

　　

尊敬的编辑：

我们怀着极大兴趣阅读了标题为"在多重比较校正前选择结果可能会夸大青少年脊柱侧凸神经影像证据"的读者来信。该来信针对我们最近发表的一项研究¹提出了意见。我们欢迎这一建设性批评，并感谢有机会作出回应。我们希望以下说明能够充分解答并澄清来信中提出的各项观点与疑虑。

来信作者正确地指出，统计分析部分描述了对在ANCOVA步骤中"确认为显著的区域"应用Benjamini-Hochberg（BH）校正程序²，并且从字面意义理解，这样的步骤将保留所有已选定的区域。我们承认，方法与结果部分对分析顺序的描述不够精确。在同行评审过程中，应评审者的要求，我们对全部110个区域皮层厚度测量值的筛选p值应用了BH校正，而我们在原始文章表4¹中报告的11个区域即为在该步骤中被保留的区域。我们亦同意，随后在同一参与者中进行的ANCOVA并不构成独立的确认³。该步骤的作用是验证组间效应并非由年龄或颅内容积所致，"确认性"一词应在此有限意义下理解。Holm校正仅应用于各区域内的配对比较，我们在原始文章表4和表5¹中展示的p值为未经校正的ANCOVA p值。

为直接回应来信作者的核心问题，我们将相同的ANCOVA模型（组别、年龄和颅内容积）应用于全部110个皮层区域，未进行事先筛选，并对完整家族实施了错误发现率（FDR）控制。17个区域显示出名义上的组间效应（p<.05），全部表现为Schroth组的皮层比传统运动组更厚。在5% FDR水平上仍有6个区域保持显著：双侧内侧额叶和枕叶皮层、右侧中央后回以及左侧内侧中央后回（q=.0001–.043）。我们在原始文章表4¹中报告的11个区域恰好即为q<.10的区域（最大q=.084）。在这11个区域中，校正后Schroth组与传统组的差异范围为0.19至0.38 mm，95%置信区间均不包含零；而Schroth组与对照组的差异范围为?0.04至0.00 mm，置信区间包含零。后者表示未检测到差异，而非等效性——正如文章中所指出的，我们并未正式检验等效性。完整结果，包括未经校正的p值、q值以及带置信区间的校正后对比值，可向通讯作者索取。关于原始文章表5¹中的纤维束成像指标，ANCOVA模型是在所调查纤维束特定变量的家族范围内进行评估的（双侧皮质脊髓束和内侧丘系指标），所报告的p值反映了这些特定通路的协变量校正后组间效应。

关于各向异性分数（FA）的意见，我们感谢有机会澄清因讨论部分一处不恰当的措辞而引起的误解。具体而言，"我们的确认性结果主要基于弥散张量指标（FA）……"一句中，括号内的缩写"(FA)"系误加；我们本意是指广义上的弥散衍生纤维束成像指标。重要的是，结果部分及原始文章表5¹中报告的统计发现清楚地表明，右侧皮质脊髓束中观察到的显著组间差异严格涉及纤维数量（ANCOVA p=.018）和纤维比例（ANCOVA p=.019），而FA被明确报告为无显著性（p=.331）。因此，FA从未被用作我们确认性发现的统计依据。尽管我们承认这一括号表述属于校对阶段的疏忽，需要加以澄清，但它并不改变或削弱实证发现，这些发现均已透明、准确地记录在案。

ORCID iDs

Ahmet Payas https://orcid.org/0000-0002-1629-9794

Mustafa Ar?k https://orcid.org/0000-0001-5596-9175