R/cutT.R

In flexrsurv: Flexible Relative Survival Analysis

cutT<-function(T, step, mult=2L, min=as.integer(mult), max=+Inf){
	# cut T in Nstep regular sub intervals such that T/Nstep ~step
	# mult : Nstep = mult * i
	mult <- as.integer(mult)
	min <- as.integer(min)
	Nstep<- as.integer(T %/% step)
	Nstep<- as.integer(((Nstep %/% mult) + 1L) * mult)
	Nstep<-ifelse(Nstep<min, min, Nstep)
	Nstep<-ifelse(Nstep>max, max, Nstep)
	data.frame(NstepT=Nstep, stepT=T/Nstep)
}


cutTfromto<-function(fromT, toT, step, mult=2L, min=as.integer(mult), max=+Inf){
	# cut T in Nstep regular sub intervals such that T/Nstep ~step
	# mult : Nstep est une multiple de mult
	mult <- as.integer(mult)
	min <- as.integer(min)
	
	Nstep<- as.integer((toT - fromT)  %/% step)
	Nstep<- as.integer(((Nstep %/% mult) + 1L) * mult)
	Nstep<-ifelse(Nstep<min, min, Nstep)
	Nstep<-ifelse(Nstep>max, max, Nstep)
	data.frame(NstepT=Nstep, stepT=(toT - fromT) /Nstep)
}