KMのPlotのマークを変える

km <- NULL
km.summary.table <- NULL
km <- survfit(Surv((intervention.days/30.4375),intervention.event==1)~1, data=FTUL1107, na.action
= na.omit, conf.type="log-log")
windows(width=7, height=7); par(lwd=1, las=1, family="sans", cex=1, mgp=c(2.5,1,0))
mar <- par("mar")
mar[1] <- mar[1] + 1 + 0.5
par(mar=mar)
opar <- par(mar = mar)
on.exit(par(opar))

#mark .time=TRUEでひげつける。pch=""か数字でマーク変更

#http ://cse.naro.affrc.go.jp/takezawa/r-tips/r/53.html

# 

plot(km, bty="l", col=1:32, lty=1, lwd=1, conf.int=FALSE, mark.time=TRUE,pch ="|", xlim=c(0,60),
xlab="intervention.days", ylab="Probability")
xticks <- axTicks(1)
n.atrisk <- nrisk(km, xticks)
axis(1, at = xticks, labels = n.atrisk, line = 3, tick = FALSE)
title(xlab = "Number at risk", line = 3, adj = 0)
summary(km)
summary.km(survfit=km)