有偏抽样下带终止时间和带信息观察时间的纵向数据分析.pdfVIP

第36卷第6期 应用数学学报 V01．36No．6 2013年11月ACTAMATHEMATICAEAPPLICATAESINICANov．，2013 有偏抽样下带终止时间和带信息 观察时间的纵向数据分析木 苗瑞 (中国科学院数学与系统科学研究院，北京100190) 李怿 (中国科学技术大学统计与金融系，合肥230026) (E-mail：wuguily@mail．ustc．edu．an) 摘要纵向数据经常出现在生物医学等实际领域，而且在许多情况下，观察时间和终止时间 是带信息的，同时可能会存在有偏抽样．本文在有偏抽样下，对于带信息观察时间和终止时间 的纵向数据，提出了—个联合建模方法，并利用借代力量方法给出了模型参数的估计，同时获 得了这些估计的渐近性质．最后，提出一种模型检验方法． 关键词有偏抽样；带信息观察；终止时间；联合建模；潜在变量；估计方程 MR(2000)主题分类62G05；62N01 中图分类0212．7 1 引言 纵向数据是对某些感兴趣的个体进行多次观察或重复测量所获得的数据，它的一 个重要特点是每个个体的各个观测值之间存在相关性．这类数据几乎出现在所有应用 科学领域，包括生命科学、医学、生态学、工程技术、环境科学和经济科学等．例如，多 次观测的癌症病人的肿瘤数；某种疾病患者的多次医疗费用；HIV感染过程中的CD4 数目；气候变化对各国天然气的需求量；年龄、性别、人种、教育和工龄对工资和福利 的影响程度；在临床试验中，为了比较新的治疗方法，对病人治疗后进行基础测量及其 后的跟踪测量，这些测量的数据都是纵向数据．文献中已经有了许多分析这类数据的方 中常用的方法进行了一个综述．但是以上这些建模方法，都有一个很强的假设：在给定 本文2013年8月21日收到．2013年9月11日收到修改稿． +国家自然科学基金11171330资助项目 962 应 用 数 学 学 报 36卷 协变量条件下，纵向响应变量与观察时间之间是相互独立的． 近些年，基于带信息观察时间的纵向数据的研究引起了广泛的关注[7-10】．例如， Lin等【7】基于一类边际回归模型提出了一个逆强度加权估计．Sun等【9】利用一个潜 在变量建立了纵向响应变量和观察时间的一般联合模型，并给出了模型参数的统计推 断．Liang等110】利用两个潜在变量对纵向响应变量和观察时间进行联合建模分析．所 有这些方法都是用来分析在没有终止时间情形下的带信息观察的纵向响应变量． 在实际中，经常存在一个终止时间，比如死亡．而且终止时间往往与纵向过程和观 察过程具有相依性．例如，病情较严重的病人经常在较短的时间死亡，而且治疗的费用 与死亡时间是密切相关的【11】．Liu等[12】说明了如果不考虑相依的终止时间可能会导 致对治疗时间所建立的模型产生有偏估计．同时，Liu等[13】也证明了如果忽略相依的 终止时间会导致对纵向医疗费用所建立的模型产生有偏估计．然而，基于带信息观察时 间和相依终止时间的纵向数据的研究非常有限．最近，Liu等[11】针对带信息观察时间 和相依终止时间的纵向数据提出了一个联合随机效应模型，在该模型中，随机效应的分 布是确定的．He等【14】提出了一类带有公共脆弱变量的模型来分析与观察时间和终止 时间相关的面板计数数据．他们在模型中需要假定一个脆弱变量服从正态分布，但实际 中，我们很难从一个特定的数据集中确定合适的脆弱变量的分布，而且错误的脆弱变量 分布的假定会导致不正确的推断结果．Sun等【15】借助两个潜在变量，对纵向响应变 量、观察时间以及终止时间建立了一般的联合模型． 另外，当不规则的抽样时间受到协变量以外的其他因素影响时，往往需要对均值反 应模型引入一个有偏抽样设计来达到更好的推断效果[7,16,17】．因此，我们迫切需要找到 一种模型能直接刻画纵向数据，而且能够同时考虑到有偏抽样，带信息观察时间和相依 终止时间． 在本文中