文档介绍:()library(survival)library(help=survival)data(aml)ami#查看你的系统中所有己装的libraries#加载一个library,也口J使用鼠标来做#seethelistofavailablefunctionsanddatasets.#加载数据ami#seethedata#《生存分析》需要的两个librarieslibrary(survival)#包含我们教材中所有的数据library(KMsurv) #包"#若想使数据中的变量能够在R中作为一个变量使用,useattach():attach(aids)infectdetach(aids)—.SurvivalObjects#由函数Surv产生一个survivalobjectSurv(time,time2,event,type)time:survivaltimetime2:endingtimeifitisintervalcensoredevent:,O=rightcensored,l=eventattime,2=leftcensored,3=:indicatingwhetheritisrightcensored,#假设一个观测值为[2,3],则R语句为Surv(time=2,time2=3,event=3,type="interval")#右删失数据attach(aml)Surv(time,status)detach(aml)#左截断右删失数据data(psych)attach(psych)<-Surv(age,age+time,death)[1](51,52](58,59](55,57](28,50](21,51+](19,47](25,57][22](29,63+](35,65+](32,67](36,76](32,71+]detach(psych)三TheKaplan-MeierEstimates#计算生存曲线:survfit,有三个重要的参数:formula,,(formula,=,="log")"plai #〃log〃 #对H(t)"logJog" #log-transformedinourbook#在新版的R里,这个"是必要的。#在新版的R里,这个T是必要的。myfit=survfit(Surv(time,status)~l,data=kidney)attach(kidney)myfit=survfit(Surv(time,status)~1)a=summary(survfit(Surv(time,status