# Running Simple DSGE Model in Dynare III – Further Analysis

Last time we analyzed the steady-state of our model. We have developed equations which we can use to calculate the steady-state values for given parameters values. Today, we use these equations to see how could changes in parameters values affect the variables values. We can use the Octave (or Matlab) script which we developed last time. We just do minor updates. We select one variable and we assign a range of value instead of single value to this variable.

Let’s take parameter $\psi$. We would expect that lower values of this parameter would mean that the household would work more and also consume more because its willingness to work would increase. Let’s test our theory. We will examine an interval between 1.0 and 1.8. We will ask Octave to generate a vector of values between these two boundaries. We will to it by this command:

psi=1.0:0.001:1.8;

The value in the middle sets a size of space between these two values. In fact, this values defines the size of the vector. The value 0.001 is small enough to provide us a smooth plot. We also substitute all *, / and ^ operators with .*, ./ and .^. The standard operators without dots are used to operations with two numbers, two vectors or two matrices. But in our case, we want to multiply (or divide or power) all values of the vector by a single value. This is completely different operation than multiplying two vectors or two matrices. That’s why we use the operators with dots.

So this is our code.

alpha=0.5;
beta=0.99;
delta=0.05;
psi=1.0:0.001:1.8;
rho=0.9;
sigma=0.01;
A=1;

X=((beta.*A.*alpha)./(beta.*(delta-1)+1)).^(1./(1-alpha));
c=((A.*(beta.*A.*alpha)./(beta.*(delta-1)+1)).^(alpha./(1-alpha)).*(1-alpha))./psi ;
n=(-c)./(X-(1-delta).*X-A.*X.^alpha);
k=X.*n;
subplot (3, 1, 1);
plot(psi,c);
xlabel("psi");
ylabel("consumption");
subplot (3, 1, 2);
plot(psi,k);
xlabel("psi");
ylabel("kapital");
subplot (3, 1, 3);
plot(psi,n);
xlabel("psi");
ylabel("labor");

We use the subplot command to generate three plots in one window – one plot for each variable. As we can see on the image below, our assumption was right. Higher values of the $\psi$ parameter result into higher consumption and higher labor supply. It also means that there will be more capital in the economy.

Let’s analyze one more parameter. I have selected $\alpha$. If we increase the value of this variable, we would increase the marginal productivity of capital and decrease the marginal productivity of labor. Consumers will probably react by increasing its savings and decreasing the time spent in work. Is this statement true? We can check it on the picture below.

So far, we have been changing the value of one parameter, leaving the other parameters constant. But can we change values simultaneously? Yes, we can? The result will be a three-dimensional plot. We will use the mesh command to render the plot. We also need the meshgrid command to generate a matrix of values which is a combination of all values from both of the vectors. In this case, the “dot operators” perform the operation on all values of the matrix.

We will change values of $\alpha$ and $\psi$.This is the final code. I rotate the plot by R command before saving it.

talpha = linspace (0.9, 0.5, 200)';
tpsi = linspace (1.2, 1.8, 200)';
[alpha, psi] = meshgrid (talpha, tpsi);

beta=0.99;
delta=0.05;
A=1;

X=((beta.*A.*alpha)./(beta.*(delta-1)+1)).^(1./(1-alpha));
c=((A.*(beta.*A.*alpha)./(beta.*(delta-1)+1)).^(alpha./(1-alpha)).*(1-alpha))./psi ;
n=(-c)./(X-(1-delta).*X-A.*X.^alpha);
k=X.*n;

mesh(talpha, tpsi, n);
xlabel("alpha");
ylabel("psi");
zlabel("labor");

Here you can see the result. As you can see, decreasing of both parameters simultaneously causes the household to work less.

Let’s try one more combination. What would do a change of $\alpha$ and $\delta$ with the capital stock? As we can see from the plot, lowering the  $\alpha$ and increasing the $\delta$ causes the capital stock to decrease.