生存分析及其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)后死亡的危險性逐漸減?。籲h3(t)為一種穩(wěn)定的危險率函數，如某些慢性病患者在穩(wěn)定期，其危險率基本不變。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和正常細胞的區(qū)別度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時，分布近似為正態(tài)分布。m1為老化現