Survival.jl icon indicating copy to clipboard operation
Survival.jl copied to clipboard

unclear std error for Kaplan Meier estimator

Open tbeason opened this issue 3 years ago • 1 comments

The standard error for the Kaplan-Meier estimator returns what I would consider to be only "half" of the correct value. Perhaps this is just a misunderstanding on my part, but I think it should return the square root of the variance of the survival function. Right now, it does not. To actually make it the standard error of the survival curve, you need to multiply the two (AKA km.survival .* km.stderr).

Should I make a PR to fix this, or was it intentional?

tbeason avatar Jul 26 '20 22:07 tbeason

Currently it computes the pointwise standard error of the log of the survivor function using Greenwood's formula. I forget why exactly I opted for that but likely it was because that's also what R does (which is what this implementation tests against).

ararslan avatar Jul 10 '22 20:07 ararslan