Henningsen and Henningsen 2011

Feb 25, 2020 · 2 min read

I read this paper because I needed to estimate a Constant Elasticity of Substitution (CES) production function. Typically, when estimating Total Factor Productivity (TFP), economists prefer the Cobb-Douglas function because it’s easy to linearize and its parameters have useful real-world interpretations. Unfortunately, for structural models like Computable General Equilibrium (CGE), estimating with Cobb-Douglas has the drawback that the substitution parameter is always 1. No modern CGE model uses a substitution elasticity of 1.

Unfortunately, CES estimation is quite difficult. Linearization is possible using Taylor series and its extensions, but if you want to use it for economic models, many restrictions are needed. The alternative is non-linear estimation, which comes in various algorithmic flavors. Plus, I needed step-by-step instructions for implementing it computationally.

{% assign author = site.data.article[page.author] %} is the perfect paper for me. They developed an R package called micEconCES. This package can perform CES estimation and documents the various algorithms available. They also survey the history of CES estimation from when CES was first proposed by Arrow and colleagues through to recent developments. Comprehensive stuff.

A quick refresher on CES: in general, the CES function takes the form:

y=γ(i=1nδixiρ)υρ\label1 y = \gamma \left(\sum_{i=1}^n \delta_i x_i^{\rho}\right)^{-\frac{\upsilon}{\rho}} \label{1}

where x represents factors of production (capital, labor, materials, etc.). Ideally, we want different substitution parameters for each pair of factors. uses an approach pioneered by Sato (1967) employing nested CES. For example, with 3 factors of production, instead of putting all factors into model \ref{1} with a single substitution parameter ρ\rho, Sato structures it as follows:

y=γ[δ(δ1x1ρ1+(1δ1)x2ρ1)ρ/ρ1+(1δ)δ2x3ρ/ρ2]υρ\label2 y = \gamma \left[\: \delta \: (\delta_1 x_1^{-\rho_1} + (1-\delta_1) x_2^{-\rho_1})^{\rho/\rho_1} + (1-\delta) \delta_2 x_3^{-\rho/\rho_2}\right]^{-\frac{\upsilon}{\rho}} \label{2}

With this nested CES method, we can have multiple substitution parameters, including the substitution between x1x_1 and x2x_2 (given by ρ/ρ1{\rho/\rho_1}), and the parameter between the nested x1x_1 and x2x_2 and x3x_3. Depending on your needs as a researcher, the model can be configured. If you’re satisfied with a single ρ\rho for all factors, you can use model \ref{1} directly. Sometimes we feel that among 3 factors, two have a different substitution relationship with the third. For instance, if your factors are capital, skilled labor, and unskilled labor, you might use model \ref{2} where x1=skilled labor,x2=unskilled laborx_1=skilled\ labor, x_2=unskilled\ labor, and x3=capitalx_3 = capital.