## bsplinePsd 0.6.0

Have included a scaling/rescaling factor within the code to remove
numerical errors when data is close to machine precision.

Have optimised the starting values for the Dirichlet process
parameters.

## bsplinePsd 0.5.0

Previous versions only allowed the user to use cubic B-spline
densities. This version allows the user to choose any degree.
Normalising the B-splines now uses a trivial integral formula.

The function gibbs_bspline can now handle odd length time series.

An S3 plot method has been included so the user can easily plot their
PSD estimate.

## bsplinePsd 0.2.0

Added an argument called k1 in the gibbs_bspline function. This
allows the user to specify the starting value for parameter k. If
well-chosen, this can speed up convergence significantly. The default is
set to 20, which works well on all of the cases I have come across. If
missing (NA), then a random integer between degree + 2 and kmax will be
selected as the starting value for k.