result = npdf_cox( Surv(time, status) ~ x, groups = family, data = weibdata2030,
K = 2, eps_conv=10^-4)
plot( result, K = 2 )
plot( result, type = "na" )
plot( result, type = "na", cols=ifelse(result$belonging==2, "purple", "black"),
na_opts = list(xlim=c(0,200), ylim=c(0,2),
xlab="Follow-up days",
ylab="Nelson-Aalen estimate"))
Run the code above in your browser using DataLab