psh_fit: Peña-Strawderman-Hollander survival estimator
In survrec: Survival Analysis for Recurrent Event Data

psh_fit

R Documentation

Peña-Strawderman-Hollander survival estimator

Description

Estimates the survival function of the inter-occurrence times of recurrent event data with the generalized product-limit estimator (PLE) of Peña, Strawderman and Hollander (2001).

Usage

psh_fit(x, tvals)

Arguments

`x`	a survival recurrent event object, see `Survr()`.
`tvals`	optional vector of times at which the survival function is also evaluated.

Details

The estimator computed by this function is the nonparametric estimator of the inter-event time survivor function under the assumption of a renewal or IID model. It generalizes the product-limit estimator to the situation where the event is recurrent. For details and the theory behind the estimator, please refer to Peña, Strawderman and Hollander (2001).

Value

A list with class survfitr containing:

n: number of units or subjects observed.
m: number of events of each subject.
failed: uncensored gap times, ordered by subject.
censored: censored gap time of each subject.
time: ordered distinct event times.
n.event: number of events at each distinct time.
AtRisk: matrix of units at risk at each distinct time (rows are subjects).
survfunc: survival estimate at each distinct time.
std.error: standard error of the survival estimate.
tvals: copy of the tvals argument.
PSHpleAttvals: survival estimate at the tvals times.

References

Peña, E.A., Strawderman, R. and Hollander, M. (2001). Nonparametric Estimation with Recurrent Event Data. Journal of the American Statistical Association 96, 1299–1315.

Examples

data(MMC)
fit <- psh_fit(Survr(MMC$id, MMC$time, MMC$event))
fit
plot(fit, conf.int = FALSE)

# compare with the gamma frailty MLE
fit <- mlefrailty_fit(Survr(MMC$id, MMC$time, MMC$event))
lines(fit, lty = 2)

# and with Wang-Chang
fit <- wc_fit(Survr(MMC$id, MMC$time, MMC$event))
lines(fit, lty = 3)

survrec documentation built on Sept. 22, 2026, 9:07 a.m.