含特定比例均匀随机删失生存数据SAS模拟实现.PDFVIP

·１０５８· 中国卫生统计２０１６年１２月第３３卷第６期 含特定比例均匀随机删失生存数据的ＳＡＳ模拟实现 １ １ ２△ 肖媛媛 　许传志　赵耐青 　　删失（ｃｅｎｓｏｒｉｎｇ）是生存分析时经常会遇到的一个 １－Ｆ（Ｙ）也服从在［０，１］上的均匀分布，即Ｕ＝ 得知， 实际问题，为生存数据的处理分析带来了一定困难。 ｅｘｐ［－Ｈ（ｔ）ｅｘｐ（′Ｘ）］服从［０，１］上的均匀分布 ０ β 针对这一问题，在实际研究中往往需要通过随机模拟 Ｕ［０，１］。 来产生含一定删失比例的特定生存分布数据，供探索 假设生存时间为Ｔ，对前式取反函数，得到： 性研究或对比性研究使用。文献检索后发现，目前国 －１ Ｔ Ｔ＝Ｈ （ｔ）［－ｌｏｇ（Ｕ）ｅｘｐ（－ Ｘ）］ ０ β 内外针对删失生存数据的模拟研究数量较少［１－３］。 基于上述思想，根据目标生存分布所对应的累积 按照其发生是否与研究对象的其他特征有关，删 风险函数表达式Ｈ 和随机模拟产生的均匀分布Ｕ，即 ０ 失分为随机删失（ｒａｎｄｏｍｃｅｎｓｏｒｉｎｇ）和非随机删失 可产生目标生存时间。Ｗｅｉｂｕｌｌ分布的累积风险函数表 （ｎｏｎｒａｎｄｏｍｃｅｎｓｏｒｉｎｇ）两类，即非信息删失（ｎｏｎｉｎ 达式为： ｆｏｒｍａｔｉｖｅｃｅｎｓｏｒｉｎｇ）和信息删失（ｉｎｆｏｒｍａｔｉｖｅｃｅｎｓｏ γ－１ Ｈ（ｔ）＝ ｔ ０ λγ ［４］ ｒｉｎｇ） 。由于非随机删失的影响因素众多，故其特征 因此可得到其生存时间Ｔ的计算公式为： 难以全面捕捉，因此目前在模拟产生含删失的生存分 ｌｏｇ（Ｕ） １／γ Ｔ＝( － ) ｅｘｐ（′Ｘ） 析数据时，大都假定为完全随机删失。 λ β 本文针对现有模拟思路和方法中存在的某些问题 在进一步产生含随机删失的生存数据时，一般的 或疑议，提出一种模拟产生含特定比例均匀随机删失 思路为：假设生存时间Ｔ和删失时间Ｃ服从相互独立 生存数据的思路，以及实现这一思路的ＳＡＳ宏程序， 的两个分布Ｓ（Ｔ）和Ｂ（Ｃ）。在模拟研究中，Ｓ（Ｔ）的参 为科学工作者开展生存分析提供便利。 数已经设定，可以顺利产生生存时间Ｔ。再根据定义 ｉ 的删失分布类型（即删失发生时间的概率分布），如指 现有方法回顾 数分布或均匀分布，结合设定的删失比例，采用