>SAS Code：
-
組間比較或者分層組間比較常用下面兩個語句：
proc lifetest data=Exposed plots=(survival(atrisk) logsurv);
time Days * Status(0);
strata Treatment; /* 按strata分層描述生存數(shù)據(jù)并進行組間比較 */
run;proc lifetest data=Exposed plots=(survival(atrisk) logsurv);
time Days * Status(0);
strata Sex / group=Treatment; /* 控制分層因素Sex后，比較組間差異 */
run;
一些回顧+注意點
*生存率的估計：K-M estimate, Fleming-Harrington estimate，Brewlow estimate
*生存率的標(biāo)準(zhǔn)誤：標(biāo)準(zhǔn)誤通過Greenwood公式進行計算
*生存率的置信區(qū)間：通過不同的轉(zhuǎn)換函數(shù)(g-transformation)構(gòu)造置信區(qū)間。通過CONFTYPE= 指定g-trans。
*生存時間的置信區(qū)間：通過Brookmeyer and Crowley法構(gòu)建
*組間比較：Wilcoxon test=Breslow test。LR (Likelihood ratio test) 假定指數(shù)分布。 Breslow test以期初人數(shù)為權(quán)重，更能檢測出早期差異。Log-Rank無權(quán)重，相對而已更能檢測出后期的差異。
圖示：Proc lifetest data=XXX plots=(LLS LS );
LLS：plots the log of negative log of estimated survivor functions versus the log of time.用以判斷Weibull 分布(呈直線)
LS：plots the negative log of estimated survivor functions versus time.用以判斷指數(shù)分布(呈直線)
>*Strata與Test
STRATA：Specify a variable or a set of variable that define the strata for the analysis. 可指定Test=(Logrank Wlicoxon...)。
TEST：Specify a list of numeric covariates to be tested for their association with the response survival time. 不常用。
proc lifetest data=Exposed plots=(survival(atrisk) logsurv);
time Days * Status(0);
Strata Cell; /* 檢驗層間一致性 /
Test Age prior ... ...; / 加入?yún)f(xié)變量，分析協(xié)變量與生存時間的相關(guān)性 */
run;
