加法风险模型的半参数估计方法及其应用.pdfVIP

第35卷总第94朋 两 北 民 族 大 学 学 报 (自然科学版) Vol_35．No2 2014年6月 JournalofNorthwestUniversity for Nati0nalities(NaturalScience) June，2O14 加法风险模型的半参数估计方法及其应用 金凌辉 ，郭丽莎 ，。 (1．武汉科技大学 城市学院公共课部，湖北 武汉 430083：2．叶l南民族大学 数学与统计学院，湖北 武汉 430074：3 武汉大学 数学与统计学院，湖北 武汉 430072) I摘 要I 加法风险模型足生存分析 的一种重要回归模型，但传统的偏似然方法不能直接用于该模型的估计， 凶此文章介绍了加法模型的一种半参数估计方法及其渐近性质，并讨论了该方法在病例一队列设计下的应用． I关键词I 生存分析；加法风险模型；伪得分；渐进性质；病例一队列设计 [中图分类号】 0212 [文献标识码】 A [文章编号】1009-2102(2014)02-0001--05 Cox¨ 1972年提出的比例危险率模型是生存分析中最重要的回归模型．但在相关研究中，人们往 往还会关注所暴露的风险差异 (riskdifference)，因此 Aalen提出了适用于此 目的的加法风险模 型．Cox和0akes 、Buckley 、Lin和Ying 等许多学者都对加法模型进行了研究，加法风险模型也 得到了广泛的应用．然而，在流行病学研究中，往往需要对数量很大的研究对象进行几年或更长时间 的随访研究，但研究时问截止时却只有少量的随访对象产生了感兴趣的病例，而其他大多数对象所能 观测到的仅为删失数据．显然，若收集所有研究对象的协变量信息用于统计分析，将付出极高的成本． 因此，我们可将加法模型用于Prentice提出的病例一队列设计 (case—cohortdesign)进行成本控制． 1 记号和模型 令 t为观测时间， 为观测截止时问，71为失效时问， 为删失时问，％=min(L ， =I(T C)， (为P维随时问改变的协变量向量，Ao(f)为基准危险率， 为P维回归参数向量，若在给定协变 量 的条件下，危险率函数 (f)满足关系 (f；z)= n(f)P “‘， (1) 则模型称为 Cox模型或 比例危险率模型，运用所谓的偏似然方法 (Partiallikelihoodapproach) 可以得到未知参数的 估计，现简介如下： 首先记X(1) (2)…X()为Xf的次序统计量， f)和z(f)(t)分别是对应于 f【)的示性函数 与协变量，记R(t)为 t时刻的风险集，则具有协变量Z㈤(t)的个体在时刻 ㈩死亡的概率为 z(x(f)IZ，(f))／∑_，∈(，)((f)ZI(f))．由式(1)可知，将分子及分母中所含的 (f)约掉后，此概率 实际为P卢j)【‘／∑，E(．) ，‘．再将所有时问点的概率都乘起来后，便得到偏似然函数 = 骢 ， [收稿 日期] 2014-05-·20 [基金项 目] 湖北省教育ft．：2013科研计划项 目(B2013260)，2012中央高校专项资金资助青年项 目(CZQ12015) [通讯作者] ● [作者简介】 金凌辉 (1979——)，男 (瑶族)，广西人，讲师，博士，主要从事生物统计学方 的研究． 一 1一 两边取对数后并求导，便可得到关于 的得分函数 ，： ∽ 一 =，1 一 ㈦ 再由 ()=0可解得 的极大偏似然估计． 若危险率函数Aft)满足关系 ；Z)= 。)+Pz(t)， (4) 则模型称为力I【法风险模型．由于