R语言-survival包

R语言实战 · 公众号 · · 2017-07-17 17:39

正文

假设我们想调查一些因素。以生存为例，我们可以比较糖尿病患者和使用二甲双胍的患者的经验注射胰岛素作为主要治疗方式。有一些证据证明二甲双胍有一个积极的特别是在癌症中，但确是由事实引起的这是一种rst线疗法:二甲双胍的患者平均年龄更小。患糖尿病的时间比使用胰岛素的人要短。年轻人“活得更久”并不是一个特别新奇的观察。

fdata =7,]
       fdata$age2        fdata$group 7) + 1*(fdata$flc.grp >9), levels=1:3, labels=c("FLC < 3.38", "3.38 - 4.71", "FLC > 4.71"))
       sfit1        plot(sfit1, mark.time=F, col=c(1,2,4), lty=1, lwd=2, xscale=365.25, xlab="Years from Sample", ylab="Survival")
      text(c(11.1, 10.5, 7.5)*365.25, c(.88, .57, .4),c("FLC < 3.38", "3.38 - 4.71", "FLC > 4.71"), col=c(1,2,4))

       temp        dd    # Select subjects
       set.seed(1978)
       select        for (j in 1:dd[2]) {
   for (k in 1:dd[3]) {
      n    for (i in 1:2) {
   indx      select[i,j,k] temp[i,j,k])))
      }
indx select[3,j,k]    }
    }
     data2      sfit2      plot(sfit2,col=c(1,2,4), lty=1, lwd=2,xscale=365.25, xlab="Years from Sample", ylab="Survival")
    lines(sfit1, col=c(1,2,4), lty=2, lwd=1,xscale=365.25)
    legend(730, .4, levels(fdata$group), lty=1, col=c(1,2,4),bty='n', lwd=2)

       tab2        tab3        rwt        fdata$rwt        sfit3        temp        temp        round(temp, 1) #show female data
       plot(sfit3, mark.time=F, col=c(1,2,4), lty=1, lwd=2,xscale=365.25, xlab="Years from Sample", ylab="Survival")
       lines(sfit3a, mark.time=F, col=c(1,2,4), lty=1, lwd=1,xscale=365.25)
       lines(sfit1, mark.time=F, col=c(1,2,4), lty=2, lwd=1, xscale=365.25)

       legend(730, .4, levels(fdata$group), lty=1, col=c(1,2,4),bty='n', lwd=2)        id         cfit         summary(cfit)$robscore
        if (exists("svykm")) { #true if the survey package is loaded
sdes dfit       }