Kaplan-MeierでNumber at riskの計算めんどくせぇとか
思っていたら、素晴らしい方がプログラムしてくれていた
リンク：http://beetama.blog14.fc2.com/blog-entry-901.html
こいつはいいね！
# LDL/HDL比での分析 L/H比が1.7で差が最大となるが、それでも有意差はでなかった
ds$Lhratio1<-ifelse(ds$Lhratio>1.7,0,1)
ds.sf<-survfit(Surv(evday, MACCE==1)~Lhratio1,data=ds)
kmplot(ds.sf,mark='',col.surv=1:4,lty.surv=1:4,lwd.surv=2,xlab="Days",ylab="event free",group.names=c("L/H>1.7","L/H<1.7"),loc.legend='bottomright')
