plot.psd {bsplinePsd} | R Documentation |
This function plots the log periodogram, log posterior median PSD, and log 90% credible region PSD. The x-axis uses angular frequency and the y-axis is plotted on the log scale. The PSD at the zero frequency is removed from the plot. If the time series is even length, the PSD at the last frequency is also removed from the plot.
## S3 method for class 'psd' plot(x, legend.loc = "topright", ylog = TRUE, ...)
x |
an object of class psd |
legend.loc |
location of legend out of "topright" (default), "topleft", "bottomright", "bottomleft". If set to NA then no legend will be produced |
ylog |
logical value (default is TRUE) to determine if PSD (y-axis) should be on natural log scale |
... |
other graphical parameters from the plot.default function |
plot of the estimate of the (log) PSD
## Not run: set.seed(12345) # Simulate AR(4) data n = 2 ^ 7 ar.ex = c(0.9, -0.9, 0.9, -0.9) data = arima.sim(n, model = list(ar = ar.ex)) data = data - mean(data) # Run MCMC with linear B-spline prior (may take some time) mcmc = gibbs_bspline(data, 10000, 5000, degree = 1) # Plot result plot(mcmc) # Plot result on original scale with title plot(mcmc, ylog = FALSE, main = "Estimate of PSD using the linear B-spline prior") ## End(Not run)