Suppose that we live in world where there is a single commodity called wheat, whose price at dateAttention conservation notice: Longish post on the details of a stochastic speculative pricing model from Economic Dynamics. It is my hope that working through this model will help me find a way forward in implementing a computational solution to the Kiyotaki-Moore (1997) model which will inform my extension of that model.

*t*is

*p*. Due to global warming, weather patterns in this world are quite volatile and as a result the harvest of wheat each period (i.e., production of wheat),

_{t}*W*, is IID random and drawn from a common density function φ.

_{t}Harvests take values in

*S≡[a,∞)*, where

*a>0*. The demand for wheat is made up by two groups: gatherers, and farmers. Gatherers' quantity demanded is

*Q*, where

_{G}=D(p)*p*is the price of fruit. It is assumed that the inverse demand function,

*D*

^{-1}(Q_{G})=p*(Q*,

_{G})*exists, is strictly decreasing and continuous, and satisfies*

*p*

*(Q*

_{G})*→∞*, as

*Q*

_{G}*→0*.

Farmers are risk neutral and have a special technology, called a barn, where they can store wheat between periods. If farmers purchase

*I*units of wheat on date

_{t}*t*, then they will be left with

*αI*, where

_{t}*α ∈ (0,1)*, next period. Risk-free interest rate is taken to be zero, so we can write the farmers' expected profits on

*I*units of wheat as:

_{t}*E*

_{t}(p_{t+1})αI_{t}- p_{t}I_{t}= (αE_{t}(p_{t+1}) - p_{t})I_{t}*αE*(1)

_{t}(p_{t+1}) - p_{t}≤ 0

There is also a logical constraint that tells us that if our farmers are maximizing expected profits, then if

Finally, market clearing condition requires that supply of wheat equal demand for wheat. Supply of wheat,

*αE*implies that

_{t}(p_{t+1}) - p_{t}< 0

*I*= 0

_{t}*(2)*

*X*, is sum of carryover from farmers,_{t}*αI*,_{t-1}*and the current harvest**W*, while demand is_{t}*D(p*_{t}) +*I*. Thus we get:_{t}*αI*+

_{t-1}*W*=

_{t}*X*=

_{t}*D(p*

_{t}) +*I*(3)

_{t}Take the initial condition for supply of wheat,

The question now becomes, how does one construct a system

I skip over the theory related to solving this system (the solution makes heavy use of Banach's fixed point theorem, consult Economic Dynamics for the details), and jump straight to the numerical solution.

Here is a plot of the pricing function when gatherers' inverse demand function is assumed to follow

This figure reproduces Figure 6.13 from Economic Dynamics. Once the pricing function has been calculated, it can be used to define a stochastic dynamical system as follows:

*X*as given._{0}∈ SThe question now becomes, how does one construct a system

*(I*_{t}, p_{t},*X*for investment, prices, and supply of wheat that satisfies (1-3)? The idea is to find a price function_{t})_{t ≥ 0}*p:S→(0,∞)*that depends*only*on the current state*X*(i.e.,_{t}*p*_{t}=*p(**X*for every_{t})*t*). The vector*(I*_{t}, p_{t},*X*for investment, prices, and supply of wheat would then evolve according to:_{t})_{t ≥ 0}*p*

_{t}=*p(*

*X*

_{t}),

*I*

_{t}*= X*-

_{t}*D(p*and

_{t}),*X*=

_{t+1}*αI*+

_{t}*W*(4)

_{t+1}For a given initial condition for supply of wheat,

*X*and a exogenous shock process_{0}∈ S*(**W*, the system (4) determines the time path of our vector_{t})_{t ≥ 1}*(I*_{t}, p_{t},*X*as a sequence of random variables. The idea then is to find a pricing function_{t})_{t ≥ 0}*p*such that (1) and (2) hold for the system (4).I skip over the theory related to solving this system (the solution makes heavy use of Banach's fixed point theorem, consult Economic Dynamics for the details), and jump straight to the numerical solution.

Here is a plot of the pricing function when gatherers' inverse demand function is assumed to follow

*D*

^{-1}(Q_{G})=p*(Q*=1/

_{G})*Q*and the random harvests

_{G}*W*

_{t}= a + c*B*where the

_{t}*B*are drawn from a beta distribution

_{t}*Beta(5,5)*.

*X*=

_{t+1}*αI(*

*X*+

_{t})*W*

_{t+1}where

I have simulated several trajectories, and marginal densities from this stochastic process.

*I(**x) = x - D(P*(x)).*I have simulated several trajectories, and marginal densities from this stochastic process.

In addition to marginal densities, I also plotted the density of the random harvest. The fact that the marginal densities of the supply of wheat converge to the density of the harvest process indicates that the farmers (who are wheat speculators) do not have any impact on the long-run supply of wheat.

For fun, I also reproduced the above analysis assuming that the harvest follows a Pareto distribution with α=2 (location and scale parameters are the same as those of the Beta distribution above). I find the plots amusing...

Note that, in this case, the marginal densities seem to converge...but not to the density of the harvest process. What does this mean? Well for this parametrization the Pareto distribution has a well defined mean, which implies that there is an "average" harvest (about 10 units in this case), but the variance is infinite! Thus while tomorrows harvest is likely to be close to 10, there is also a non-negligible probability that it could be 10,000. The mean is

*very much*influenced by these large harvests. The median of the distribution, however, is about 7.5. It is almost as if the speculators, because they wake up to below average harvests on most mornings, believe that the best thing to do is try to store wheat and sell it next period.Any other interpretations? The code is posted to my repository (you will want cpdynam.py and pSolve.py amongst other dependencies).

For those of you following my research...I am playing around with the idea of using the above methods to try to solve for the price of land in the

*Credit Cycles*model. The first thing I will need to endogenize the production side of the above economy. Any strong beliefs about whether this would be a useful way to proceed?

_{}
## No comments:

## Post a Comment