data source: “ArMorr DBP DNA 12 23 2015.xlsx”

dbp1 <- survfit(Surv(fu, Death)~DBP1, data=dat)
ggkmTable(dbp1, ystratalabs=c("G", "GT", "T"), timeby=50, main="DBP 1")

dbp2 <- survfit(Surv(fu, Death)~DBP2, data=dat)
ggkmTable(dbp2, ystratalabs=c("A", "GA", "G"), timeby=50, main="DBP 2")

dbp3 <- survfit(Surv(fu, Death)~DBP3, data=dat)
ggkmTable(dbp3, ystratalabs=c("A", "AG", "G"), timeby=50, main="DBP 3")

dbp4 <- survfit(Surv(fu, Death)~DBP4, data=dat)
ggkmTable(dbp4, ystratalabs=c("C", "CT", "T"), timeby=50, main="DBP 4")

dbp5 <- survfit(Surv(fu, Death)~DBP5, data=dat)
ggkmTable(dbp5, ystratalabs=c("A", "AG", "G"), timeby=50, main="DBP 5")

dbp6 <- survfit(Surv(fu, Death)~DBP6, data=dat)
ggkmTable(dbp6, ystratalabs=c("C", "CT", "T"), timeby=50, main="DBP 6")

dbp7 <- survfit(Surv(fu, Death)~DBP7, data=dat)
ggkmTable(dbp7, ystratalabs=c("A", "AG", "G"), timeby=50, main="DBP 7")

dbp8 <- survfit(Surv(fu, Death)~DBP8, data=dat)
ggkmTable(dbp8, ystratalabs=c("C", "TC", "T"), timeby=50, main="DBP 8")

dbp9 <- survfit(Surv(fu, Death)~DBP9, data=dat)
ggkmTable(dbp9, ystratalabs=c("C", "TC", "T"), timeby=50, main="DBP 9")