Lspec: logspline estimation of a spectral distribution
Plots a spectral density function,
line spectrum, or spectral distribution from a model fitted with lspec
## S3 method for class 'lspec' plot(x, what = "b", n, add = FALSE, xlim, ylim, xlab = "", ylab = "", type, ...)
x |
|
what |
what should be plotted: b (spectral density and line spectrum superimposed), d (spectral density function), l (line spectrum) or p (spectral distribution function). |
n |
the number of equally spaced points at which to plot the fit; default is |
add |
indicate that the plot should be added to an existing plot. |
xlim |
X-axis plotting limits: default is c(0,pi), except when what = "p", when the default is c(-pi,pi). |
ylim |
Y-axis plotting limits. |
xlab,ylab |
axis labels. |
type |
plotting type; default is |
... |
all regular plotting options are passed on. |
If what = "p"
the plotting range cannot extend beyond the interval [-pi,pi].
Charles Kooperberg clk@fredhutch.org.
Charles Kooperberg, Charles J. Stone, and Young K. Truong (1995). Logspline Estimation of a Possibly Mixed Spectral Distribution. Journal of Time Series Analysis, 16, 359-388.
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.
data(co2) co2.detrend <- lm(co2~c(1:length(co2)))$residuals fit <- lspec(co2.detrend) plot(fit)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.