https://github.com/cran/emplik
Tip revision: 19969178b30432d77fffc9a178392fb0b06c7718 authored by Mai Zhou on 03 July 2022, 02:20:02 UTC
version 1.2
version 1.2
Tip revision: 1996917
kmc.clean.R
kmc.clean <- function(Xtime, delta){
##TASK: 1. sort Time
## 2. delete the smallest censored obs. so that the first obs. is uncen! the last also uncen.
n <- length(Xtime)
if( all(delta == 0) ) stop('All obs. are censored?')
dataOrder <- order(Xtime, -delta)
kmc.time <- Xtime[dataOrder]
delta <- delta[dataOrder] ### changed 10/2018
#### tmp <- sort(kmc.time,index.return=TRUE)
#### kmc.time <- kmc.time[tmp$ix]
#### delta <- delta[tmp$ix]
delta[n] <- 1
FirstUnCenLocation <- which(delta==1)[1]
if (FirstUnCenLocation==n) stop('Only one uncensored point.')
if (FirstUnCenLocation!=1){
delta <- delta[FirstUnCenLocation:n]
kmc.time <- kmc.time[FirstUnCenLocation:n]
}
list(kmc.time=kmc.time, delta=delta)
}