Abstract

Explaining changes in productivity involves explaining changes in output and input quantities. Several economic models can be used for this purpose. This paper considers a model that accounts for weather and output price uncertainty. Changes in productivity are then explained in two steps. First, a stochastic production frontier model is used to decompose a proper productivity index into measures of technical progress, environmental change, technical efficiency change, scale-and-mix efficiency change, and changes in statistical noise. Second, a system of input demand equations is used to further decompose the measure of scale-and-mix efficiency change into a measure of technical progress, a measure of input price change, various measures of changes in expectations, and a measure of changes in allocative efficiency and statistical noise. The methodology is applied to U.S. agricultural data. The effects of weather and climate change on agricultural productivity are found to be small relative to the effects of changes in input prices.

1 Introduction

Most economists consider measures of productivity change to be measures of output quantity change divided by measures of input quantity change (e.g. Jorgenson and Griliches, 1967, p.250; Schreyer, 2001, p.11; O’Donnell, 2018, p.11). Measuring changes in output and input quantities (and therefore productivity) is an index number problem. One of the more worrying features of applied research on productivity is that most authors use output and input quantity indexes that do not satisfy a set of basic axioms from index theory. Indexes that do not satisfy at least one axiom include the well-known Fisher, Törnqvist, Malmquist, Elteto-Koves-Szulc (EKS), and Caves-Christensen-Diewert (CCD) indexes. The axioms that are violated most frequently are a transitivity axiom and a proportionality axiom. This has serious implications: indexes that violate these two axioms will say that outputs and inputs have increased and/or decreased when they may in fact have done the opposite, or perhaps not changed at all. This paper measures changes in output and input quantities using what O’Donnell (2016, 2018) calls proper indexes. Proper quantity indexes satisfy six basic axioms from index theory, including transitivity and proportionality.

If measures of productivity change are considered to be measures of output quantity change divided by measures of input quantity change, then explaining changes in productivity necessarily involves explaining changes in output and input quantities. Economists have many models that can be used for this purpose. For example, a widely-used model is one in which firms are assumed to be price takers in both output and input markets, and where managers are assumed to choose output and input quantities to maximize profits. Economists who use these types of models typically assume that market prices and characteristics of production environments are known at the time production decisions are made. These are unrealistic assumptions, at least in industries like agriculture. This paper considers an economic model that accounts for both weather and output price uncertainty. Estimating the model involves estimating a stochastic production frontier and a system of input demand equations. The paper shows how the estimated parameters of the model can be used to assess the effects of weather and climate change on total factor productivity (TFP).

The basic economic idea behind the paper is that changes in weather and climate potentially affect agricultural inputs and outputs (and therefore productivity) through two channels: first, realizations of weather variables potentially affect the outputs that can be produced using predetermined inputs; and, second, expectations about weather and climate potentially affect the input and planned output choices of managers.

Most previous studies that have examined the first channel have focused on measures of partial factor productivity (PFP): Anand and Khetarpal (2015), for example, use a biophysical simulation model to examine the effects of changes in surface air temperatures on wheat yield per hectare; and Nastis et al. (2012) use a production function model to estimate the effects of changes in temperature and precipitation on land productivity. And while several studies have looked at the effects of weather variables on TFP, most have used TFP indexes that have poor axiomatic properties: Salim and Islam (2010), for example, use a vector error correction (VEC) model to estimate the effects of changes in precipitation on Törnqvist index numbers; and Hughes et al. (2011) use a stochastic production frontier model to construct a climate effects index, which they then use to deflate Fisher index numbers. Only a handful of studies appear to have used proper TFP indexes: Sabasi and Shumway (2018), for example, use a seemingly unrelated regression (SUR) model to estimate the effects of changes in temperature and precipitation on Lowe index numbers; and Njuki et al. (2018, 2020) use stochastic production frontier models to estimate the effects of changes in temperature and precipitation on multiplicative index numbers. This paper goes a step further by considering the second channel: it considers the way in which expectations about weather and climate potentially affect input and planned output choices. No previous studies appear to have examined this second channel.

The structure of the paper is as follows. Section 2 deals with the problem of measuring TFP change.1 It starts by using artificial data to demonstrate the properties of various indexes. It then uses U.S. agricultural data to demonstrate that there are real-world situations where the choice of index matters. Section 3 considers the first channel by which weather and climate potentially affect productivity. It starts by writing the relationship between observed outputs, observed inputs and observed weather variables in the form of a stochastic production frontier model. It then uses the estimated parameters of the model to decompose a proper TFP index into measures of technical progress, environmental change, technical efficiency change, scale-and-mix efficiency change, and changes in statistical noise. Section 4 considers the second channel. It starts by writing the relationship between observed input quantities, observed input prices, expected output prices, and expected weather variables in the form of a system of input demand equations. It then uses the estimated parameters of the system to further decompose the measure of scale-and-mix efficiency change into a measure of technical progress, a measure of input price change, various measures of changes in expectations, and a measure of changes in allocative efficiency and statistical noise. Section 4 summarizes the paper and offers some concluding remarks.

2 Measuring TFP change

Measures of TFP change are measures of total output quantity change (i.e. output quantity indexes) divided by measures of total input quantity change (i.e. input quantity indexes). Computing output and input quantity index numbers is a matter of assigning numbers to baskets of outputs and inputs. Measurement theory says such numbers cannot be assigned in an arbitrary way. Instead, they must be assigned in such a way that the relationships between the numbers reflect the relationships between the baskets (Tal, 2016). This section uses artificial data on apples and oranges to illustrate this basic principle. It then uses U.S. agricultural data to demonstrate that there are situations where the choice of index matters.

2.1 Comparing apples and oranges

Consider the baskets of apples and oranges presented in Table 1. This table also presents quantity index numbers that have been computed using six different indexes: the MEW index is a type of multiplicative index that gives the products equal weight; the MOLS index is a multiplicative index that uses OLS parameter estimates as weights; the geometric Young (GY) index is a multiplicative index that uses average value shares as weights; the binary Törnqvist (BT) index is a ‘superlative’ index that uses observation-varying value shares as weights; the chained Törnqvist (CT) index is obtained by treating the observations as time-series observations and chaining pairs of BT indexes; and the CCD index is obtained by treating the observations as cross-section observations and taking geometric averages of sets of BT indexes. Technical details concerning each of these indexes can be found in O’Donnell (2018, Ch.3).2R code for computing the numbers in Table 1 can be found in Figure A1 in the Appendix.

Table 1.

Quantity index numbers.

MEWMOLSGYBTCTCCD
A =111111
graphic
B =1.7321.2091.8601.581*1.581*1.688*
graphic
C =2.4492.7972.3862.5983.150*2.700
graphic
D =1.7321.2091.8602.215*2.508*2.118*
graphic
E =22221.886*1.884*
graphic
MEWMOLSGYBTCTCCD
A =111111
graphic
B =1.7321.2091.8601.581*1.581*1.688*
graphic
C =2.4492.7972.3862.5983.150*2.700
graphic
D =1.7321.2091.8602.215*2.508*2.118*
graphic
E =22221.886*1.884*
graphic

Notes: MEW = Multiplicative with Equal Weights; MOLS = Multiplicative with OLS weights; GY = Geometric Young; BT = binary Törnqvist; CT = chained Tornqvist; CCD = Caves-Christensen-Diewert.

*Incoherent (and not because of rounding errors).

Table 1.

Quantity index numbers.

MEWMOLSGYBTCTCCD
A =111111
graphic
B =1.7321.2091.8601.581*1.581*1.688*
graphic
C =2.4492.7972.3862.5983.150*2.700
graphic
D =1.7321.2091.8602.215*2.508*2.118*
graphic
E =22221.886*1.884*
graphic
MEWMOLSGYBTCTCCD
A =111111
graphic
B =1.7321.2091.8601.581*1.581*1.688*
graphic
C =2.4492.7972.3862.5983.150*2.700
graphic
D =1.7321.2091.8602.215*2.508*2.118*
graphic
E =22221.886*1.884*
graphic

Notes: MEW = Multiplicative with Equal Weights; MOLS = Multiplicative with OLS weights; GY = Geometric Young; BT = binary Törnqvist; CT = chained Tornqvist; CCD = Caves-Christensen-Diewert.

*Incoherent (and not because of rounding errors).

The MEW, MOLS, and GY indexes are proper indexes. Among other things, this means they yield numbers that are consistent with measurement theory. Observe, for example, that the MEW index number in row B of Table 1 is the same as the number in row D, reflecting the fact that basket B contains the same number of apples and oranges as basket D; and the GY index number in row E is twice as big as the number in row A, reflecting the fact that basket E contains twice as many apples and oranges as basket A.

The BT, CT, and CCD indexes are not proper indexes. Consequently, they do not yield numbers that are consistent with measurement theory. Observe, for example, that the BT index number in row B differs from the number in row D, even though baskets B and D contain the same numbers of apples and oranges; the CT index number in row C is 3.15 times greater than the number in row A, even though basket C contains less than 3 times as many apples and oranges as basket A; and the CCD index number in row E is less than twice the number in row A, even though basket E contains exactly twice as many apples and oranges as basket A. The basic problem with BT, CT, and CCD quantity indexes is that they measure changes in value shares as well as quantities. Of course, if there is no variation in value shares, then they will only measure changes in quantities, as required; in that case, BT, CT, and CCD index numbers will be equal to GY index numbers.

2.2 TFP change in U.S. agriculture

This subsection uses different indexes to make time-series and cross-section comparisons of output, input, and TFP change in U.S. agriculture. The dataset comprises observations on the prices and quantities of four inputs (capital, land, labor, and materials) and three outputs (livestock, crops, and other outputs) in forty-eight states over the years from 1961 to 2004. The data were assembled by the Economic Research Service (ERS) of the U.S. Department of Agriculture (USDA). Details concerning the data can be found in Ball et al. (1997).

Several sets of index numbers are summarized in Fig. 1. The panels on the left-hand side of this figure present GY, BT, and CT indexes of output, input, and TFP change in Alabama from 1961 to 2004. The panels on the right-hand side present GY, BT, and CCD indexes of output, input, and TFP change across six states in 1961. Observe from the panels on the left-hand side that different indexes sometimes tell quite different stories about output, input, and/or TFP change over time: in panel (c), for example, the GY index tells us that input use in Alabama in 1983 was 9% lower in than it had been in 1961, while the CT index tells us that it was 5% higher. The panels on the right-hand side reveal that different indexes can also tell quite different stories about output, input, and/or TFP change across space: in panel (f), for example, the GY index tells us that Oklahoma was 3% less productive than Alabama in 1961, while the BT and CCD indexes tell us that Oklahoma was at least 11% more productive.

Output, Input, and TFP change in selected years and states (Alabama in 1961 = 1).
Figure 1.

Output, Input, and TFP change in selected years and states (Alabama in 1961 = 1).

2.3 The choice of index

There are at least three reasons why economists should choose proper indexes. First, all proper indexes satisfy important axioms from index theory (e.g. transitivity and proportionality). Second, and relatedly, proper quantity indexes always yield numbers that are consistent with measurement theory (i.e. the patterns in the numbers always mirror the patterns in the quantities). Third, any given proper index can be used to make comparisons across both time and space; in contrast, the CT index, for example, can only be used to make (erroneous) comparisons across time, while the CCD index can only be used to make (erroneous) comparisons across space.

In practice, the choice of proper index is largely a matter of taste. In the case of output indexes, for example, O’Donnell (2018) suggests that ‘analysts who want to minimize the amount of variation in the index numbers’ should use a multiplicative index with weights that are computed within a regression framework (p.97); the GY index ‘should be used by analysts who regard revenue shares as [appropriate] measures of relative value (e.g. analysts who might otherwise use a Tornqvist, chained Tornqvist or CCD index)’ (p.96); a primal index should be used by ‘analysts who regard marginal rates of transformation as [appropriate] measures of relative value (e.g. analysts who might otherwise use a generalized Malmquist index)’ (p.97); and a benefit-of-the-doubt (BOD) ‘index should be used by analysts who believe measures of relative value should vary from one output comparison to the next. It can also be used by analysts who have no information about output prices, revenue shares, or production technologies’ (p.99).

The aim of this paper is to estimate the effects of weather and climate change on U.S. agricultural productivity. The first step is to measure changes in outputs, inputs, and TFP. GY indexes will be used for this purpose, mainly because they are proper indexes that combine nicely with the double-log econometric models used in Sections 3 and 4. For a precise definition of these indexes, let qit = (q1it, …, qNit) and xit = (x1it, …, xMit) denote the output and input vectors of state i in period t. The GY indexes that compare the outputs and inputs of state i in period t with the outputs and inputs of state k in period s are
(1)
(2)
where |$\bar{r}_1, \dots , \bar{r}_N$| are average revenue shares and |$\bar{s}_1, \dots , \bar{s}_M$| are average cost shares. The associated GY index that compares the TFP of state i in period t with the TFP of state k in period s is
(3)
These indexes can be traced back at least as far as O’Donnell (2016). They were used to produce the GY index numbers presented in Fig. 1. GY TFP index numbers for all states in all periods are presented in Table A3 in the Appendix. The remaining sections of this paper will use a combination of economic theory and statistical methods to explain variations in these numbers.

3 Explaining changes in TFP

If TFP change is defined as a measure of total output quantity change divided by a measure of total input quantity change, then explaining changes in TFP necessarily involves explaining changes in outputs and inputs. This paper assumes that managers choose outputs and inputs in two stages: first, at the beginning of the production period, managers choose inputs and planned outputs to maximize expected profits in the face of uncertainty about output prices and one or more characteristics of the production environment (e.g. rainfall); and, second, after inputs have been chosen, managers seek to maximize the outputs that can be produced using their chosen inputs in whatever production environment is realized. This section focuses on the second stage. If this second-stage behavior is true, then the relationship between outputs, inputs, and environmental variables can be written in the form of a stochastic production frontier model in which the explanatory variables are exogenous. This section uses this type of model to explain variations in the GY TFP index numbers reported in Section 2.2. To do this, the USDA input and output data used in that section is supplemented with observations on three characteristics of the production environment: DD830 measures the number of degree days between 8C and 30C between March and August; DD30 measures the number of degree days above 30C between March and August; and PREC measures total precipitation in inches between March and August. These weather variables were the only environmental variables that were available; data on these variables were also assembled and supplied by the USDA.

3.1 The stochastic production frontier model

Let Pt(xit, zit) denote the set of outputs that can be produced using the input vector xit in period t in a production environment characterized by the vector zit = (z1it, …, zJit). Under weak regularity conditions, this set can be represented by the following period-and-environment-specific output distance function:
(4)
This function is linearly homogeneous in outputs; this implies that |$D^{t}_O(x_{it},q_{it},z_{it})=q_{1it}D^{t}_O(x_{it},q^{*}_{it},z_{it})$| where |$q^{*}_{it}\equiv q_{it}/q_{1it}$| denotes a vector of normalized outputs. Equivalently, after some simple algebra, the relationship between outputs, inputs, and environmental variables can be written as
(5)
where |$u_{it}\equiv -\ln D_O^t(x_{it},q_{it},z_{it})\ge 0$| is an output-oriented technical inefficiency effect.3 Unfortunately, the functional form of the output distance function is unknown. In practice, we must replace the negative of the logarithm of the output distance function with an arbitrary approximating function. The relationship between outputs, inputs, and environmental variables then becomes
(6)
where f t(.) is the approximating function and |$v_{it}=-\ln D^{t}_O(x_{it},q^{*}_{it},z_{it})-f^t(x_{it},q^{*}_{it},z_{it})$| is a variable that accounts for potential functional form errors and other sources of statistical noise. An approximating function that provides for a neat decomposition of GY TFP index numbers is
(7)
where αi is a fixed effect that accounts for nonstochastic time-invariant characteristics of the production environment (e.g. soil type and elevation), dhit is a trend variable that takes the value t in decade h and 0 beforehand (and 10 afterwards), zjit is the jth element of zit, and |$Q(q_{it})=\prod _{n=1}^{N}q^{\bar{r}_n}_{nit}$| is a GY measure of total output. This is the approximating function used in this paper. Equations (6) and (7) together imply that
(8)
This is a stochastic production frontier model with the same basic structure as the model of Aigner et al. (1977, Eq.7).4 Unfortunately, the slope parameters can only be given an economic interpretation if we make some assumptions about the nature of statistical noise: if we assume that vit does not depend on dhit, for example, then λh can be interpreted as a measure of technical progress in decade h; if we assume that vit does not depend on zjit, then δj can be interpreted as an elasticity that measures the percent change in output due to a one percent increase in the jth environmental variable; and if we assume that vit does not depend on xmit, then βm can be interpreted as an elasticity that measures the percent increase in output due to a one percent increase in the mth input. This paper makes these assumptions and interprets the parameters accordingly.

3.2 Estimation

The unknown parameters in stochastic frontier models are most often estimated using the method of maximum likelihood (ML). The are two problems with the ML approach. First, there is no satisfactory way of incorporating inequality information into the estimation process; this is because binding inequality constraints will lead to parameter estimates that lie on the constraint boundary, with standard errors of zero (implying we know the values of the parameters with certainty). Second, we cannot make valid finite-sample inferences concerning nonlinear functions of the parameters; this because all the theory underpinning ML estimation is asymptotic. For these reasons, this paper estimates the parameters in (8) using a Bayesian approach. Bayesian estimation involves summarizing the information we have about the unknown parameters in the form of a joint posterior probability density function (pdf). One source of information is the data; this so-called sample information is summarized in the form of the usual likelihood function (i.e. the function that is maximized in the ML approach). Other sources of information include economic theory and common sense; this so-called nonsample information is summarized in the form of a prior pdf. The likelihood function and the prior pdf are combined using Bayes’s theorem to form the joint posterior pdf.

In practice, the exact form of the likelihood function depends on the assumed distributions of the noise and inefficiency effects. In this paper, the likelihood function is formed under the assumption that vit is an independent half-normal random variable with precision parameter h, and that uit is an independent exponential random variable with scale parameter λ. These are the two most common assumptions made in the Bayesian stochastic frontier literature. The prior pdf is a proper pdf of the form
(9)
where θ is a vector containing all the unknown parameters in (8), fN(θ|0, B) denotes a multivariate normal pdf with mean vector zero and covariance matrix B, fG(h|1, c) denotes a gamma pdf with shape parameter one and rate parameter c, I(.) is an indicator function that takes the value one if the argument is true (and 0 otherwise), and R denotes the admissible region of θ. In this empirical application, R simply says that δ2 must be nonpositive and all other slope coefficients except δ3 must be nonnegative. The prior hyperparameters are set to c = 0.0001, τ = 0.9 and B = c−1I where I is an identity matrix. These values mean that, apart from the inequality constraints, the prior contains almost no information about the unknown parameters.

Evaluating characteristics of marginal posterior pdfs for the individual parameters (e.g. means, variances) involves evaluating multiple integrals. In practice, this involves sampling from the joint posterior pdf. In this paper, samples of size one million were drawn using the Markov Chain Monte Carlo (MCMC) sampling package of Plummer (2019). The thinning interval for monitors was set to 100. The convergence of the chains was confirmed both visually and using the convergence diagnostic of Geweke (1992). Relevant R code is presented in Figure A2 in the Appendix.

Bayesian point estimates of the first five intercept parameters and all the slope parameters are presented in Table 2. Estimates of all forty-eight intercept parameters are presented in Table A1 in the Appendix. These tables also report associated estimated standard errors and lower and upper bounds of 95% highest posterior density (HPD) intervals. HPD intervals are the Bayesian counterparts to confidence intervals. By construction, the point and interval estimates in Table 2 are all consistent with prior expectations. Given the assumptions made earlier about the nature of statistical noise, they have straightforward interpretations: the point estimates of α1 and α4, for example, indicate that fixed characteristics of the production environment in California are |$46\%$| more favorable for agricultural production than fixed characteristics of the production environment in Alabama; the estimates of λ1, …, λ5 indicate there was technical progress in every decade, and that the rate of technical progress reached a maximum of 2.3% per annum in the 1980s; the estimates of β1, …, β5 indicate that, all other things being equal, increasing any input will lead to an increase in output, and that the production frontier exhibits decreasing returns to scale (the elasticity of scale is estimated to be 0.858); the estimate of δ1 indicates that an increase in the number of so-called ‘good degree days’ will increase output; the estimate of δ2 indicates that an increase in the number of ‘bad degree days’ will decrease output; and the estimate of δ3 indicates that an increase in precipitation will increase output.

Table 2.

Estimated parameters of the stochastic frontier model.

Est.St. Err.2.5%97.5%
α1AL1.6670.3091.0632.189
α2AR1.9630.3131.3482.495
α3AZ1.9160.3031.3012.428
α4CA2.4260.3291.7762.991
α5CO1.7260.3111.1132.251
λ1t in the 1960s0.0060.0010.0030.008
λ2t in the 1970s0.0030.0010.0020.005
λ3t in the 1980s0.0230.0010.0210.025
λ4t in the 1990s0.0070.0010.0050.009
λ5t in the 2000s0.0090.0030.0040.014
β1Capital0.1550.0200.1190.195
β2Land0.0180.0120.0010.045
β3Labor0.1050.0090.0860.123
β4Materials0.5810.0130.5560.606
δ1DD8300.0190.0160.0010.058
δ2DD30−0.0170.002−0.021−0.013
δ3PREC0.0060.009−0.0110.022
Est.St. Err.2.5%97.5%
α1AL1.6670.3091.0632.189
α2AR1.9630.3131.3482.495
α3AZ1.9160.3031.3012.428
α4CA2.4260.3291.7762.991
α5CO1.7260.3111.1132.251
λ1t in the 1960s0.0060.0010.0030.008
λ2t in the 1970s0.0030.0010.0020.005
λ3t in the 1980s0.0230.0010.0210.025
λ4t in the 1990s0.0070.0010.0050.009
λ5t in the 2000s0.0090.0030.0040.014
β1Capital0.1550.0200.1190.195
β2Land0.0180.0120.0010.045
β3Labor0.1050.0090.0860.123
β4Materials0.5810.0130.5560.606
δ1DD8300.0190.0160.0010.058
δ2DD30−0.0170.002−0.021−0.013
δ3PREC0.0060.009−0.0110.022
Table 2.

Estimated parameters of the stochastic frontier model.

Est.St. Err.2.5%97.5%
α1AL1.6670.3091.0632.189
α2AR1.9630.3131.3482.495
α3AZ1.9160.3031.3012.428
α4CA2.4260.3291.7762.991
α5CO1.7260.3111.1132.251
λ1t in the 1960s0.0060.0010.0030.008
λ2t in the 1970s0.0030.0010.0020.005
λ3t in the 1980s0.0230.0010.0210.025
λ4t in the 1990s0.0070.0010.0050.009
λ5t in the 2000s0.0090.0030.0040.014
β1Capital0.1550.0200.1190.195
β2Land0.0180.0120.0010.045
β3Labor0.1050.0090.0860.123
β4Materials0.5810.0130.5560.606
δ1DD8300.0190.0160.0010.058
δ2DD30−0.0170.002−0.021−0.013
δ3PREC0.0060.009−0.0110.022
Est.St. Err.2.5%97.5%
α1AL1.6670.3091.0632.189
α2AR1.9630.3131.3482.495
α3AZ1.9160.3031.3012.428
α4CA2.4260.3291.7762.991
α5CO1.7260.3111.1132.251
λ1t in the 1960s0.0060.0010.0030.008
λ2t in the 1970s0.0030.0010.0020.005
λ3t in the 1980s0.0230.0010.0210.025
λ4t in the 1990s0.0070.0010.0050.009
λ5t in the 2000s0.0090.0030.0040.014
β1Capital0.1550.0200.1190.195
β2Land0.0180.0120.0010.045
β3Labor0.1050.0090.0860.123
β4Materials0.5810.0130.5560.606
δ1DD8300.0190.0160.0010.058
δ2DD30−0.0170.002−0.021−0.013
δ3PREC0.0060.009−0.0110.022

Associated with the parameter estimates reported in Table 2 are estimates of output-oriented technical efficiency (OTE). Point estimates for the first five states in both 1961 and 2004 are presented in Table 3. Estimates for all states in all years are presented in Table A2 in the Appendix. These tables also report associated estimated standard errors and lower and upper bounds of 95% HPD intervals. Among other things, these estimates indicate that Alabama farmers were slightly less technically efficient in 2004 than they had been in 1961. The average level of technical efficiency across all states in all years is estimated to be 0.959.

Table 3.

Selected estimates of OTE.

StateYearEst.St. Err.2.5%97.5%
AL19610.9790.0190.9300.999
AR19610.9130.0470.8150.992
AZ19610.9710.0240.9110.999
CA19610.9210.0460.8240.994
CO19610.9650.0280.8950.999
AL20040.9570.0330.8790.998
AR20040.9840.0150.9441.000
AZ20040.9790.0190.9300.999
CA20040.9780.0200.9260.999
CO20040.9690.0260.9060.999
StateYearEst.St. Err.2.5%97.5%
AL19610.9790.0190.9300.999
AR19610.9130.0470.8150.992
AZ19610.9710.0240.9110.999
CA19610.9210.0460.8240.994
CO19610.9650.0280.8950.999
AL20040.9570.0330.8790.998
AR20040.9840.0150.9441.000
AZ20040.9790.0190.9300.999
CA20040.9780.0200.9260.999
CO20040.9690.0260.9060.999
Table 3.

Selected estimates of OTE.

StateYearEst.St. Err.2.5%97.5%
AL19610.9790.0190.9300.999
AR19610.9130.0470.8150.992
AZ19610.9710.0240.9110.999
CA19610.9210.0460.8240.994
CO19610.9650.0280.8950.999
AL20040.9570.0330.8790.998
AR20040.9840.0150.9441.000
AZ20040.9790.0190.9300.999
CA20040.9780.0200.9260.999
CO20040.9690.0260.9060.999
StateYearEst.St. Err.2.5%97.5%
AL19610.9790.0190.9300.999
AR19610.9130.0470.8150.992
AZ19610.9710.0240.9110.999
CA19610.9210.0460.8240.994
CO19610.9650.0280.8950.999
AL20040.9570.0330.8790.998
AR20040.9840.0150.9441.000
AZ20040.9790.0190.9300.999
CA20040.9780.0200.9260.999
CO20040.9690.0260.9060.999

3.3 Decomposing GY TFP index numbers

O’Donnell (2018, Sect. 8.5.2) explains how the parameters of stochastic frontier models can be used to decompose proper TFP index numbers into economically meaningful components. This section uses the results reported in Section 3.2 to decompose the GY TFP index numbers reported earlier in Section 2.2. The easiest way forward is to take the antilogarithms of both sides of (8) and use the definition of the dependent variable in that equation to rewrite it as
(10)
A similar equation holds for state k in period s. These two equations can be substituted into Equation (3) to yield the following decomposition:
(11)
Our earlier assumptions about the nature of statistical noise mean that the first term in square brackets on the right-hand side can be viewed as an output-oriented technology index (OTI) (i.e. a measure of technical progress), the second term can be viewed as an output-oriented environment index (OEI) (i.e. a measure of environmental change), and the third term can be viewed as an output-oriented scale-and-mix efficiency index (OSMEI) (i.e. a measure of change in economies of scale and substitution). The fourth term is an output-oriented technical efficiency index (OTEI) (i.e. a measure of technical efficiency change), and the last term is a statistical noise index (SNI) (i.e. a measure of changes in omitted variables and other sources of statistical noise).

Figure 2 presents a decomposition of TFP change in selected states. The TFPI numbers depicted in the top-left-hand panel are the GY numbers presented earlier in panel (e) of Fig. 1. The OTI, OEI, and OSMEI numbers were computed by using the estimated slope parameters in Table 2 and the estimated intercept parameters in Table A1 in the Appendix to evaluate the first three terms on the right-hand side of equation (11). The OTEI numbers were computed using the point estimates of OTE reported in Table 3 and in Table A2 in the Appendix. The SNI numbers were computed as residuals. Results for all forty-eight states are presented in Table A3 in the Appendix.

TFP change in selected states from 1961 to 2004 (AL in 1961 = 1).
Figure 2.

TFP change in selected states from 1961 to 2004 (AL in 1961 = 1).

The results reported in Fig. 2 and Table A3 indicate that the main drivers of TFP change over time were technical progress and scale-and-mix efficiency change. In the case of Alabama, for example, TFP in 2004 was 1.776 times higher than it had been in 1961. This increase can be broken down as follows: TFPI = OTI × OEI × OSMEI × OTEI × SNI = 1.517 × 0.996 × 1.328 × 0.977 × 0.906 = 1.776. This decomposition indicates that, in Alabama, (i) technical progress provided for a 51.7% increase in TFP; (ii) changes in the production environment had virtually no impact on TFP; (iii) improvements in scale-and-mix efficiency led to a 32.8% increase in TFP; (iv) lower technical efficiency led to a 2.7% fall in TFP; and (v) changes in omitted variables and other sources of statistical noise accounted for a 9.4% fall in TFP.

The results reported in Fig. 2 and Table A3 also indicate that the main drivers of TFP change across states were differences in production environments and levels of scale-and-mix efficiency. In 1961, for example, Texas was 6.7% less productive than Alabama. This difference can be broken down as follows: TFPI = OTI × OEI × OSMEI × OTEI × SNI = 1 × 1.230 × 0.780 × 0.997 × 0.976 = 0.933. The reason the OTI component takes the value one is that Equation (8) plausibly does not allow for differences in rates of technical progress across states. The remaining components indicate that, in 1961, (i) a better production environment allowed Texas farmers to be 23% more productive than farmers in Alabama; (iii) Texas farmers were 22% less scale-and-mix efficient than farmers in Alabama; (iv) Texas farmers were only 0.3% less technically efficient than farmers in Alabama; and (v) omitted variables and other sources of statistical noise accounted for only 2.4% of the difference in TFP.

Finally, it is possible to decompose the environmental change component of TFP change (i.e. the OEI) into the product of an output-oriented fixed effects index (OFEI) and an output-oriented temperature and precipitation index (OTPI). The OFEI measures the direct effects of time-invariant environmental variables on TFP; mathematically, the index that compares the effects in state i with the effects in state k is exp (αi − αk). The OTPI measures the direct effects of changes in temperature and precipitation on TFP; mathematically, the index that compares the effects in state i in period t with the effects in state k in period s is |$\prod _{j=1}^{J}(z_{jit}/z_{jks})^{\delta _j}$|⁠. To illustrate, let us reconsider the finding that, in 1961, a better production environment allowed Texas farmers to be 23% more productive than farmers in Alabama. This difference can be broken down as follows: OEI = OFEI × OTPI = 1.265 × 0.973 = 1.23. This indicates that time-invariant characteristics of the production environment (e.g. soil type and terrain) allowed Texas farmers to be 26.5% more productive than farmers in Alabama, and that, in 1961, differences in temperature and precipitation accounted for only 2.7% of the difference in TFP. Results for other states and other years also indicate that changes in weather had a relatively small effect on TFP.

4 Explaining changes in OSME

Explaining changes in OSME involves explaining changes in input quantities. This section uses a system of input demand equations to decompose the OSME index numbers reported in Section 3.3. To do this, the USDA data used in previous sections were supplemented with observations on expected aggregate output prices and weather variables. Expected aggregate output prices were measured using one-period lagged values of GY output price index numbers, and expected values of weather variables were measured using arithmetic averages of observed weather variables over the previous 10 years. To the extent that they measure average atmospheric conditions over a long period of time, changes in these expected weather variables can be viewed as measures of climate change.

4.1 The input demand functions

In Section 3 it was assumed that inputs are chosen to maximize expected profits in the face of uncertainty about output prices and one or more characteristics of the production environment. If this assumption is true, then the mth input demand function can be written in the form
(12)
where αmi is an unobserved fixed effect that again accounts for nonstochastic time-invariant characteristics of the production environment (e.g. terrain), |$p^e_{it}$| is the expected aggregate output price, |$z^e_{jit}$| is the expected value of the jth environmental variable, wkit is the kth input price, and emit represents a combination of allocative inefficiency (i.e. failure to choose input combinations that maximize expected profits) and statistical noise. Again, the slope parameters can only be given a meaningful interpretation if we make some assumptions about the error term: if we assume that emit does not depend on t, for example, then λm can be interpreted as an input-specific measure of technical progress; if we assume that emit does not depend on |$p^e_{it}$|⁠, then ϕm can be interpreted as an elasticity that measures the percent change in the demand for input m due to a one percent increase in the expected aggregate output price; if we assume that emit does not depend on |$z^e_{jit}$|⁠, then δmj can be interpreted as an elasticity that measures the percent change in the demand for input m due to a one percent increase in |$z^e_{jit}$|⁠; and if we assume that emit does not depend on whit, then ξmh can be interpreted as an elasticity that measures the percent change in the demand for input m due to a one percent increase in the price of input h. Again, this paper makes these assumptions and interprets the parameters accordingly.

4.2 Estimation

In this paper, each of the input demand equations defined by (12) is estimated separately using a Bayesian approach. The likelihood function for the mth input demand equation is formed under the assumption that emit is distributed as an independent |$N(0,\sigma _m^2)$| random variable. The prior pdf is chosen to be a truncated multivariate normal distribution that once again contains almost no information about the parameters: all it says is that ϕm must be nonnegative, and that λm and ξmm must be nonpositive. Again, characteristics of marginal posterior pdfs are evaluated using MCMC sampling. Again, samples of size one million were drawn using the sampling package of Plummer (2019). The thinning interval for monitors was again set to 100. The convergence of the chains was again confirmed both visually and using the convergence diagnostic of Geweke (1992). R code for estimating the parameters of the input demand equation for capital is presented in Figure A3 in the Appendix.

Bayesian point estimates of the slope parameters in each input demand equation are presented in Table 4. Estimates of the intercept parameters are presented in Tables A4 to A7 in the Appendix. These tables also report associated standard errors and lower and upper bounds of 95% HPD intervals. All of the estimates reported in these tables are plausible: the estimates of λ1, …, λ4, for example, indicate that technical progress was capital-, land-, and labor-saving; the estimates of ϕ1, …, ϕ4 indicate that increases in expected output prices lead to increases in the demand for all inputs; and the estimates of ξ11, ξ22, ξ33, and ξ44 indicate that, all other things being equal, an increase in the price of any input leads to a fall in the demand for that input.

Table 4.

Estimated slope parameters of the input demand equations.

Est.St. Err.2.5%97.5%
Demand for capital:
 λ1t−0.0050.001−0.007−0.004
 ϕ1E(output price)0.0170.0150.0000.056
 ξ11Price of capital−0.2520.026−0.302−0.202
 ξ12Price of land0.1450.0080.1310.160
 ξ13Price of labor−0.1650.011−0.187−0.143
 ξ14Price of materials0.2540.0250.2030.302
 δ11E(DD830)−1.4650.117−1.644−1.221
 δ12E(DD30)−0.1230.015−0.153−0.093
 δ13E(PREC)−0.1800.051−0.278−0.085
Demand for land:
 λ2t−0.0050.000−0.005−0.004
 ϕ2E(output price)0.0030.0030.0000.012
 ξ21Price of capital0.1710.0170.1380.204
 ξ22Price of land−0.1200.005−0.130−0.110
 ξ23Price of labor0.0030.007−0.0120.017
 ξ24Price of materials−0.0560.016−0.087−0.026
 δ21E(DD830)−0.0630.110−0.2560.120
 δ22E(DD30)−0.0560.011−0.077−0.035
 δ23E(PREC)0.2850.0340.2180.350
Demand for labor:
 λ3t−0.0080.001−0.009−0.006
 ϕ3E(output price)0.0130.0130.0000.046
 ξ31Price of capital0.1350.0330.0700.200
 ξ32Price of land−0.1550.010−0.174−0.137
 ξ33Price of labor−0.2520.015−0.281−0.223
 ξ34Price of materials0.2590.0320.1970.320
 δ31E(DD830)0.1750.115−0.0180.392
 δ32E(DD30)−0.0110.019−0.0480.028
 δ33E(PREC)0.1480.0680.0170.277
Demand for materials:
 λ4t0.0000.0000.0000.000
 ϕ4E(output price)0.0040.0040.0000.016
 ξ41Price of capital0.1490.0290.0920.206
 ξ42Price of land0.0250.0080.0090.040
 ξ43Price of labor0.1190.0120.0960.142
 ξ44Price of materials−0.2970.021−0.338−0.256
 δ41E(DD830)−0.2410.324−1.0700.099
 δ42E(DD30)−0.1410.023−0.181−0.087
 δ43E(PREC)0.0890.067−0.0430.225
Est.St. Err.2.5%97.5%
Demand for capital:
 λ1t−0.0050.001−0.007−0.004
 ϕ1E(output price)0.0170.0150.0000.056
 ξ11Price of capital−0.2520.026−0.302−0.202
 ξ12Price of land0.1450.0080.1310.160
 ξ13Price of labor−0.1650.011−0.187−0.143
 ξ14Price of materials0.2540.0250.2030.302
 δ11E(DD830)−1.4650.117−1.644−1.221
 δ12E(DD30)−0.1230.015−0.153−0.093
 δ13E(PREC)−0.1800.051−0.278−0.085
Demand for land:
 λ2t−0.0050.000−0.005−0.004
 ϕ2E(output price)0.0030.0030.0000.012
 ξ21Price of capital0.1710.0170.1380.204
 ξ22Price of land−0.1200.005−0.130−0.110
 ξ23Price of labor0.0030.007−0.0120.017
 ξ24Price of materials−0.0560.016−0.087−0.026
 δ21E(DD830)−0.0630.110−0.2560.120
 δ22E(DD30)−0.0560.011−0.077−0.035
 δ23E(PREC)0.2850.0340.2180.350
Demand for labor:
 λ3t−0.0080.001−0.009−0.006
 ϕ3E(output price)0.0130.0130.0000.046
 ξ31Price of capital0.1350.0330.0700.200
 ξ32Price of land−0.1550.010−0.174−0.137
 ξ33Price of labor−0.2520.015−0.281−0.223
 ξ34Price of materials0.2590.0320.1970.320
 δ31E(DD830)0.1750.115−0.0180.392
 δ32E(DD30)−0.0110.019−0.0480.028
 δ33E(PREC)0.1480.0680.0170.277
Demand for materials:
 λ4t0.0000.0000.0000.000
 ϕ4E(output price)0.0040.0040.0000.016
 ξ41Price of capital0.1490.0290.0920.206
 ξ42Price of land0.0250.0080.0090.040
 ξ43Price of labor0.1190.0120.0960.142
 ξ44Price of materials−0.2970.021−0.338−0.256
 δ41E(DD830)−0.2410.324−1.0700.099
 δ42E(DD30)−0.1410.023−0.181−0.087
 δ43E(PREC)0.0890.067−0.0430.225
Table 4.

Estimated slope parameters of the input demand equations.

Est.St. Err.2.5%97.5%
Demand for capital:
 λ1t−0.0050.001−0.007−0.004
 ϕ1E(output price)0.0170.0150.0000.056
 ξ11Price of capital−0.2520.026−0.302−0.202
 ξ12Price of land0.1450.0080.1310.160
 ξ13Price of labor−0.1650.011−0.187−0.143
 ξ14Price of materials0.2540.0250.2030.302
 δ11E(DD830)−1.4650.117−1.644−1.221
 δ12E(DD30)−0.1230.015−0.153−0.093
 δ13E(PREC)−0.1800.051−0.278−0.085
Demand for land:
 λ2t−0.0050.000−0.005−0.004
 ϕ2E(output price)0.0030.0030.0000.012
 ξ21Price of capital0.1710.0170.1380.204
 ξ22Price of land−0.1200.005−0.130−0.110
 ξ23Price of labor0.0030.007−0.0120.017
 ξ24Price of materials−0.0560.016−0.087−0.026
 δ21E(DD830)−0.0630.110−0.2560.120
 δ22E(DD30)−0.0560.011−0.077−0.035
 δ23E(PREC)0.2850.0340.2180.350
Demand for labor:
 λ3t−0.0080.001−0.009−0.006
 ϕ3E(output price)0.0130.0130.0000.046
 ξ31Price of capital0.1350.0330.0700.200
 ξ32Price of land−0.1550.010−0.174−0.137
 ξ33Price of labor−0.2520.015−0.281−0.223
 ξ34Price of materials0.2590.0320.1970.320
 δ31E(DD830)0.1750.115−0.0180.392
 δ32E(DD30)−0.0110.019−0.0480.028
 δ33E(PREC)0.1480.0680.0170.277
Demand for materials:
 λ4t0.0000.0000.0000.000
 ϕ4E(output price)0.0040.0040.0000.016
 ξ41Price of capital0.1490.0290.0920.206
 ξ42Price of land0.0250.0080.0090.040
 ξ43Price of labor0.1190.0120.0960.142
 ξ44Price of materials−0.2970.021−0.338−0.256
 δ41E(DD830)−0.2410.324−1.0700.099
 δ42E(DD30)−0.1410.023−0.181−0.087
 δ43E(PREC)0.0890.067−0.0430.225
Est.St. Err.2.5%97.5%
Demand for capital:
 λ1t−0.0050.001−0.007−0.004
 ϕ1E(output price)0.0170.0150.0000.056
 ξ11Price of capital−0.2520.026−0.302−0.202
 ξ12Price of land0.1450.0080.1310.160
 ξ13Price of labor−0.1650.011−0.187−0.143
 ξ14Price of materials0.2540.0250.2030.302
 δ11E(DD830)−1.4650.117−1.644−1.221
 δ12E(DD30)−0.1230.015−0.153−0.093
 δ13E(PREC)−0.1800.051−0.278−0.085
Demand for land:
 λ2t−0.0050.000−0.005−0.004
 ϕ2E(output price)0.0030.0030.0000.012
 ξ21Price of capital0.1710.0170.1380.204
 ξ22Price of land−0.1200.005−0.130−0.110
 ξ23Price of labor0.0030.007−0.0120.017
 ξ24Price of materials−0.0560.016−0.087−0.026
 δ21E(DD830)−0.0630.110−0.2560.120
 δ22E(DD30)−0.0560.011−0.077−0.035
 δ23E(PREC)0.2850.0340.2180.350
Demand for labor:
 λ3t−0.0080.001−0.009−0.006
 ϕ3E(output price)0.0130.0130.0000.046
 ξ31Price of capital0.1350.0330.0700.200
 ξ32Price of land−0.1550.010−0.174−0.137
 ξ33Price of labor−0.2520.015−0.281−0.223
 ξ34Price of materials0.2590.0320.1970.320
 δ31E(DD830)0.1750.115−0.0180.392
 δ32E(DD30)−0.0110.019−0.0480.028
 δ33E(PREC)0.1480.0680.0170.277
Demand for materials:
 λ4t0.0000.0000.0000.000
 ϕ4E(output price)0.0040.0040.0000.016
 ξ41Price of capital0.1490.0290.0920.206
 ξ42Price of land0.0250.0080.0090.040
 ξ43Price of labor0.1190.0120.0960.142
 ξ44Price of materials−0.2970.021−0.338−0.256
 δ41E(DD830)−0.2410.324−1.0700.099
 δ42E(DD30)−0.1410.023−0.181−0.087
 δ43E(PREC)0.0890.067−0.0430.225

4.3 Decomposing OSME index numbers

This section uses the parameter estimates reported in Table 4 (and Tables A4 to A7 in the Appendix) to decompose the OSME index numbers reported in Section 3.3. Recall that those numbers were computed by evaluating the third term on the right-hand side of Equation (11). Substituting Equation (12) (and a similar equation for firm k in period s) into that term yields the following:
(13)
Given the assumptions made earlier about emit, the first term on the right-hand side can be viewed as an input-oriented technology index (ITI), the second term can be viewed as an expected environment index (EEI), the third term can be viewed as an expected output price index (EPI), the fourth term can be viewed as an input price index (WI), and the last term can be viewed as an allocative efficiency and statistical noise index (AESNI).

Figure 3 presents a decomposition of OSME change in selected states. The OSME index numbers depicted in the six panels are the same numbers that were depicted earlier in Fig. 2, but using a different vertical scale. The ITI, EEI, EPI, and WI numbers were computed by using the parameter estimates reported in Table 4 (and Tables A4 to A7 in the Appendix) to evaluate the first four terms on the right-hand side of Equation (13). The AESNI numbers were computed as residuals. Results for all forty-eight states are presented in Table A8 in the Appendix.

OSME change in selected states from 1961 to 2004 (AL in 1961 = 1).
Figure 3.

OSME change in selected states from 1961 to 2004 (AL in 1961 = 1).

The results reported in Fig. 3 (and Table A8 in the Appendix) indicate that the main drivers of OSME change over time were technical progress and input price change. In the case of Alabama, for example, OSME in 2004 was 32.8% higher than it had been in 1961. This increase can be broken down as follows: OSMEI = ITI × EEI × EPI × WI × AESNI = 1.067 × 1.000 × 0.998 × 1.212 × 1.029 = 1.328. This decomposition indicates that, in Alabama, (i) technical progress provided for a 6.7% increase in OSME; (ii) changes in expected environmental conditions had no impact on OSME; (iii) changes in expected output prices had virtually no impact on OSME; (iv) changes in input prices led to a 21.2% increase in OSME; and (v) changes in allocative efficiency and statistical noise led to a 2.9% increase in OSME.

The results reported in Fig. 3 (and Table A8 in the Appendix) also indicate that the main drivers of OSME change across states were differences in input prices and expectations about environmental conditions. In 1961, for example, Texas farmers were 22% less scale-and-mix efficient than farmers in Alabama. This difference can be broken down as follows: OSMEI = ITI × EEI × EPI × WI × AESNI = 1 × 0.750 × 1.000 × 1.029 × 1.010 = 0.780. The reason the ITI component takes the value one is that, again, the model does not allow rates of technical progress to vary across states. The remaining components indicate that, in 1961, (i) differences in expected environmental conditions led Texas farmers to be 25% less scale-and-mix efficient than farmers in Alabama; (iii) differences in expected output prices made no difference to levels of scale-and-mix efficiency; (iv) differences in input prices incentivized Texas farmers to be 2.8% more scale-and-mix efficient than farmers in Alabama; and (v) levels of scale-and-mix efficiency in Texas were 4% higher than levels in Alabama due to differences in allocative efficiency and statistical noise.

Finally, it is possible to decompose the expected environmental change component of OSME change (i.e. the EEI) into the product of an input-oriented fixed effects index (IFEI) and an expected temperature and precipitation index (ETPI). The IFEI measures the effects of time-invariant environmental variables on OSME; mathematically, the IFEI that compares state i with state k is |$\prod _{m=1}^{M}\exp (\alpha _{mi}-\alpha _{mk})^{(\beta _m-\bar{s}_m)}$|⁠. The ETPI measures of the effects of expectations about temperature and precipitation on OSME; mathematically, the ETPI that compares expectations in state i in period t with expectations in state k in period s is |$\prod _{m=1}^{M}\prod _{j=1}^J (z^e_{jit}/z^e_{jks})^{\delta _{mj}(\beta _m-\bar{s}_m)}$|⁠. To illustrate, let us reconsider the finding that, in 1961, levels of OSME in Texas were 25% lower than levels in Alabama due to differences in expected environmental conditions. This difference can be broken down as follows: EEI = IFEI × ETPI =0.734 × 1.022 = 0.750. This indicates that time-invariant characteristics of the production environment (e.g. soil type and terrain) led Texas farmers to be 26.6% less scale-and-mix efficient than farmers in Alabama, and that, in 1961, differences in expectations about temperature and precipitation accounted for only 2.2% of the difference in OSME. Results for other states and other years also indicate that changes in long-term average temperature and precipitation (i.e. changes in climate) have a relatively small effect on OSME.

5 Summary and conclusion

Measuring productivity change involves measuring changes in output and input quantities. Economists tend to use chained Fisher (CF) and chained Törnqvist (CT) indexes to measure changes in quantities over time (see, for example, Schreyer, 2001, p.83). They also tend to use Elteto-Koves-Szulc (EKS) and Caves-Christensen-Diewert (CCD) indexes to measure changes in quantities across space (see, for example, Coelli et al. 2005, p.117). This paper measured changes in output and input quantities (and therefore productivity) in U.S. agriculture using geometric Young (GY) indexes. There are three reasons why GY indexes are preferred to CF, CT, EKS, and CCD indexes: first, they are proper indexes in the sense that they satisfy basic axioms from index theory; second, and relatedly, they yield numbers that are consistent with measurement theory; and, third, they can be used to make valid comparisons across both time and space.

Measuring productivity change is relatively straightforward. Explaining it is much more complicated because it involves explaining changes in output and input quantities. This paper explained changes in output and input quantities (and therefore productivity) using an economic model in which farmers choose outputs and inputs in two stages: first, at the beginning of the production period, farmers choose inputs and planned outputs to maximize expected profits in the face of uncertainty about output prices and one or more characteristics of the production environment (e.g. temperature); and, second, after inputs have been chosen, managers seek to maximize the outputs that can be produced using their chosen inputs in whatever production environment is realized. These behavioural assumptions gave rise to a stochastic production frontier model and a system of input demand equations.

This paper used the estimated parameters of the stochastic frontier model to decompose GY TFP index numbers into measures of technical progress, environmental change, technical efficiency change, scale-and-mix efficiency change, and changes in statistical noise. The main drivers of TFP change over time were found to be technical progress (i.e. the discovery of new techniques for transforming inputs into outputs) and scale-and-mix efficiency change (i.e. changes in how well farmers captured economies of scale and substitution). The main drivers of TFP change across states were found to be environmental change (i.e. changes in variables that are physically involved in the production process but never within the control of farmers) and, again, scale-and-mix efficiency change. The environmental change component was further decomposed into the product of time-invariant effects and weather effects. The results indicated that changes in weather had a relatively small effect on TFP. The policy implications are that there are two main ways in which governments can improve U.S. agricultural productivity: by improving rates of technical progress, and by increasing levels of scale and mix efficiency. Governments can improve rates of technical progress by conducting, or funding others to conduct, more R&D. To increase levels of scale-and-mix efficiency, they must first identify the drivers of scale-and-mix efficiency change.

This paper used the estimated parameters of the input demand equations to decompose the measure of scale-and-mix efficiency change into a measure of technical progress, a measure of input price change, various measures of changes in expectations, and a measure of changes in allocative efficiency and statistical noise. The main drivers of scale-and-mix efficiency change over time were found to be technical progress and input price change. The main drivers of TFP differences across states were found to be differences in input prices and expectations about environmental variables. The expected environmental change component was further decomposed into the product of time-invariant effects and expected weather effects. The results indicated that changes in expectations about temperature and precipitation (i.e. changes in climate) had a relatively small effect on scale and mix efficiency. The policy implications are that there are two main ways in which governments can improve scale and mix efficiency: by improving rates of technical progress, and by changing input prices. Many governments routinely change the prices of capital and labor by changing interest rates and the minimum wage, and they often change the prices of land and materials through taxes and subsidies.

Finally, this paper found that a small proportion of measured productivity change in U.S. agriculture could be attributed to changes in statistical noise. Statistical noise is a combination of functional form errors, omitted variable errors, measurement errors, and included variable errors. An avoidable source of statistical noise in this paper is errors in the data supplied by the USDA. Among other things, the USDA measures changes in output and input quantities by dividing revenue and cost indexes by output and input price indexes. So-called implicit quantity indexes are not proper indexes, and they do not yield numbers that are consistent with measurement theory. Economists at the USDA might want to investigate whether the use of proper quantity indexes is feasible, and, if so, whether it would lead them to change their views on spatial and temporal differences in agricultural output and input quantities (and TFP).

Footnotes

1.

Measures of TFP are measures of total output quantity divided by measures of total input quantity. Measures of multifactor productivity (MFP) and partial factor productivity (PFP) can be viewed as measures of TFP in which one or more inputs are given weight of zero. For this reason, and to avoid repetition, this paper focuses on measures of TFP.

2.

For readers who would like to double-check the computations, the prices of the apples (resp. oranges) in baskets A, B, C, D, and E are 1, 6, 5, 1, and 1 (resp. 1, 1, 2, 6, and 5).

3.

The output-oriented technical efficiency of state i in period t is OTEt(xit, qit, zit) = exp (−uit).

4.

The dependent variable in the Aigner et al. (1977) model is an output, not the logarithm of an output. This has implications for the interpretation of the technical inefficiency effect. For details, see O’Donnell (2018, p.326, footnote 1)

References

Aigner
D.
,
Lovell
C.
,
Schmidt
P.
(
1977
).
“Formulation and Estimation of Stochastic Frontier Production Function Models”
,
Journal of Econometrics
,
6
:
21
37
.

Anand
A.
,
Khetarpal
S.
(
2015
).
“Impact of Climate Change on Agricultural Productivity”
, in
Bahadur
B.
,
Raham
M.
,
Sahijram
L.
,
Krishnamurthy
K.
(eds).
Plant Biology and Biotechnology, vol. I: Plant Diversity, Organization, Function and Improvement
.
Springer India
.

Ball
V.
,
Bureau
J.-C.
,
Nehring
R.
et al. (
1997
).
“Agricultural Productivity Revisited”
,
American Journal of Agricultural Economics
,
79
:
1045
63
.

Coelli
T.
,
Rao
D.
,
O’Donnell
C.
et al. (
2005
).
An Introduction to Efficiency and Productivity Analysis
,
New York
:
Springer
, 2nd edn.

Geweke
J.
(
1992
).
“Evaluating the Accuracy of Sampling-based Approaches to Calculating Posterior Moments.”
, in
Bernado
J.
,
Berger
J.
,
Dawid
A.
,
Smith
A.
(eds).
Bayesian Statistics 4
.
Clarendon Press
.

Hughes
N.
,
Lawson
K.
,
Davidson
A.
et al. (
2011
).
“Productivity Pathways: Climate-adjusted Production Frontiers for the Australian Broadacre Cropping Industry”
,
Research Report 11.5
,
Australian Bureau of Agricultural and Resource Economics and Sciences
.

Jorgenson
D.
,
Griliches
Z.
(
1967
).
“The Explanation of Productivity Change”
,
The Review of Economic Studies
,
34
:
249
83
.

Nastis
S.
,
Michailidis
A.
,
Chatzitheodoridis
F.
(
2012
).
“Climate Change and Agricultural Productivity”
,
African Journal of Agricultural Research
,
7
:
4885
93
.

Njuki
E.
,
Bravo-Ureta
B.
,
Cabrera
V.
(
2020
).
“Climatic Effects and Total Factor Productivity: Econometric Evidence for Wisconsin Dairy Farms”
,
European Review of Agricultural Economics
,
47
:
1276
301
.

Njuki
E.
,
Bravo-Ureta
B. E.
,
O’Donnell
C.
(
2018
).
“A New Look at the Decomposition of Agricultural Productivity Growth Incorporating Weather Effects”
,
PLOS One
,
13
:
e0192432
.

O’Donnell
C.
(
2016
).
“Using Information about Technologies, Markets, and Firm Behaviour to Decompose a Proper Productivity Index”
,
Journal of Econometrics
,
190
:
328
40
.

O’Donnell
C.
(
2018
).
Productivity and Efficiency Analysis: An Economic Approach to Measuring and Explaining Managerial Performance
,
Singapore
:
Springer Nature
.

Plummer
M.
(
2019
).
rjags: Bayesian Graphical Models Using MCMC
,
R package version 4-10
.

Sabasi
D.
,
Shumway
C.
(
2018
).
“Climate Change, Health Care Access and Regional Influence on Components of U.S. Agricultural Productivity”
,
Applied Economics
,
50
:
6149
64
,
DOI:10.1080/00036846.2018.1489504
.

Salim
R.
,
Islam
N.
(
2010
).
“Exploring the Impact of R&D and Climate Change on Agricultural Productivity Growth: The Case of Western Australia”
,
Australian Journal of Agricultural and Resource Economics
,
54
:
561
82
.

Schreyer
P.
(
2001
).
OECD Manual: Measuring Productivity: Measurement of Aggregate and Industry-Level Productivity Growth
,
Paris
:
Organisation for Economic Development and Cooperation
.

Tal
E.
(
2016
).
“Measurement in science”
, in
Zalta
E. N.
(ed.)
The Stanford Encyclopedia of Philosophy
.
Metaphysics Research Lab, Stanford University
, winter 2016 edn.

This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]

Supplementary data