Compute the Reduced Sample estimator of a survival time distribution function, from histogram data
reduced.sample(nco, cen, ncc, show=FALSE, uppercen=0)
vector of counts giving the histogram of uncensored observations (those survival times that are less than or equal to the censoring time)
vector of counts giving the histogram of censoring times
vector of counts giving the histogram of censoring times for the uncensored observations only
number of censoring times greater than the rightmost histogram breakpoint (if there are any)
Logical value controlling the amount of detail returned by the function value (see below)
If show = FALSE
, a numeric vector giving the values of
the reduced sample estimator.
If show=TRUE
, a list with three components which are
vectors of equal length,
Reduced sample estimate of the survival time c.d.f.
numerator of the reduced sample estimator
denominator of the reduced sample estimator
This function is needed mainly for internal use in spatstat, but may be useful in other applications where you want to form the reduced sample estimator from a huge dataset.
Suppose
If the number of observations cen
of all censoring times obs[k]
counts the number of values
(breaks[k],breaks[k+1]]
for [breaks[1],breaks[2]]
for nco
of all uncensored times,
i.e. those kaplan.meier
.
The return value rs
is the reduced-sample estimator
of the distribution function rs[k]
is the reduced sample estimate of F(breaks[k+1])
.
The value is exact, i.e. the use of histograms does not introduce any
approximation error.
Note that, for the results to be valid, either the histogram breaks
must span the censoring times, or the number of censoring times
that do not fall in a histogram cell must have been counted in
uppercen
.