1. 導(dǎo)入數(shù)據(jù)
webuse drugtr,clear
*webuse -- Use dataset from Stata website
stset,clear
*上一行命令導(dǎo)入進(jìn)來的數(shù)據(jù)是“生存時(shí)間數(shù)據(jù)”，會(huì)比最原始的數(shù)據(jù)多四列遭京，這一行命令的目的是重置為原始數(shù)據(jù)杰扫，之所以這樣做元践，是為了顯示數(shù)據(jù)處理的全部過程
codebook drug
codebook studytime
codebook died
codebook age
*codebook -- Describe data contents
2. 將原始數(shù)據(jù)轉(zhuǎn)換為生存數(shù)據(jù)
stset studytime,failure(died==1)
*指定生存時(shí)間為studytime列萨咳，終止(failure)事件為died列族展，其中值為1表示終止
stsum
stdescribe
sts list
*以上幾個(gè)為描述命令
3. KM生存曲線
sts graph,by(drug)
*畫KM生存曲線
sts graph if age < 50,by(drug)
sts graph,by(drug) risktable
*生存曲線下方加表格
4. log-rank檢驗(yàn)
sts test drug
*檢驗(yàn)兩組差異是否顯著窝稿，默認(rèn)為log-rank檢驗(yàn)
5. Cox比例風(fēng)險(xiǎn)模型
stcox drug
stcox drug age
*cox比例風(fēng)險(xiǎn)模型从绘，可以解讀為控制另一個(gè)（些）變量，當(dāng)前變量對(duì)發(fā)生終點(diǎn)事件風(fēng)險(xiǎn)的影響
6. PH假定的檢驗(yàn)
estat phtest
*Tests of proportional-hazards assumption撵幽，零假設(shè)為納入Cox回歸模型的變量滿足PH假定
stphplot, by(drug)adjust(age)
*by(drug)adjust(age)表示控制age變量看drug對(duì)發(fā)生終點(diǎn)事件的影響
*從圖形看PH假定灯荧，若兩條線近似平行，可以認(rèn)為是成立的
