Plots an oldlogspline density, distribution function, hazard
function or survival function
from
a logspline density that was fitted using
the 1992 knot deletion algorithm.
The 1997 algorithm using knot
deletion and addition is available using the logspline function.
Usage
## S3 method for class 'oldlogspline'
plot(x, n = 100, what = "d", xlim, xlab = "", ylab = "",
type = "l", add = FALSE, ...)
Arguments
x
logspline object, typically the result of logspline.
n
the number of equally spaced points at which to plot the density.
what
what should be plotted:
"d" (density), "p" (distribution function), "s" (survival
function) or "h" (hazard function).
xlim
range of data on which to plot. Default is from the 1th to the 99th percentile of
the density, extended by 10% on each end.
xlab,ylab
labels plotted on the axes.
type
type of plot.
add
should the plot be added to an existing plot.
...
other plotting options, as desired
Details
This function produces a plot of a oldlogspline fit at n equally
spaced points roughly covering the support of the density. (Use
xlim=c(from,to) to change the range of these points.)
Charles Kooperberg and Charles J. Stone. Logspline density estimation
for censored data (1992). Journal of Computational and Graphical
Statistics, 1, 301–328.
Charles J. Stone, Mark Hansen, Charles Kooperberg, and Young K. Truong.
The use of polynomial splines and their tensor products in extended
linear modeling (with discussion) (1997). Annals of Statistics,
25, 1371–1470.