生存分析及其R程序学习教案

1、会计学1生存分析生存分析(fnx)及其及其R程序程序第一页，共42页。第1页/共42页第二页，共42页。第2页/共42页第三页，共42页。第3页/共42页第四页，共42页。第4页/共42页第五页，共42页。第5页/共42页第六页，共42页。第6页/共42页第七页，共42页。第7页/共42页第八页，共42页。第8页/共42页第九页，共42页。第9页/共42页第十页，共42页。ttnttntnttTttTtPthtt )()()(lim)(lim)(00第10页/共42页第十一页，共42页。nh1(t)是一种上升的曲线，表示危险率随时间变化而增加，如急性白血病患者治疗(zhlio)无效其危险率随时

2、间呈增加趋势；nh2(t)的曲线为下降趋势，表示危险率逐渐减小，如意外事故造成的外伤经有效治疗(zhlio)后死亡的危险性逐渐减小；nh3(t)为一种稳定的危险率函数，如某些慢性病患者在稳定期，其危险率基本不变。n 第11页/共42页第十二页，共42页。第12页/共42页第十三页，共42页。tdtthtS0)(exp)(第13页/共42页第十四页，共42页。dttdStSdttSdntSndttdntnttnttntnttTttTtPthtt)()(1)()0()()0(1)()(1)()()(limlim)(00)(ln)()(1)()(1)(tStdStSdtdttdStSdtth)(ex

3、p1)(exp)(dtthdtthtS第14页/共42页第十五页，共42页。第15页/共42页第十六页，共42页。如果(rgu)风险比与时间无关，及任何时刻，两组的风险比值是相等的，则称为比例风险（PH），否则称为非比例风险。第16页/共42页第十七页，共42页。第17页/共42页第十八页，共42页。第18页/共42页第十九页，共42页。第19页/共42页第二十页，共42页。第20页/共42页第二十一页，共42页。序号序号确诊后确诊后月数月数期内截期内截尾人数尾人数期内死亡期内死亡人数人数期初观期初观察人数察人数校正校正人数人数死亡死亡概率概率生存生存概率概率生存率生存率标准误标准误第21页/

4、共42页第二十二页，共42页。第22页/共42页第二十三页，共42页。第23页/共42页第二十四页，共42页。第24页/共42页第二十五页，共42页。第25页/共42页第二十六页，共42页。序号序号生存生存时间时间（t)期初期初人数人数（n)死亡数死亡数（d） 甲组甲组 乙组乙组期初人数（n1)死亡数（d1)理论死亡数（E1）期初人数（n2)死亡数（d2)理论死亡数（E2）=d1-E1=w(d1-E1) (w=n)第26页/共42页第二十七页，共42页。变量名变量名描述描述CASE编号INST机构代码TX实验代码（1：标准，2：处理）GRADE和正常细胞的区别度COND身体状况SITE病变位置

5、T.STAGE癌症T分期N.STAGE癌症N分期ENTRY.DT进入实验日期STATUS删失（0：右删失，1：死亡）TIME如未删失则是存活时间否则是最后记录时间第27页/共42页第二十八页，共42页。aa-table(COND) aamydata-my,-c(1,11)library(survival)sd-survfit(Surv(TIME,as.numeric(STATUS)1,data=mydata)plot(sd,lty=1,col=“red”,ylab=“s(t)”,xlab=“t”,main=“survival”) #生存曲线fit0，决定了生存时间分布的形状(xngzhun)，m越小，分布越偏，当m=3.57时，分布近似为正态分布。m1为老化现