jackknife: Compute jackknife pseudo values.

Description Usage Arguments Details Note Author(s) References See Also Examples

Description

Compute jackknife pseudo values.

Usage

1
jackknife(object, times, cause, keepResponse = FALSE, ...)

Arguments

object

Object of class "prodlim".

times

Time points at which to compute pseudo values.

cause

For competing risks the cause of failure.

keepResponse

If TRUE add the model response, i.e. event time, event status, etc. to the result.

...

not used

Details

Compute jackknife pseudo values based on marginal Kaplan-Meier estimate of survival, or based on marginal Aalen-Johansen estimate of cumulative incidence.

Note

The R-package pseudo does a similar job, and appears to be a little faster in small samples, but much slower in large samples. See examples.

Author(s)

Thomas Alexander Gerds <tag@biostat.ku.dk>

References

Andersen PK & Perme MP (2010). Pseudo-observations in survival analysis Statistical Methods in Medical Research, 19(1), 71-99.

See Also

prodlim

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
## pseudo-values for survival models

d=SimSurv(20) 
f=prodlim(Hist(time,status)~1,data=d) 
jackknife(f,times=c(3,5))

## in some situations it may be useful to attach the
## the event time history to the result
jackknife(f,times=c(3,5),keepResponse=TRUE)

# pseudo-values for competing risk models
d=SimCompRisk(10) 
f=prodlim(Hist(time,event)~1,data=d) 
jackknife(f,times=c(3,10),cause=1)
jackknife(f,times=c(3,10,17),cause=2)


Search within the prodlim package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.