Take the 2-minute tour ×
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It's 100% free, no registration required.

I am doing time series data analysis by state space methods. With my data the stochastic local level model totally outperformed the deterministic one. But the deterministic level and slope model gives better results than with stochastic level and stochastic/deterministic slope. Is this something usual? All methods in R require initial values, and I read somewhere that fitting an ARIMA model first and taking values from there as initial values for state space analysis is one way; possible? or any other proposition? I should confess here that I'm totally new to state space analysis.

share|improve this question
Please provide the example. Now it is unclear what is your actual problem. –  mpiktas Feb 28 '11 at 11:08
Do you mean an exponential smoothing state space model? What R packages are you using? –  Zach May 7 '11 at 22:27
are you trying to compare models, or do you want to select a model? –  naught101 Apr 29 '12 at 9:10
Firstly, as already mentioned, it is not clear what your actual problem is. You write that A outperformed B and B gives better results than A. This is confusing. Secondly, the "forecast" R package has some automatic time series methods. They include: auto.arima(), ets(), tbats() and bats(). –  power Nov 30 '12 at 9:55
Can you explain what you mean when you say 'outperformed' and 'gives better results than'? –  Glen_b Dec 25 '12 at 22:56
show 1 more comment

1 Answer

To answer your first question. Yes, all is possible. It is not usual or unusual. You should let the data tell you what the correct model is. Try augmenting the model further with seasonals, cycles, and explanatory regressors if possible.

You should not only be comparing the Akaike Information Criterion (AIC) to compare models, but also checking to see that residuals (irregular term) are normal, homoskedastic, and independent (Ljung-Box test). If you can find a model that has all of these desirable properties. This should be your preferred model (it is likely that a model with all these properties will have the best AIC).

Although the initial values will affect which maximum point of the log-likelihood function is found, if your model is well specified, it shouldn't vary too much and there should be an obvious candidate for the best model with the best initial values. I do a lot of this type of analysis in Matlab and I found the best way to find initial values is just to play around for a bit. It can be tedious but it works out well in the end.

share|improve this answer
add comment

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.