4.13 多个单因素两水平设计定量资料Meta分析

4.13.1 问题与数据

【例4-18】为了应用Meta分析确定国内成人哮喘患者血清总IgE含量的平均范围，以便临床哮喘诊断、分型及预后估计时参考，宋士德等人通过计算机光盘检索哮喘研究文献2000余篇，其中有关成人哮喘患者血清总IgE含量测定研究的文献102篇。再从中按照该Meta分析的设计要求（纳入Meta分析的条件略）进行筛选。经筛选，得到21项符合要求的、来自国内多种期刊、分布于全国各地的研究结果，如表4-21所示。请对该资料进行Meta分析。

表4-21 21项国内成年哮喘患者血清总IgE水平（IU/ml）的研究结果

资料来源：宋士德等. 国人哮喘患者血清IgE水平的META分析. 数理医药学杂志，1998，11(3)：232-234.

4.13.2 对数据结构的分析

研究者共收集了21项符合纳入标准的研究，对每一个项目而言，其设计类型都是单因素两水平设计。

4.13.3 分析目的与统计分析方法的选择

研究者希望通过对这21项研究及其结果进行综合分析和再评价，从而得出国内成人哮喘患者血清总IgE含量的平均范围。合适的统计分析方法是多个单因素两水平设计定量资料的Meta分析。

4.13.4 SAS程序

对本资料进行单因素两水平设计一元定量资料meta分析，SAS程序名为SASTJFX4_18. SAS。

程序中分别用n_Ti、xbar_Ti和s_Ti表示各研究中实验组的样本含量、样本均值和样本标准误差；分别用n_Ci、xbar_Ci和s_Ci表示各研究中对照组的样本含量、样本均值和样本标准误差；用s_combined表示各研究的混合标准误差；用di表示各研究的效应（此处用标准均差表示）；用Ni表示各研究的样本含量；用wi表示采用固定效应模型时的权重；用widi表示wi与di的乘积；用widi2表示wi与di的平方的乘积；用number_of_study表示纳入研究的数量；用Sum_of_Ni_temp表示Ni累积的中间结果，用Sum_of_Ni表示Ni累积的最终结果；用Sum_of_wi_temp、Sum_of_wi2_temp、Sum_of_widi_temp和Sum_of_widi2_temp分别表示wi、wi2、widi和widi2累积的中间结果，用Sum_of_wi、Sum_of_wi2、Sum_of_widi和Sum_of_widi2分别表示wi、wi2、widi和widi2的累积的最终结果；用d_combined和d_bar_hat分别表示采用固定效应模型和随机效应模型时的合并效应；用w_bar、Square_of_s_w、Square_of_tao_hat和wi_bar分别表示w、s2w、τˆ2和wi；用DF表示自由度；用Q和P分别表示异质性检验的χ2值和P值；I2为异质性检验的另一个统计量；用Sum_of_wi_bar_temp和Sum_of_wi_bardi_temp分别表示wi_bar和wi_bardi累积的中间结果，用Sum_of_wi_bar和Sum_of_wi_bardi分别表示wi_bar和wi_bardi累积的最终结果；用V_hat_d_bar_hat表示Vˆ ˆ(d)；用FLOW、FUP及LOW、UP分别表示采用固定效应模型和采用随机效应模型时合并效应的95%置信区间的下限和上限。

值得一提的是，程序中多个变量是用下画线将几个单词连接在一起形成的。

4.13.5 主要分析结果及解释

异质性检验的结果为：Q=χ2=541.309，df=21−1=20，P=0<0.05，I2=0.96305>0.5，因此认为各研究同质性很差，应采用随机效应模型进行综合分析，结果表明：平均效应尺度（SMD）为3.80倍标准误差，其95% 置信区间为3.05～4.54倍标准误差。由于置信区间位于0的右侧，所以21项研究的综合结果表明，哮喘对成人血清IgE水平的影响有统计学意义。专业结论为：哮喘对成人血清IgE水平具有严重影响，成人患哮喘后血清IgE水平显著上升，平均上升水平为3.80倍标准误差，其95%置信区间为3.05～4.54倍标准误差。