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

来源：大连理工大学学报杂志2018年第6期北京时间：

作者：张中文,王晓光,[宋立新]

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

关键词：伽马脆弱模型B-样条变系数模型惩罚部分似然估计黄金搜索算法

基金资助：国家自然科学基金资助项目(114710651137107781502891)；全国统计科学研究项目(2018LY57)；山东省统计科研重点课题资助项目(KT16244)；山东省医药卫生科技发展计划资助项目(2016WS0009)；山东省软科学研究计划项目(2018RKB14103).