半变系数伽马脆弱模型惩罚部分似然估计

为了更好地分析对数风险函数与协变量之间复杂的非线性关系，提出一种半变系数伽马脆弱模型并给出其估计方法．首先，应用B-样条将半变系数伽马脆弱模型近似转化为线性伽马脆弱模型，然后运用惩罚部分似然法估计转化后模型的线性参数，随后采用近似轮廓似然法并运用黄金搜索算法估计随机效应的参数；在通过迭代获得转化后的线性系数以及随机效应参数的估计以后，运用B-样条得到变系数函数的估计．经蒙特卡罗模拟研究发现，该方法可以给出协变量的线性参数以及变系数函数较为精准、稳定的估计，是分析协变量对于风险率影响的有效方法．最后，应用所提出的方法分析了NCCTG肺癌数据