<p>One way to expand the model is to allow for different income classes as stated in assumption (1) below. Let the model be completely standard with all the usual assumption of the monocentric city model as laid out in Brueckner (1987) the paper you cite.</p>
<p>I will only present the basic idea of the expanded model. Hence I do not give proof of existence of equilibria, proof of uniqueness, nor a comparative static analysis or numerical procedures to solve for equilibrium given values of exogenous variables.</p>
<p>I will allow myself to assume Cobb-Douglas preferences with residents utility given as
<span class="math-container">$$U(c,h) = A\left(\frac{c}{1-\alpha}\right)^{1-\alpha}\left(\frac{h}{\alpha}\right)^\alpha,$$</span>
where <span class="math-container">$c$</span> is composite good and <span class="math-container">$h$</span> is square meters of housing.</p>
<p>It is well known that the Marshall demands are <span class="math-container">$c^\star=(1-\alpha)I/p_c$</span> and <span class="math-container">$h^\star=\alpha I/p_h$</span> with <span class="math-container">$I$</span> being income. Plugged into utility function results in the value function</p>
<p><span class="math-container">$$V(p_c,p_h,I) = \frac{AI}{p_h^\alpha p_c^{1-\alpha}} = \frac{I}{p_h^\alpha},$$</span></p>
<p>having assumed that <span class="math-container">$p_c=1$</span> and <span class="math-container">$A=1$</span>. The value function tells us that the max price a consumer can pay for a unit of housing <span class="math-container">$h$</span> if she is to attain the utility level <span class="math-container">$u$</span> is</p>
<p><span class="math-container">$$\frac{I}{p_h^\alpha} = u \Leftrightarrow p(I,u) = \left(\frac{I}{u}\right)^\frac{1}{\alpha},$$</span>
found by inverting the value function to express price as a function of income and utility (which offcourse is the well known <strong>bid-function</strong>).</p>
<p>The next step is to find the <strong>bid-lot size</strong> which is the amount of housing consumed by a resident attaining the utility level <span class="math-container">$u$</span> having income <span class="math-container">$I$</span>. To derive this I will simply use that the Marshall demand <span class="math-container">$h^\star_M(p,I)$</span> with bid-function <span class="math-container">$p(I,u)$</span> substituted for price is the bid-lot size <span class="math-container">$h_b(u,I)$</span> [same holds offcourse for the Hicks demand]. Hence the bid lot size is given as</p>
<p><span class="math-container">$$h^\star_M(p(I,u),I) = h_b(I,u) = \frac{\alpha I}{p(I,u)} = \alpha I^{1-1/\alpha}u^{1/\alpha},$$</span></p>
<p>which is decreasing in income because when consumer gets higher income and utility attained is held constant the price offered for housing is increase and hence less housing consumed. Importantly this is not a description of equilibrium behavior, what is being said is that if the consumer has higher income she <em>can</em> offer a higher price for a unit of housing while retaining utility level <span class="math-container">$u$</span> but that will require her to substitute away from housing over to consumption of <span class="math-container">$c$</span>.</p>
<p>I now add the linear geometry and linear transport costs of standard monocentric city expositions. Hence the city expands from center at <span class="math-container">$d=0$</span> to potentially infinity <span class="math-container">$d \in [0,\infty)$</span> and transport costs are <span class="math-container">$t(d) = t d$</span> such that income net transport is <span class="math-container">$I(d) = w-td$</span>.</p>
<p>I now introduce the extending assumption that</p>
<p><strong>Assumption (1):</strong> There are two income classes <span class="math-container">$j=\{1,2\}$</span> with incomes exogenously given <span class="math-container">$w_2&gt;w_1$</span> and size of population income segments exogenously given as <span class="math-container">$N_1$</span> and <span class="math-container">$N_2$</span>.</p>
<p>and for ease of exposition I will let housing be exogenously</p>
<p><strong>Assumption (2):</strong> At location <span class="math-container">$d$</span> the housing is given as <span class="math-container">$H(d)$</span>,</p>
<p>the second assumption is not so restrictive as it may appear (but that is easier to illustrate later).</p>
<p>Given the wages I have income net transport for the two classes <span class="math-container">$I_1(d) = w_1 - td$</span> and <span class="math-container">$I_2(d)=w_2 - td$</span> which are plugged into bid function <span class="math-container">$p(u,I)$</span> and bid lot size <span class="math-container">$h_m(u,I)$</span> to get the main result that</p>
<p><span class="math-container">$$(3)\ \ \ p_j(u_j,d) = \left(\frac{w_j - td}{u_j}\right)^{1/\alpha} \ \ \ \ h_j^m(u_j,d) =\alpha (w_j - td)^{1-1/\alpha}u_j^{1/\alpha},$$</span></p>
<p>which is really not new in the sense that it is completely like the model with homogenous residents as in the Brueckner (1987) model just adding Cobb-Douglas preferences.</p>
<p>To illustrate the use of these functions select parameter values for <span class="math-container">$(w_1,w_2,\alpha,t)$</span>. Assume existence of equilibrium with equilibrum utilities <span class="math-container">$u^\star_1$</span> and <span class="math-container">$u_2^\star$</span> and select values for these (<span class="math-container">$u_2^\star &gt; u_1^\star$</span> rich get higher utlity in equilibrium). Also select a value for the alternative use of land <span class="math-container">$r_a$</span> referred to as the land rent (not to be confused by the rent paid for land by producers of housing if that sector was endogenized in the model).</p>
<p>In the figure below the bid function <span class="math-container">$p_j(u_j,d)$</span> of the high and low income are illustrated. Residents will live where they bid the highest value for housing. The city expands to where the bid-rents of either income class is higher than the land rent.</p>
<blockquote>
<blockquote>
<p><a href="https://i.sstatic.net/LbePz.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/LbePz.jpg" alt="monocentric city two income classes" /></a></p>
</blockquote>
</blockquote>
<p>Two observations are made:</p>
<p><strong>Observation (1): The low income resident will live in the center of the city (this requires housing to be normal good)</strong></p>
<p><strong>Observation (2): There are perfect income stratification with low income in one area and high income in another area</strong></p>
<p>To close the model the housing market must be in equilibrium. Let <span class="math-container">$[0,\hat d]$</span> be the locations where the low income reside. Then for <span class="math-container">$d \in [0,\hat d]$</span> housing demand is given as <span class="math-container">$h_{j}^m(u_{j},d)$</span> per person and units of housing is <span class="math-container">$H(d)$</span> so population density at <span class="math-container">$d$</span> must be <span class="math-container">$H(d)/h_{j}^m(u_{j},d)$</span> and integrating over <span class="math-container">$[0,\hat d]$</span> this must be equal to <span class="math-container">$N_1$</span> the number of low income residents I have exogenously stipulated to live in the city. So equilibrium conditions includes that</p>
<p><span class="math-container">$$N_1 = \int_0^{\hat d} \frac{H(t)}{h_{1}^m(u_{1},t)} dt$$</span></p>
<p>and</p>
<p><span class="math-container">$$N_2 = \int_{\hat d}^{\bar d} \frac{H(t)}{h_{2}^m(u_{2},t)} dt.$$</span></p>
<p>As is hopefully apparent from the figure above it is the case that: Given two utilities <span class="math-container">$(u_1,u_2)$</span> the bid functions determine <span class="math-container">$\hat d$</span> and <span class="math-container">$\bar d$</span> - where <span class="math-container">$\bar d$</span> is the value of <span class="math-container">$d$</span> where <span class="math-container">$p_2(u_2,d) = r_a$</span> - and given these values the integrals can be calculated to get <span class="math-container">$N_1(u_1,u_2)$</span> and <span class="math-container">$N_2(u_1,u_2)$</span>. When the equilibrium is unique these equations</p>
<p><span class="math-container">$$N_1 = N_1(u_1,u_2) \ \ \ N_2=N_2(u_1,u_2),$$</span></p>
<p>will be invertible such that equilibrium utilities can be found from the exogenously given <span class="math-container">$N_1$</span> and <span class="math-container">$N_2$</span>.</p>
<p>Proof of existence and uniqueness of equilibrium for a model with <span class="math-container">$K$</span> income classes and standard restrictions on the utility functions (not assuming Cobb Douglas) are given in</p>
<p>John Hartwick (1976) <em>Comparative Statics of a Residential Economy with
Several Classes</em> Journal of Economic Theory 13, 396-413.</p>
<p>which also include comparative statics of the model.</p>
<p>Proof of existence and uniqueness of equilibrium for a model is also given in</p>
<p>Fujita (1989) <em>Urban economic theory: Land use and city size</em></p>
<p>which also offers an algorithm for solving for the equilibrium.</p>
<p>See also</p>
<p>Puga and Duranton (2015) <em>Urban Land Use</em></p>
<p>for discussion of main results of this and other extensions and empirical relevance of the models (but not detailed exposition).</p>
<p><strong>Note on housing supply:</strong> Instead of assuming housing to be exogenously given I could have followed Brueckner (1989) and assumed land to exogenously given <span class="math-container">$L(d)$</span> and then housing be produced by constant returns to scale production under perfect competition using capital and land with the capital supply being infinitely price elastic and having price <span class="math-container">$r_K$</span> in which case housing supplied is a function of <span class="math-container">$d$</span> through <span class="math-container">$p(d)$</span>. So housing supply is a function of the price <span class="math-container">$H(p)$</span> but at any distance <span class="math-container">$d$</span> the price residents are willing to give is set by the bid function which is a function of <span class="math-container">$d$</span> hence either way we get <span class="math-container">$H(d)$</span>.</p>
<p>Here is some really fast written R-code implementing the solution algorithm in Fujita (1989):</p>
<pre><code>library(rootSolve)

a &lt;- 0.35
t &lt;- 200
w &lt;- c(6000,12000)
N &lt;- c(500,500)
J &lt;- 2
x &lt;- rep(1,J+1)

u &lt;- rep(170,J)


# Function to make bid-max lot size for agent of type j
# with income and utility (w_j,u_j). 
# Bid-max lot size is a function of distance to center d
# given type specifics (w_j,u_j).

make_g &lt;- function(u_j,w_j,t,a)
    {
        g &lt;- function(d)
            {
                nominator &lt;- a*u_j^(1/a)
                denominator &lt;- (w_j-t*d)^((1-a)/a)*.6*pi*d
                bid_lot &lt;- nominator/denominator
                # return should be 1/bid_lot
                # but denominator can be 0 so do not calculate it
                # and use it as denominator instead 
                denominator &lt;- ifelse((w_j-t*d)&gt;0,denominator,0)
                return(denominator/nominator)
            }
        return(g)
    }


# Solve the integral of g with respect to r from x_L = x_j to x_U = x_j+1
# such that the integral equals N_j 

# Initiate with j=1 so x_l is set to 0 and u_j is some guess
j &lt;- 1
x_L &lt;- x[j]
u_j &lt;- u[j]
    
    # Get exogenous income and population size for relevant segment
    # of population.
    w_j &lt;- w[j]
    N_j &lt;- N[j]

    g &lt;- make_g(u_j,w_j,t,a)
    v &lt;- function(x_U) {return(N_j-integrate(g,x_L,x_U)<span class="math-container">$value)} 
    x[j+1] &lt;- uniroot(v,c(x_L+0.0001,500))$</span>root

# Initiate j=2
j &lt;- 2
    
    # Find u_2 given I_1,I_2 at x and u_1
    u[j] &lt;- u[j-1]*(w[j] - t*x[j])/(w[j-1]-t*x[j])

x_L &lt;- x[j]
u_j &lt;- u[j]

    # Get exogenous income and population size for relevant segment
    # of population.
    w_j &lt;- w[j]
    N_j &lt;- N[j]

    g &lt;- make_g(u_j,w_j,t,a)
    v &lt;- function(x_U) {return(N_j-integrate(g,x_L,x_U)<span class="math-container">$value)} 
    x[j+1] &lt;- uniroot(v,c(x_L+0.0001,500))$</span>root



# Calculate price at x[j+1]
j &lt;- 3  
p_edge &lt;- ((w[j-1] - t*x[j])/u[j-1])^(1/a)
x
u
p_edge

p_max &lt;- (w[1]/u[1])^(1/a)+10
plot(x[1]:(x[J+1]+1),type=&quot;n&quot;,ylim=c(0,p_max),ylab=&quot;housing price&quot;,xlab=&quot;distance to city center&quot;)
colors &lt;- c(&quot;blue&quot;,&quot;red&quot;)
for (j in 1:J)
    {
        x_L &lt;- x[j]
        x_U &lt;- x[j+1]

        d &lt;- seq(0,x[J+1],length.out=100)
        u_j &lt;- u[j] 
        p &lt;- ((w[j] - t*d)/u_j)^(1/a)
        points(d,p,type=&quot;l&quot;,col=colors[j])
    }

abline(v=x[2],col=&quot;grey&quot;)
abline(h=p_edge,col=&quot;black&quot;)
legend(12,20000,c(&quot;bid function (low income)&quot;,&quot;bid function (high income)&quot;,&quot;sorting threshold&quot;,
&quot;land rent&quot;),
col=c(&quot;blue&quot;,&quot;red&quot;,&quot;grey&quot;,&quot;black&quot;),cex=0.7,bty=&quot;n&quot;,pch=24)
#savePlot(&quot;monocentric.jpg&quot;,type=&quot;jpg&quot;)
</code></pre>

CC BY-SA 4.0
2021-05-03T14:30:37.780Z
2021-05-04T18:51:12.853Z
2021-05-04T18:51:12.853Z
Jesper Hybel
<p>I am exploring the generalized treatment of the Muth-Mills monocentric city model covered in Brueckner's <a href="http://www.socsci.uci.edu/%7Ejkbrueck/course%20readings/handbook%20chapter.pdf" rel="nofollow noreferrer">article</a>. The model is based on an assumption of homogenous income levels across the city. However, a more realistic assumption would be to assume the presence of different income groups, which theoretically results in different bid-rent curves. How might I modify the model given in the article to reflect this?</p>

<urban-economics>
Extending the Muth-Mills Model to Different Income Groups?
kékszajkók
# Extending the Muth-Mills Model to Different Income Groups?

(Asked by: username_0 on 2021-04-29)

I am exploring the generalized treatment of the Muth-Mills monocentric city model covered in Brueckner's [article](http://www.socsci.uci.edu/%7Ejkbrueck/course%20readings/handbook%20chapter.pdf). The model is based on an assumption of homogenous income levels across the city. However, a more realistic assumption would be to assume the presence of different income groups, which theoretically results in different bid-rent curves. How might I modify the model given in the article to reflect this?

**Tags:** urban-economics

---

## 1 answers

### Answer by username_1 on 2021-05-03. Score: 6 | Accepted answer

One way to expand the model is to allow for different income classes as stated in assumption (1) below. Let the model be completely standard with all the usual assumption of the monocentric city model as laid out in Brueckner (1987) the paper you cite.
I will only present the basic idea of the expanded model. Hence I do not give proof of existence of equilibria, proof of uniqueness, nor a comparative static analysis or numerical procedures to solve for equilibrium given values of exogenous variables.
I will allow myself to assume Cobb-Douglas preferences with residents utility given as
$$U(c,h) = A\left(\frac{c}{1-\alpha}\right)^{1-\alpha}\left(\frac{h}{\alpha}\right)^\alpha,$$
where $c$ is composite good and $h$ is square meters of housing.
It is well known that the Marshall demands are $c^\star=(1-\alpha)I/p_c$ and $h^\star=\alpha I/p_h$ with $I$ being income. Plugged into utility function results in the value function
$$V(p_c,p_h,I) = \frac{AI}{p_h^\alpha p_c^{1-\alpha}} = \frac{I}{p_h^\alpha},$$
having assumed that $p_c=1$ and $A=1$. The value function tells us that the max price a consumer can pay for a unit of housing $h$ if she is to attain the utility level $u$ is
$$\frac{I}{p_h^\alpha} = u \Leftrightarrow p(I,u) = \left(\frac{I}{u}\right)^\frac{1}{\alpha},$$
found by inverting the value function to express price as a function of income and utility (which offcourse is the well known bid-function).
The next step is to find the bid-lot size which is the amount of housing consumed by a resident attaining the utility level $u$ having income $I$. To derive this I will simply use that the Marshall demand $h^\star_M(p,I)$ with bid-function $p(I,u)$ substituted for price is the bid-lot size $h_b(u,I)$ [same holds offcourse for the Hicks demand]. Hence the bid lot size is given as
$$h^\star_M(p(I,u),I) = h_b(I,u) = \frac{\alpha I}{p(I,u)} = \alpha I^{1-1/\alpha}u^{1/\alpha},$$
which is decreasing in income because when consumer gets higher income and utility attained is held constant the price offered for housing is increase and hence less housing consumed. Importantly this is not a description of equilibrium behavior, what is being said is that if the consumer has higher income she can offer a higher price for a unit of housing while retaining utility level $u$ but that will require her to substitute away from housing over to consumption of $c$.
I now add the linear geometry and linear transport costs of standard monocentric city expositions. Hence the city expands from center at $d=0$ to potentially infinity $d \in [0,\infty)$ and transport costs are $t(d) = t d$ such that income net transport is $I(d) = w-td$.
I now introduce the extending assumption that
Assumption (1): There are two income classes $j=\{1,2\}$ with incomes exogenously given $w_2>w_1$ and size of population income segments exogenously given as $N_1$ and $N_2$.
and for ease of exposition I will let housing be exogenously
Assumption (2): At location $d$ the housing is given as $H(d)$,
the second assumption is not so restrictive as it may appear (but that is easier to illustrate later).
Given the wages I have income net transport for the two classes $I_1(d) = w_1 - td$ and $I_2(d)=w_2 - td$ which are plugged into bid function $p(u,I)$ and bid lot size $h_m(u,I)$ to get the main result that
$$(3)\ \ \ p_j(u_j,d) = \left(\frac{w_j - td}{u_j}\right)^{1/\alpha} \ \ \ \ h_j^m(u_j,d) =\alpha (w_j - td)^{1-1/\alpha}u_j^{1/\alpha},$$
which is really not new in the sense that it is completely like the model with homogenous residents as in the Brueckner (1987) model just adding Cobb-Douglas preferences.
To illustrate the use of these functions select parameter values for $(w_1,w_2,\alpha,t)$. Assume existence of equilibrium with equilibrum utilities $u^\star_1$ and $u_2^\star$ and select values for these ($u_2^\star > u_1^\star$ rich get higher utlity in equilibrium). Also select a value for the alternative use of land $r_a$ referred to as the land rent (not to be confused by the rent paid for land by producers of housing if that sector was endogenized in the model).
In the figure below the bid function $p_j(u_j,d)$ of the high and low income are illustrated. Residents will live where they bid the highest value for housing. The city expands to where the bid-rents of either income class is higher than the land rent.





Two observations are made:
Observation (1): The low income resident will live in the center of the city (this requires housing to be normal good)
Observation (2): There are perfect income stratification with low income in one area and high income in another area
To close the model the housing market must be in equilibrium. Let $[0,\hat d]$ be the locations where the low income reside. Then for $d \in [0,\hat d]$ housing demand is given as $h_{j}^m(u_{j},d)$ per person and units of housing is $H(d)$ so population density at $d$ must be $H(d)/h_{j}^m(u_{j},d)$ and integrating over $[0,\hat d]$ this must be equal to $N_1$ the number of low income residents I have exogenously stipulated to live in the city. So equilibrium conditions includes that
$$N_1 = \int_0^{\hat d} \frac{H(t)}{h_{1}^m(u_{1},t)} dt$$
and
$$N_2 = \int_{\hat d}^{\bar d} \frac{H(t)}{h_{2}^m(u_{2},t)} dt.$$
As is hopefully apparent from the figure above it is the case that: Given two utilities $(u_1,u_2)$ the bid functions determine $\hat d$ and $\bar d$ - where $\bar d$ is the value of $d$ where $p_2(u_2,d) = r_a$ - and given these values the integrals can be calculated to get $N_1(u_1,u_2)$ and $N_2(u_1,u_2)$. When the equilibrium is unique these equations
$$N_1 = N_1(u_1,u_2) \ \ \ N_2=N_2(u_1,u_2),$$
will be invertible such that equilibrium utilities can be found from the exogenously given $N_1$ and $N_2$.
Proof of existence and uniqueness of equilibrium for a model with $K$ income classes and standard restrictions on the utility functions (not assuming Cobb Douglas) are given in
John Hartwick (1976) Comparative Statics of a Residential Economy with
Several Classes Journal of Economic Theory 13, 396-413.
which also include comparative statics of the model.
Proof of existence and uniqueness of equilibrium for a model is also given in
Fujita (1989) Urban economic theory: Land use and city size
which also offers an algorithm for solving for the equilibrium.
See also
Puga and Duranton (2015) Urban Land Use
for discussion of main results of this and other extensions and empirical relevance of the models (but not detailed exposition).
Note on housing supply: Instead of assuming housing to be exogenously given I could have followed Brueckner (1989) and assumed land to exogenously given $L(d)$ and then housing be produced by constant returns to scale production under perfect competition using capital and land with the capital supply being infinitely price elastic and having price $r_K$ in which case housing supplied is a function of $d$ through $p(d)$. So housing supply is a function of the price $H(p)$ but at any distance $d$ the price residents are willing to give is set by the bid function which is a function of $d$ hence either way we get $H(d)$.
Here is some really fast written R-code implementing the solution algorithm in Fujita (1989):
library(rootSolve)

a <- 0.35
t <- 200
w <- c(6000,12000)
N <- c(500,500)
J <- 2
x <- rep(1,J+1)

u <- rep(170,J)


# Function to make bid-max lot size for agent of type j
# with income and utility (w_j,u_j). 
# Bid-max lot size is a function of distance to center d
# given type specifics (w_j,u_j).

make_g <- function(u_j,w_j,t,a)
    {
        g <- function(d)
            {
                nominator <- a*u_j^(1/a)
                denominator <- (w_j-t*d)^((1-a)/a)*.6*pi*d
                bid_lot <- nominator/denominator
                # return should be 1/bid_lot
                # but denominator can be 0 so do not calculate it
                # and use it as denominator instead 
                denominator <- ifelse((w_j-t*d)>0,denominator,0)
                return(denominator/nominator)
            }
        return(g)
    }


# Solve the integral of g with respect to r from x_L = x_j to x_U = x_j+1
# such that the integral equals N_j 

# Initiate with j=1 so x_l is set to 0 and u_j is some guess
j <- 1
x_L <- x[j]
u_j <- u[j]
    
    # Get exogenous income and population size for relevant segment
    # of population.
    w_j <- w[j]
    N_j <- N[j]

    g <- make_g(u_j,w_j,t,a)
    v <- function(x_U) {return(N_j-integrate(g,x_L,x_U)$value)} 
    x[j+1] <- uniroot(v,c(x_L+0.0001,500))$root

# Initiate j=2
j <- 2
    
    # Find u_2 given I_1,I_2 at x and u_1
    u[j] <- u[j-1]*(w[j] - t*x[j])/(w[j-1]-t*x[j])

x_L <- x[j]
u_j <- u[j]

    # Get exogenous income and population size for relevant segment
    # of population.
    w_j <- w[j]
    N_j <- N[j]

    g <- make_g(u_j,w_j,t,a)
    v <- function(x_U) {return(N_j-integrate(g,x_L,x_U)$value)} 
    x[j+1] <- uniroot(v,c(x_L+0.0001,500))$root



# Calculate price at x[j+1]
j <- 3  
p_edge <- ((w[j-1] - t*x[j])/u[j-1])^(1/a)
x
u
p_edge

p_max <- (w[1]/u[1])^(1/a)+10
plot(x[1]:(x[J+1]+1),type="n",ylim=c(0,p_max),ylab="housing price",xlab="distance to city center")
colors <- c("blue","red")
for (j in 1:J)
    {
        x_L <- x[j]
        x_U <- x[j+1]

        d <- seq(0,x[J+1],length.out=100)
        u_j <- u[j] 
        p <- ((w[j] - t*d)/u_j)^(1/a)
        points(d,p,type="l",col=colors[j])
    }

abline(v=x[2],col="grey")
abline(h=p_edge,col="black")
legend(12,20000,c("bid function (low income)","bid function (high income)","sorting threshold",
"land rent"),
col=c("blue","red","grey","black"),cex=0.7,bty="n",pch=24)
#savePlot("monocentric.jpg",type="jpg")


---

