COEFbothscale(l, plotclustonly = FALSE, ...)
findstysols
,
which contains the results of an optimization
to find solutions that correspond to stationary
series which are the time-varying linear combination
of two locally sTRUE
then only produce the hierarchical
clustering plot.epscale
and epclust
respectively.findstysols
runs multiple optimization
runs (it can do only one is required) to try and find (multiple)
solutions to the problem of finding costationary time series.
It does this by running the optimization from multiple random
starts. The solutions (if they exist) to these optimizations
are not necessarily unique, and it is useful to find as many
genuine different solutions as possible. However, even if there
was one unique (mathematical) solution, due to different random
starts, the vagaries of the optimizer, and the fact that the
hypothesis test for stationarity is not a deterministic test,
it is usually the case
the many solutions very close to, but not identical to, the
unique solutions might be found. However, they would morally
be thought to be the same solution.In the case of multiple actual solutions, then the optimization algorithm might find many solutions close to the real ones and, of course, it might completely miss a real one (or not be near to a real one). However, the situation often arises where, after multiple starts, one ends up with clusters of computed solutions that group themselves around the `real' solution.
In other words, the solution vectors can be viewed as a multivariate data set where the cases correspond to the results of different optimization runs and the variables correspond to the coefficients of the time-varying linear combinations.
Both multidimensional scaling (cmdscale
) and
hiearchical clustering (hclust
) are used to determine
possible clusterings of solutions. Then, representative members
from these clusters can be further investigated with a function
such as LCTSres
COEFscale
,findstysols
, LCTSres
#
# To be run directly after the example in findstysols
#
COEFbothscale(tmp)
Run the code above in your browser using DataLab