Dear All

Does any one know if the following survival analysis exists, and if so, is it available in R?

I am trying to parametrically model survival where some events happen at time = 0. I am particularly interested in generating the 95% confidence interval around the survival function over its full range (including time=0). For example, the data below are for time (Time) to component failure (Status =1) for a set of 10 new components. Failure is recorded in 9 components, the first 3 of which failed to start (Status=1, Time=0).

> Time
 [1]  0  0  0  8 16 17 17 22 23 28
> Status
 [1] 1 1 1 1 1 1 1 1 1 0

I know how to do this in two parts (using a parametric survival model where Time>0 to develop confidence intervals over the range Time>0, and using the proportion of failures to start to generate a binomial confidence interval for Time=0). But I am hoping there is a single modelling approach - in this case it would suggest likelihood of survival at Time=0 is <1, rather than the models I am familiar with that indicate survival probabilityat Time=0 is always 1.

Any suggestion would be appreciated

Dr Terry Beutel
Senior Scientist
Agri-Science Queensland


********************************DISCLAIMER**************...{{dropped:15}}

