1、例如,为了研究两种降压药物对血压的控制效果是否存在差异,研究者会对两个人群服药后在不同时间点记录血压值,然后评价降压效果。或者对两组动物分别施加两种干预,连续记录多个时间点的结局,然后比较两种干预的效果。这种设计可以用如下示意图表示:另外,有时研究只需要收集一个时间点的数据,但是一个研究对象会提供多个部位的数据点。例如,研究者想评价冠心病患者在冠脉搭桥术后应用阿司匹林是否可以有效降低患者血管的再堵塞,评价的方法是术后1年做冠脉造影观察血管是否堵塞,但是每个患者可能会在同一次手术中对多条冠状动脉血管进行搭桥,因此有的患者可能会贡献多组数据。以上两种设计,不管是临床试验还是动物试验都非常常见,它的
2、特点在于数据间非独立,同一个体间数据具有相关性。对于这样的设计类型,该如何分析呢?今天我们来介绍另外一种非常好的方法广义估计方程()。既可以处理连续型结局变量也可以处理分类型结局变量,它实际上代表了一种模型类别,即在传统模型的基础上对相关性数据进行了校正,可以拟合回归、泊松回归、回归、一般线性回归等广义线性模型。本文将以阿司匹林预防冠脉搭桥后血管再堵塞为例介绍运用进行的操作方法。以下为数据格式:表1. 数据格式每名患者贡献数据量不等。如编号为1的患者只对一根血管进行了搭桥手术,编号为2的患者则有两根血管进行搭桥手术。表2. 变量赋值(注:本例中数据纯属虚构,分析结果不能产生任何结论。性别为待调
3、整变量。)二、分析方法1. 数据录入首先在变量视图( )中新建上述表2中变量,然后在数据视图( )中录入数据。2. 选择 3. 选项设置 在模块中,将患者选入 框,搭桥血管选入 框。在 的下拉列表中选择同一患者贡献数据的相关性结构。其有5个选项,分别为:各数据间独立,即同一患者贡献的数据可以视为来源于不同个体。(1):自相关,用于不同时间点的数据,相邻时间点相关性最大,时间间距越大相关性越小。等相关,各数据点的相关性均相等。相邻相关,仅临近的1个数据相关。不限定相关结构,由数据本身决定。本例中,同一患者搭桥的血管互相之间应等同看待,两两间相关性应相等,因此选择等相关更合适。其他维持默认。 在
4、模块中设置模型类型在此界面有6种模型类型可选: :结局变量为连续型变量结局变量为有序分类变量结局变量为计数结局变量为二分类变量结局变量为分布变量自定义以上6大类类型下还需选择连接函数。本例中结局变量为二分类变量,回归应用更加普遍,因此我们选择 来拟合校正数据相关性的回归模型。 在模块指定结局变量将变量选入 ,在 中选择最小值作为参考值(本例中0为最小值,代表未发生结局)。 在模块中将要纳入模型的和变量选入中。应放入连续型变量,中放入分类变量,但是本例中两自变量均为二分类变量,放入处也可。 在模块中将和选入框内,还可纳入交互项,本例中不纳入。 在模块中除默认外,另外勾选 以获得值,然后点击。三、
5、结果解读结果中首先给出的是一些模型基本信息,在此略过。这里仅介绍模型参数估计值结果:从结果可以看到,变量的系数值(B)为-1.077,(B)为0.341,即值,0.046,说明术后使用阿司匹林的患者1年后发生血管再堵塞的风险是使用安慰剂组的0.341倍且有统计学意义,即阿司匹林具有保护作用。另外我们纳入校正的性别变量 值为0.692,但0.531未达显著水平。如果我们不考虑数据间的相关性,直接采用传统的回归模型,那么将得到以下结果:从结果中可以看到,尽管的值和模型中相差不大,但其P值却不再显著。从本例中可以看到,如果在数据分析中方法选用错误将可能得不到正确结论。四、结果汇总冠脉搭桥术后患者应用
6、阿司匹林可以有效降低发生血管堵塞的风险。服用阿司匹林组患者发生血管堵塞风险是服用安慰剂组患者的0.341倍(0.046)。五、总结和拓展对于文首提及的两个例子,常见的分析方法有重复测量方差分析和多水平模型,但是重复测量方差分析要求结局变量为连续变量,不适用于分类变量。多水平模型处理相关性数据时非常灵活,结局可以为连续变量,也可以为分类变量,应用非常广泛。的应用似乎没有前两种广泛,但其具有非常好的特性。既可以处理连续型结局变量也可以处理分类型结局变量,其优势在于,即使设定的数据间相关结构和实际有偏差,在样本量较大时其模型估计参数仍然具有无偏性。模型中自变量系数估计值准确性的论证相较于多水平模型更加充分,因此部分研究者更加推荐使用。虽然模型中数据间相关结构指定错误时模型系数也具有无偏性,但前提是样本量要足够大。在同样的样本量下,正确的相关结构更可能得到准确的参数估计和置信区间,因此应仔细分析数据类型,选择最可能正确的相关结构。关于的内容较为复杂,在此无法详述,对有兴趣的读者或需要使用模型的读者,推荐阅读专业书籍和文献。
copyright@ 2008-2023 冰点文库 网站版权所有
经营许可证编号:鄂ICP备19020893号-2