Dear Michel
Perhaps this question is not well thought out, but…
I want to use Bayesian estimation to deal with model with structural breaks. Namely, I used to have the following model (simplified version):
model(linear);
y = gamma*y(-1) + delta*x + e;
x = alpha*x(-1) + beta*y + u;
end;
y and x are observable
Now, I believe that x has a structural break at time t0, I add one extra data series which is zero untill t0 and one ater that, call it It. In description I list it among exogenous observable variables
the model becomes
model;
y = gamma*y(-1) + (delta+z_d*It)*x + e;
x = (alpha+z_a*It)*x(-1) + (beta+z_b*It)*y + u;
end;
I obviously add three new parameters and do not call model linear any more.
If I run dynare, it checks steady state, plots priors and then says that “inner matrix dimensions must agree”. Is there anything wrong with this definition? How can I fix it?
PS. I think I know: the problem is with definition of a product. Clearly, I need a “dot product”. The dynare does not seem to recognise it, if I simply write
model;
y = gammay(-1) + deltax+z_d*It.*x+ e;
x = alphax(-1)+z_aIt.x(-1) + betay+z_b*It.*y + u;
end;
I also tried to define separately xt = x.*It; but it did not help. I am afraid I cannot make x.*It as a new observable variable as I need to multiply It on unobservable variables too. Is the any way to explain DYNARE what it is?