2020-01-06单基因生存曲线

数据格式

> head(dangenen)                          

id                          status   time         FAMB3H-AS1

 TCGA-B6-A408      0     7455         up

 TCGA-E9-A1N6      1     678         up

 TCGA-BH-A18P      1     921      down

 TCGA-A2-A0CK      0    7455       down

 TCGA-AC-A3TM      0   7455     down

 TCGA-A2-A0SW      1   1365       down

a=dangenen[,"FAMB3H"]=="down"

diff=survdiff(Surv(time, status) ~a,data = dangenen)

summary(diff)

pValue=1-pchisq(diff$chisq,df=1)

if(pValue<0.001){

  pValue=signif(pValue,4)

  pValue=format(pValue, scientific = TRUE)

}else{

  pValue=round(pValue,3)

}

fit <- survfit(Surv(time, status) ~ a, data = dangenen)

pdf(file="survival.pdf",

    width=6,

    height=6)

plot(fit,

    lwd=2,

    col=c("red","blue"),

    xlab="Time (year)",

    mark.time=T,

    ylab="Survival rate",)

legend("topright",

      c("High expression","Low expression"),

      lwd=2,

      col=c("red","blue"))

dev.off()

summary(fit)