指定删失比例的生存数据模拟及R实现

作者：蔡丽馨; 仲子航; 杨旻; 于全骥; 周佳薇; 倪森淼; 于浩^*; 柏建岭^*

来源：中国卫生统计, 2022, 39(01): 143-148.

• 生存数据
• 模拟
• 删失比例
• Cox比例风险模型

目的 探讨生成具有指定删失比例的模拟数据的新方法，并编写完整的模拟数据产生的R代码，方便使用。方法 基于Cox比例风险模型框架分别推导考虑和不考虑协变量情况下的删失参数的闭式表达式或核密度估计，并基于R 4.0.2软件编写函数模拟产生满足指定删失比例的生存数据。通过1000次模拟数据来验证生成数据的删失比例是否与指定的删失比例一致。结果 函数CenDatNoCov、CenDatBin、CenDatNorm和CenDatMixed可分别返回不考虑协变量，协变量为二项分布、正态分布以及混合分布情况下的删失数据，其中后三者内嵌的函数CensProbBin, CensProbNorm, CensProbMixed可分别得到相应情况下的删失参数值。模拟生成数据的删失比例与指定的删失比例近似相等。结论 本研究所提供的方法及编写的R函数可有效生成指定删失比例的数据。