Dear All, I'm in trouble with state space models, as I find them difficult to implement. I was able to set up a few with just one time-variant unobservable variable but I'm not sure on what I'm doing when there are more. Let's say we depart from the example of the Phillips curve on pages 345-346 of Gretl manual, where the inflation rate (INFQ) depends on the unemployment rate (URX). In the example, the intercept is time-invariant and the coefficient for URX follows a random walk. Let's expand the model such that there is a second independent variable Effective Exchange Rate (EER) explaining INFQ. Additionally, I want all three coefficients to follow an AR(1) process. The image attached shows this new model. How hard would that be to implement these changes on the script depicted on manual page 347 (as follows):

On Fri, 8 Oct 2021, Sven Schreiber wrote: > As explained below equation (36.5), what the listing 36.4 does is to > treat the time-varying coefficent as the latent state, using the > (apparently implicit) default of F == 1 == statemat. > (@jack: It seems that that default for statemat is not documented. If > there is no such default, then I don't know what's going on.) In that example, F is set to 1 because it is passed as an argument to ksetup(), and is never changed.

On Fri, 8 Oct 2021, Sven Schreiber wrote: > Am 08.10.2021 um 13:44 schrieb Riccardo (Jack) Lucchetti:

>> In that example, F is set to 1 because it is passed as an argument to >> ksetup(), and is never changed. > > Ah of course, that makes sense!

> One follow-up question though: What about the line with B.depvar = > {INFQ}? A search for "depvar" in the state-space chapter doesn't turn up > anything useful. And isn't the first arg to ksetup() already specifying > that? A Kalman bundle can have extra keys, apart from the "special" ones such as statevar, obsymat, etc. In that example depvar is used to have an object you can subtract b0 from, just for convenience.

On Thu, 7 Oct 2021, F.R.Costa wrote: > Dear All, > > I'm in trouble with state space models, as I find them difficult to implement. I was able to set up a few with just one time-variant unobservable variable but I'm not sure on what I'm doing when there are more. Let's say we depart from the example of the Phillips curve on pages 345-346 of Gretl manual, where the inflation rate (INFQ) depends on the unemployment rate (URX). In the example, the intercept is time-invariant and the coefficient for URX follows a random walk. > > Let's expand the model such that there is a second independent variable Effective Exchange Rate (EER) explaining INFQ. Additionally, I want all three coefficients to follow an AR(1) process. The image attached shows this new model. How hard would that be to implement these changes on the script depicted on manual page 347 (as follows): > Speaking from personal experience, the state-space approach to time-varying parameter models is often difficult in practice because for large-ish models the identification issues become quite relevant and maximising the loglikelihood is quite difficult (although the EM algorithm may help). There's an interesting non-parametric alternative that I find myself using quite often recently that was put forward by Giraitis, Kapetanios and several co-authors and is implemented in the gretl "ketvals" package. I'm attaching an example script which shows the two alternatives on some simulated data, where the simulated coefficients are smoothed versions of AR(1) processes. As you can see, both approaches reconstruct the histories for both coefficients relatively well, but the non-parametric approach is much faster. ------------------------------------------------------- Riccardo (Jack) Lucchetti Dipartimento di Scienze Economiche e Sociali (DiSES) Università Politecnica delle Marche (formerly known as Università di Ancona) r.lucchetti(a)univpm.it http://www2.econ.univpm.it/servizi/hpp/lucchetti -------------------------------------------------------

Are there other possible reserved Kalman keys which do not show up in this particular example, or is this it?

On Thu, 7 Oct 2021, F.R.Costa wrote: > Dear All, > > I'm in trouble with state space models, as I find them difficult to implement. I was able to set up a few with just one time-variant unobservable variable but I'm not sure on what I'm doing when there are more. Let's say we depart from the example of the Phillips curve on pages 345-346 of Gretl manual, where the inflation rate (INFQ) depends on the unemployment rate (URX). In the example, the intercept is time-invariant and the coefficient for URX follows a random walk. > > Let's expand the model such that there is a second independent variable Effective Exchange Rate (EER) explaining INFQ. Additionally, I want all three coefficients to follow an AR(1) process. The image attached shows this new model. How hard would that be to implement these changes on the script depicted on manual page 347 (as follows): > Speaking from personal experience, the state-space approach to time-varying parameter models is often difficult in practice because for large-ish models the identification issues become quite relevant and maximising the loglikelihood is quite difficult (although the EM algorithm may help). There's an interesting non-parametric alternative that I find myself using quite often recently that was put forward by Giraitis, Kapetanios and several co-authors and is implemented in the gretl "ketvals" package. I'm attaching an example script which shows the two alternatives on some simulated data, where the simulated coefficients are smoothed versions of AR(1) processes. As you can see, both approaches reconstruct the histories for both coefficients relatively well, but the non-parametric approach is much faster. ------------------------------------------------------- Riccardo (Jack) Lucchetti Dipartimento di Scienze Economiche e Sociali (DiSES) Università Politecnica delle Marche (formerly known as Università di Ancona) r.lucchetti(a)univpm.it http://www2.econ.univpm.it/servizi/hpp/lucchetti -------------------------------------------------------