Kaplan-Meier and Nelson-Aalen estimates with and without survfit function from survival package

Data reading:


rm(list=ls(all=TRUE)) # wipe out everything else
require(survival)
data = aml[aml$x == "Maintained",]
data

Y-axis on both sides in Kaplan-Meier Survival Curve

Here is an example of how to draw Y-axis on both sides in Kaplan-Meier Survival Curve, with same label, axis-titles and positions. This should work for the general plots as well.