Description Usage Arguments Details Value Author(s) References Examples

The method provides one step ahead imputed values for tied censored observations.

1 2 | ```
imputeYn.extra(Y, delta, hc.Yn, method = "km.TPQ", trans.sprob=NULL,
stime2=NULL, sprob2=NULL, trace=F)
``` |

`Y` |
response. Typically the logarithmic of the survival time. |

`delta` |
status; it includes value 1 for uncensored and value 0 for censored subject. |

`hc.Yn` |
set of the lifetimes for the last largest censored observations. |

`method` |
one of "it.PDQ" (repeated predicted difference quantity or called iterative (Khan and Shaw, 2013a)), "km.TPQ" (Kaplan-Meier trend predicted quantity or called "extrapolation" (Khan and Shaw, 2012a)). Default is "km.TPQ". |

`trans.sprob` |
use only for "km.TPQ". Transformation for the survival probabilities. This transformation is needed if survival probability versus survival time plot is not linear. It takes "exp" for exponential transformation, and any value for respective power transformation. Default is NULL. |

`stime2` |
use only for "km.TPQ". Survival times after necessary transformation if needed in order to obtain a linear relationship between the survival probability and survival time. Default is NULL. |

`sprob2` |
use only for "km.TPQ". Survival probability after necessary transformation if needed in order to obtain a linear relationship between the survival probability and survival time. Default is NULL. |

`trace` |
If TRUE then Kaplan-Meier survival plots will be printed for both data- the original and the data with imputed values. Default is FALSE. |

The method is developed for imputing the largest censored observations (can be seen often for heavy censored data) if Kaplan-Meier weights are involved in modeling. For example, if weighted least squares is used then the extra imputing methods will provide one step ahead of the Efron's (1967) tail correction approach. These methods satisfy the Efron's tail correction, achieve less biased estimates, and impute the largest censored observations that are also tied. Furthermore the two methods satisfy the right censoring assumption. If there is heavy censoring toward the right for right censored data, that is found as typical case for many areas like indistry, clinical trials etc, then this function provides the sensible imputations for those censored observations that are the largest and also tied observations. Details are discussed in Khan and Shaw (2013a).

It includes sorted lifetimes, censoring indicators, sorted lifetimes after imputation, censoring indicators after imputation, censored lifetimes for the Y(n)+ observations, imputed lifetimes for the Y(n)+ observations, the survival times, the survival probabilities, the survival times after transformation, the survival probabilities after transformation, the transformation used, and trace.

Hasinur Rahaman Khan and Ewart Shaw

Efron, B. (1967). The two sample problem with censored data. In Proceedings of the fifth Berkeley symposium on mathematical statistics and probability, Vol. 4, p. 831-853.

Khan and Shaw. (2013a). On Dealing with Censored Largest Observations under Weighted Least Squares. CRiSM working paper, Department of Statistics, University of Warwick, UK, No. 13-07. Also available in http://arxiv.org/abs/1312.2533.

Khan and Shaw (2013b). Variable Selection with The Modified Buckley-James Method and The Dantzig Selector for High-dimensional Survival Data. Proceedings 59th ISI World Statistics Congress, 25-30 August 2013, Hong Kong, p. 4239-4244.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
## For Channing House data (heavy censored data)##
## Not run: require(package="boot")
## Not run: time.ch<-channing[1:97,]$time #for male
## Not run: delta.ch<-channing[1:97,]$cens # for male
## Not run: hc.Yn.m<-rep(137,19) # there are 19 last largest censored male each has 137 lifetime
## Not run: imp.PDQ<-imputeYn.extra(time.ch, delta.ch, hc.Yn=hc.Yn.m,
method="it.PDQ", trace=T)
## End(Not run)
## Not run: imp.PDQ
## Not run: imp.TPQ<-imputeYn.extra(time.ch, delta.ch, hc.Yn=hc.Yn.m,
method = "km.TPQ", trace=T)
## End(Not run)
## Not run: imp.TPQ
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.