I have specified error = "greenwood" in the survfit function: summary(survfit(Surv(time, status) ~ 1,error = "greenwood", data=aml)) The output does not appear to give variance, but standard error. Preferably I would like variance to be output, but is it okay to convert std errror to variance, and how can I do this correctly? Thanks Murray