Next Article in Journal
Harvester Maintenance Resource Scheduling Optimization, Based on the Combine Harvester Operation and Maintenance Platform
Next Article in Special Issue
Technical Efficiency of Rice Production in the Upper North of Thailand: Clustering Copula-Based Stochastic Frontier Analysis
Previous Article in Journal
Potential of Silicon to Improve Biological Control of Fall Armyworm, Spodoptera frugiperda on Maize
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measurement and Analysis of Contribution Rate for China Rice Input Factors via a Varying-Coefficient Production Function Model

1
College of Mathematics and Informatics, South China Agricultural University, Guangzhou 510642, China
2
Key Laboratory of Smart Agricultural Technology in Tropical South China, Ministry of Agriculture and Rural Affairs, P.R.China, Guangzhou 510642, China
3
School of Information Engineering, Zhujiang College, South China Agricultural University, Guangzhou 510900, China
4
College of Engineering, South China Agricultural University, Guangzhou 510642, China
*
Author to whom correspondence should be addressed.
Agriculture 2022, 12(9), 1431; https://doi.org/10.3390/agriculture12091431
Submission received: 19 August 2022 / Revised: 5 September 2022 / Accepted: 7 September 2022 / Published: 9 September 2022

Abstract

:
To explore the internal driving force of the growth of rice yield per unit area in China, a model based on varying-coefficient production function is proposed in this study, which comes from the idea that the constant elasticity parameters in the Cobb-Douglas production function can be extended to functional forms. Applying such model to economic growth analysis, on the one hand, the dynamic contribution rate of each input factor can be measured, and, on the other hand, the contribution rate of the input factor can be decomposed into net factor contribution rate and interaction factor contribution rate, thus expanding the explanatory ability of growth rate equation. Using such model, the output elasticity of capital and labor in China’s rice yield growth are calculated from 1978 to 2020, and the dynamic characteristics of the contribution rate of capital, labor and generalized technological progress are analyzed. Next, the capital contribution rate is decomposed according to the composition of the total capital. The results show that: (1) The capital elasticity and labor elasticity are indeed not constant in different years. In China, from 1978 to 2020 the value of capital elasticity was between 0.3209 to 0.3589, with a mean of 0.3437, and the value of labor elasticity was between −0.1759 to −0.1640, with a mean of −0.1730. (2) Natural disasters do affect capital elasticity and labor elasticity in rice production. (3) When the annual proportion of crop disasters increases, the contribution rate of interaction between capital and natural disaster (KDR) value is negative, whereas the contribution rate of interaction between labor and natural disaster (LDR) value is positive. (4) Compared with 1978, the generalized technological progress contribution rate (GTPR) of the rice yield growth in China from 1979 to 2020 shows a declining trend in fluctuations, whereas the total capital contribution rate (TKR) shows a rising trend in fluctuations and the total labor contribution rate (TLR) is relatively stable in the same period. Since 2000, capital investment has become the main factor for the rice yield growth per unit area in China, of which machinery, chemical fertilizer, seed and pesticide are the four most important input factors.

1. Introduction

Rice is a food consumed regularly and is vital for the food security of over half the world’s population [1]. FAOSTAT [2] shows that from 1978 to 2020, the world’s rice harvested area increased from 143, 503, 572 to 164, 192, 164 ha, and the production increased from 385, 208, 660 to 756, 743, 722 tones. As a result, rice production plays an increasingly important role in world food security.
China is the largest rice producer and consumer in the world [3]. Rice is the staple food for over 65% of the Chinese population and so plays a significant role in Chinese food security [4,5]. In recent years, however, the cost of rice production in China has been rising, which continuously reduces margins, limits the development of the rice industry, reduces the enthusiasm of rice farmer for planting, and even affects food security [6]. Increasing rice yield is the main objective of rice research, the core goal of which is to improve the rice yield per unit area [7,8]. Although research on rice breeding and cultivation has greatly improved the potential rice yield, the actual rice yield is often lower than its theoretical value [9,10]. It is closely relevant to the technological efficiency and technological progress of rice production. Therefore, measuring the technological progress and technological efficiency of rice production is an important research area for economic analysis in the rice industry [11]. Particularly, from a more practical and operational perspective, researchers are more concerned about what factors play a key role in increasing rice yield per unit area. Generally speaking, capital, labor and technological progress are considered to be the three main input factors for the increase of rice yield per unit area. Therefore, measuring the contribution rate of input factors has become the basic task of the economic analysis of rice industry.
In the existing literature, there are many methods to estimate the contribution rate of technological progress, and the mainstream models include Cobb-Douglas (abbreviated to C-D) production function [12], data envelopment analysis [13], frontier production function [14], total factor productivity [15], etc. Data envelopment analysis is a nonparametric method. Using data envelopment analysis to measure the contribution rate of technological progress, the production frontier is measured by the method of mathematical programming, and the weight of input and output units is taken as the calculation variable, which avoids the defect of strong subjectivity in the other methods. Therefore, its advantages lie in the improvement of digital quantization and measurement accuracy. However, the disadvantage of data envelopment analysis is that this method has poor stability and is extremely sensitive to input and output data. The frontier production function model is a parametric method that has a wide range of applications in measuring the contribution rate of technological progress. Time series data can be used for long-term measurement, and time series and cross-sectional mixed data can also be used for short-term calculation. However, the model has the disadvantages of large sample demand and a complex estimation process. Total factor productivity is a method that is generally used to reflect the changing trend of the contribution rate of long-term technological progress [16]. C-D production function method is one of the earliest mathematical models used to measure economic growth. Due to the characteristics of simple calculation, easy operation and wide application range, it is commonly used to measure the contribution rate of technological progress. When using the C-D production function to measure the contribution rate of technological progress, it is necessary to estimate the output elasticity coefficient of each input factor in the C-D production function. However, these output elasticity coefficients are estimated to be fixed constants, representing the average elasticity over the research period. Using them directly can only adapt to the analysis of the phase changes in the whole study period, they cannot adapt to the analysis of the dynamic changes phase by phase.
The main purpose of this study is to: (1) propose a varying-coefficient production model, and (2) measure the output elasticity of input factors and analyze the dynamic characteristics of the contribution rate of input factors using the proposed model.

2. Materials and Methods

In this Section, we firstly give the construction process of the varying-coefficient production function model. Secondly, we introduce how to solve the varying-coefficient regression model with two different covariate variables by using an averaged estimate method. Then, bandwidth selection is discussed. Finally, the measurement methods of contribution rate of input factors are derived.

2.1. Proposed Model

Productivity analysis is a main tool adopted to explore the source of economic growth, and the measurement of the growth rate of total factor productivity is the central content of productivity analysis [17,18]. Among this, the main challenge is to determine the contribution rate of different input factors in economic growth quantitatively. Production function is a mathematical expression that describes the dependency relationship between a certain combination of production factors and the possible maximum output in the production process [19]. It is typically used to analyze the process of economic growth and measure the contribution rate of various economic growth factors to economic growth. The general form is
Y = f ( K , L , A )
where Y, K, L and A represent real output, capital investment, labor input and technical level, respectively. Since Solow’s famous work in 1957 [20] on measuring technological progress, also known as total factor productivity, the C-D production function has become the most widely used form of studying technological progress, due to its simple structure, clear economic significance and easy estimation. It can be expressed as
Y = A K α L β
where α and β represent the elasticity of capital investment and labor force, respectively.
In model (2), α and β are fixed constants that both reflect the average elasticity over the whole research period. Using them directly can only adapt to the analysis of the phase changes in the whole study period and cannot adapt to the analysis of the dynamic changes phase by phase. However, objectively speaking, various input factors in real production are constantly changing, especially with the continuous advancement of the market-oriented reform process, the impact of various production factors on output growth should show a dynamic change law. Therefore, the production function model with variable output elasticity is worth studying [19].
The varying-coefficient model is a useful generalization of the classical linear model, originating from Shumway’s monograph it was then systematically studied by Cleveland, Hastie and Tibshirani [21,22,23]. Currently, the widely used model is the one proposed by Cai et al. [24], and its specific form is written as follows:
Y = j = 1 p a j ( U ) X j + ε
where Y denotes the dependent variable, ( U ,   X 1 , , X p ) denotes independent variables, a j ( · ) ( j = 1 ,   , p ) are some unknown functions, ε is the random error satisfying E ( ε | U ,   X 1 , , X p ) = 0 and V a r ( ε | U ,   X 1 , , X p ) = σ 2 ( U ) . Such a model arises naturally when one wishes to examine how regression coefficients are changed with certain specific independent variables such as time, temperature, area and more.
Based on varying-coefficient modeling techniques, Ahmad et al. [25], Luo et al. [26] and Zhang et al. [19,27,28] introduced the varying-coefficient model to estimate the variable elasticity of the capital and labor force. In Luo’s literature, the model is supposed as Y = K α ( t ) L β ( t ) , and is also written as ln Y = α ( t ) ln K + β ( t ) ln L + ε , where t is time, ε is the random error, α ( t ) and β ( t ) denote the elasticity coefficient function of the capital and the labor force, respectively. The model is quite simple but unreasonable, because the technical level A(t) is not included in the model. Of course, the model cannot be assumed to be in the form of Y = A ( t ) K α ( t ) L β ( t ) , because of it is indiscernibility in statistics, as discussed in Luo’s literature [26]. Therefore, in the articles of Zhang et al. [19,27], the model is written as ln Y = i = 1 m γ i Z i + α ( t ) ln K + β ( t ) ln L + ε , where ln A = i = 1 m γ i Z i . Of which, it is assumed that the technical level can be expressed as an exponential linear combination of a set of controllable variables. In essence, it is a semi-parametric varying-coefficient regression model.
In the above models, it is assumed that the output elasticity of capital and labor is a nonparametric smooth function that changes with time. Although they can reflect the time-varying characteristics of elasticity, they cannot explain why they change with time. Since “time” is a concept defined by people, it is not the fundamental reason for the change of capital elasticity and labor elasticity. In the practice of rice production, it can be considered that the technical level of rice production in two adjacent years is roughly the same. Assuming that the capital and labor input of a certain region are the same for the two adjacent years, theoretically, the rice yield per unit area should be similar, but the actual rice yield per unit area will be different. Sometimes, significant differences exist, the reasons for which are very complex and not completely clear. We believe that natural disasters may be one of the main causes. Moreover, exploring the influence of climate change on rice yield has received extensive attention [29,30]. Therefore, we assume that the output elasticity is a nonparametric smoothing function dependent on natural disaster factors (such as drought, flood, storm, low temperature, cold damage and so on). Although it is assumed that natural disasters are independent of time in the later part of this study, the numerical value of natural disasters in different years changes with time. Therefore, there are still time-varying characteristics in output elasticity.
To establish an appropriate production function model for rice yield growth analysis, the following assumptions will be made.
(1)
Natural disasters are an important factor that affect the technical efficiency of rice production. Therefore, they are also a covariant factor that affects the changes of rice capital elasticity and labor elasticity. A natural disaster is represented by the annual proportion of crop disaster Z. It is assumed that the output elasticity is a smooth function dependent on the proportion of crop disaster.
(2)
Generally, the technical level, A, is neutral, which essentially reflects the impact of all other factors except for capital and labor inputs on output growth. Therefore, it is assumed that the technical level, A, does not depend on the change due to natural disasters, but only on the change of time, reflecting its dynamic characteristics.
(3)
Supposing the natural disaster factor and the time factor are independent. Such an assumption suggests that the occurrence of natural disasters is completely random.
Based on the above three assumptions, the following varying-coefficient production function model is constructed:
Y = A ( t ) K α ( Z ) L β ( Z )
Take the natural logarithm on both sides and add a random error term to obtain:
ln Y = A 0 ( t ) + α ( Z ) ln K + β ( Z ) ln L + ε
where Y, K, L, Z and t represent rice yield, capital investment, labor input, annual proportion of crop disaster and time, respectively. A 0 ( t ) = ln A ( t ) reflects generalized technical progress. ε is the random error, satisfying that the mean is 0 and the variance is σ 2 .

2.2. Estimation Method

In model (5), the intercept function A 0 ( t ) and the coefficient functions α ( Z ) and β ( Z ) are two types of function with different smoothing variables. The estimation method will be different from the traditional varying-coefficient model where all coefficient functions share a single smoothing variable in a model, due to that the commonly used kernel estimation methods or the local polynomial estimation method are not applicable here. At present, there are two approaches to solve the problem. One is the marginal integration technique [31,32,33], the other is the smooth backfitting [34,35].
In this study, the idea of averaged estimation [31,32] will be adopted to construct the estimation method, which is essentially a marginal integration technique. The basic procedure is described as follows:
Supposing that { Y ( t i ) , K ( t i ) , L ( t i ) , Z ( t i ) , t i ,   i = 1 , , n } is a sample from model (5), A 0 ( t ) , α ( Z ) and β ( Z ) all have a Lipschitz continuous second derivative, then A 0 ( t ) can be locally approximated by a linear function at a neighborhood of t i , in which t i belongs to the support of t , and α ( Z ) , β ( Z ) can be locally approximated by a linear function at a neighborhood of Z j = Z ( t j ) , in which Z j belongs to the support of Z . The expressions are written as follows:
A 0 ( t ) A 0 ( t i ) + A 0 ( t i ) ( t t i ) = def A 0 + A 1 ( t t i ) ,
α ( Z ) α ( Z j ) + α ( Z j ) ( Z Z j ) = def α 0 + α 1 ( Z Z j ) ,
β ( Z ) β ( Z j ) + β ( Z j ) ( Z Z j ) = def β 0 + β 1 ( Z Z j ) .
where A 0 ( t i ) denotes the derivative of A 0 ( t ) at t i , α ( Z j ) and β ( Z j ) denote the derivatives of α ( Z ) and β ( Z ) at Z j = Z ( t j ) . We minimize
i = 1 n j = 1 n { ln Y ( t i ) [ A 0 + A 1 ( t t i ) ] [ α 0 + α 1 ( Z Z j ) ] ln K ( t i ) [ β 0 + β 1 ( Z Z j ) ] ln L ( t i ) } 2 H h 1 ( t t i ) s = 2 3 H h s ( Z Z j ) ,
where H h ( ) = h 1 H ( / h ) , H ( ) is a bounded, nonnegative, compactly supported symmetric about zero and Lipschitz continuous density function. This study uses the Epanechnikov kernel function H ( x ) = 0.75 ( 1 x 2 ) + . h i ( i = 1 , 2 , 3 ) is a sequence of positive numbers, and is called bandwidth.
Let
Y i * = ln Y ( t i ) , K i * = ln K ( t i ) , L i * = ln L ( t i ) , Y * = ( Y 1 * , , Y n * )
θ = ( A 0 , A 1 , α 0 , α 1 , β 0 , β 1 ) T ,   W i , j = H h 1 ( t t i ) s = 2 3 H h s ( Z Z j ) , W = d i a g ( W 1 , 1 , , W n , n ) ,
X = ( 1 1 t 1 t i t n t i K 1 * K n * ( Z 1 Z j ) K 1 * ( Z n Z j ) K n * L 1 * L n * ( Z 1 Z j ) L 1 * ( Z n Z j ) L n * ) .
where the superscript T represents the transpose of a vector or a matrix, and W = d i a g ( W 1 , 1 , , W n , n )   is a diagonal matrix with W 1 , 1 , , W n , n being its diagonal entries.
Based on the local polynomial estimation and the principle of least squares theory, solving the minimization of Equation (9), we can get
θ ˜ ( t i , Z j ) = ( X T W X ) 1 X T W Y * .
In this study, θ ˜ is called an initial value of θ . That is, let A 0 ˜ ( t i ) , α ˜ ( Z j ) , β ˜ ( Z j ) denote an initial value of A 0 ( t ) , α ( Z ) , β ( Z ) , respectively, and then
A 0 ˜ ( t i ) = A 0 ˜ = ( 1 , 0 , 0 , 0 , 0 , 0 ) θ ˜ ( t i , Z j ) ,
α ˜ ( Z j ) = α 0 ˜ = ( 0 , 0 , 1 , 0 , 0 , 0 ) θ ˜ ( t i , Z j ) ,
β ˜ ( Z j ) = β 0 ˜ = ( 0 , 0 , 0 , 0 , 1 , 0 ) θ ˜ ( t i , Z j ) .
According to the idea of averaged estimation, the estimators of the intercept function A 0 ( t ) and the coefficient functions α ( Z ) , β ( Z ) are further defined as follows:
A 0 ^ ( t i ) = ( 1 , 0 , 0 , 0 , 0 , 0 ) 1 n j = 1 n θ ˜ ( t i , Z j ) ,
α ^ ( Z j ) = ( 0 , 0 , 1 , 0 , 0 , 0 ) 1 n i = 1 n θ ˜ ( t i , Z j ) ,
β ^ ( Z j ) = ( 0 , 0 , 0 , 0 , 1 , 0 ) 1 n i = 1 n θ ˜ ( t i , Z j ) .

2.3. Bandwidth Selection

According to the local polynomial estimation theory, the bandwidth has a great influence on the estimation of the coefficient functions. Generally, if the bandwidth is small, the bias of the estimation of coefficient function is also small, but the variance is large. On the contrary, if the bandwidth is larger, the variance of the estimation of coefficient function is smaller, but the bias is larger. In model (5), due to the intercept function A 0 ( t ) and coefficient functions α ( Z ) and β ( Z ) depend on different smoothing covariates respectively, it is difficult to get a good estimation by selecting a single bandwidth. The existing literature shows that in order to make the estimations having good properties, such as asymptotic normality and more, the bandwidth needs to meet the following constraints [31,32,33]:
1)
n i = 1 p h i , n ;
2)
h i = c i n 1 5 , n h j 5 0 , j i ,   i , j = 1 , , p , and c i is a constant.
However, so far, how to select the optimal bandwidth in varying-coefficient models with different smoothing covariates remain a problem that needs further research. In this study, we adopted Zhang’s suggestions [32]: Supposing our aim is to estimate the coefficient function α ( Z j ) , let h j = c j n 1 5 , the bandwidths corresponding to other smoothing variables are set as h i = c i n 3 5 ( p 1 ) , i j , where i , j = 1 , , p ,   c i   a n d   c j are undetermined constants, p is the number of smoothing variables. In this study, p = 3. These undetermined constants, c i   and c j , are determined using the cross-validation method.

2.4. Measurement Methods on Contribution Rate of Input Factors

Firstly, the growth rate equation under the varying-coefficient production function model is derived. Supposing that { Y ( t i ) , K ( t i ) , L ( t i ) , Z ( t i ) , t i ,   i = 1 , , n }   is a time series sample, the model (5) will be written as
ln Y ( t ) = A 0 ( t ) + α ( Z ( t ) ) ln K ( t ) + β ( Z ( t ) ) ln L ( t ) + ε ( t ) , t = t i , i = 1 , , n .
where Z ( t ) is a virtual function used to derive the calculation formulas for measuring the contribution rates of input factors. In model (17), taking the derivative with respect to t on both sides, we can obtain
d Y ( t ) Y ( t ) d t = d A 0 ( t ) d t + d α ( Z ( t ) ) d Z ( t ) d Z ( t ) d t ln K ( t ) + α ( Z ( t ) ) d K ( t ) K ( t ) d t     + d β ( Z ( t ) ) d Z ( t ) d Z ( t ) d t ln L ( t ) + β ( Z ( t ) ) d L ( t ) L ( t ) d t + d ε ( t ) d t .
Equation (18) is called the time point growth rate equation. In practice, an interval growth rate equation is needed due to the fact that the sample is often time series data. In other words, the differentials in Equation (18) needs to be approximately replaced by the corresponding increment. In general, considering the period [ t , t + 1 ] , define
d t = Δ t = t + 1 t = 1 ,
d Y ( t ) = Δ Y ( t ) = Y ( t + 1 ) Y ( t ) ,
d K ( t ) = Δ K ( t ) = K ( t + 1 ) K ( t ) ,
d L ( t ) = Δ L ( t ) = L ( t + 1 ) L ( t ) ,
d A 0 ( t ) = Δ A 0 ( t ) ,
d α ( Z ( t ) ) d Z ( t ) d Z ( t ) d t = d α ( Z ( t ) ) d t = Δ α ( Z ( t ) ) = α ( Z ( t + 1 ) ) α ( Z ( t ) ) ,
d β ( Z ( t ) ) d Z ( t ) d Z ( t ) d t = d β ( Z ( t ) ) d t = Δ β ( Z ( t ) ) = β ( Z ( t + 1 ) ) β ( Z ( t ) ) ,
d ε ( t ) = Δ ε ( t ) .
In Equation (18), α ( Z ( t ) ) and β ( Z ( t ) ) are the point elasticity of K ( t ) and L ( t ) with respect to Y ( t ) , respectively. When the time point growth rate equation is replaced by an interval growth rate equation, the point elasticities in Equation (18) need to be replaced by the arc elasticities. There are two methods for estimating the arc elasticities on the period [ t , t + 1 ] as follows.
Method 1: Let α ˜ ( Z ( t ) ) = [ α ( Z ( t ) ) + α ( Z ( t + 1 ) ) ] / 2 , β ˜ ( Z ( t ) ) similar definition.
Method 2: Let α ˜ ( Z ( t ) ) = w α ( Z ( t ) ) + ( 1 w ) α ( Z t + 1 ) , β ˜ ( Z ( t ) ) similar definition.
where w denotes the weight, α ˜ ( Z ( t ) ) and β ˜ ( Z ( t ) ) are the arc elasticities of K ( t ) and L ( t ) with respect to Y ( t ) , respectively. Obviously, method 1 is a special case of method 2, in which w = 0.5 .
In Equation (18), due to the fact that the values of α ( Z ) and β ( Z ) are estimated, the estimation errors of α ( Z ) and β ( Z ) will be inevitable. Therefore, the error item d ε ( t ) exists naturally. According to Solow’s idea, the error term is included in the so-called “residual value”, we define Δ A 0 ˜ ( t ) = Δ A 0 ( t ) + Δ ε ( t ) and this is used to reflect the share of output growth contribution by all factors except K ( t ) and L ( t ) . That is, the so-called generalized technological progress.
To sum up, the interval growth rate equation corresponding to Equation (18) is written as
Δ A 0 ˜ ( t ) = Δ Y ( t ) Y ( t ) [ Δ α ( Z ( t ) ) l n K ( t ) + α ˜ ( Z ( t ) ) Δ K ( t ) K ( t ) + Δ β ( Z ( t ) ) l n L ( t ) + β ˜ ( Z ( t ) ) Δ L ( t ) L ( t ) ] .
where Δ A 0 ˜ ( t ) is called generalized technological progress, and Δ Y ( t ) Y ( t ) is the rate of output growth. α ˜ ( Z ( t ) ) Δ K ( t ) K ( t ) is called net capital contribution, which reflects the increase of output caused purely by an increase in capital input; Δ α ( Z ( t ) ) l n K ( t ) is called interaction between capital and natural disaster, which reflects the change of capital contribution under the influence of natural disasters; and Δ α ( Z ( t ) ) l n K ( t ) + α ˜ ( Z ( t ) ) Δ K ( t ) K ( t ) is called total capital contribution, which comprehensively reflects the contribution of capital to output growth under the influence of certain natural disasters. Obviously, when suffering from a serious natural disaster, Δ α ( Z ( t ) ) < 0   and Δ α ( Z ( t ) ) l n K ( t ) < 0 will appear. At this time, the interaction of capital and natural disaster has a negative impact on output growth, the total capital contribution will be smaller than the net capital contribution and vice versa. Similarly, β ˜ ( Z ( t ) ) Δ L ( t ) L ( t ) is called net labor contribution, which reflects the increase of output caused purely by an increase in labor input; Δ β ( Z ( t ) ) l n L ( t ) is called interaction between labor and natural disaster, which reflects the change of labor contribution under the influence of natural disasters; and Δ β ( Z ( t ) ) l n L ( t ) + β ˜ ( Z ( t ) ) Δ L ( t ) L ( t ) is called total labor contribution, which comprehensively reflects the labor contribution to output growth under the influence of certain natural disasters.
According to Equation (19), the calculation formulas are defined as follows:
The net capital contribution rate (abbreviated to NKR)
NKR = [ α ˜ ( Z ( t ) ) Δ K ( t ) K ( t ) / Δ Y ( t ) Y ( t ) ] × 100 %
The contribution rate of interaction between capital and natural disaster (abbreviated to KDR)
KDR = [ Δ α ( Z ( t ) ) l n K ( t ) / Δ Y ( t ) Y ( t ) ] × 100 %
The total capital contribution rate (abbreviated to TKR)
TKR = NKR + KDR
The net labor contribution rate (abbreviated to NLR)
NLR = [ β ˜ ( Z ( t ) ) Δ L ( t ) L ( t ) / Δ Y ( t ) Y ( t ) ] × 100 %
The contribution rate of interaction between labor and natural disaster (abbreviated to LDR)
LDR = [ Δ β ( Z ( t ) ) l n L ( t ) / Δ Y ( t ) Y ( t ) ] × 100 %
The total labor contribution rate (abbreviated to TLR)
TLR = NLR + LDR
The generalized technological progress contribution rate (abbreviated to GTPR)
GTPR = 1 TKR TLR

3. Empirical Analysis and Results

3.1. Data Source

In this study, we will use the national rice production data of China for empirical analysis. The data is extracted from China Statistical Yearbook (1979–2021) and China Agricultural Product Cost-Benefit Compilation (1979–2021).
The data used in this study includes rice yield per unit area Y i (Kg•hm−2), material and service cost per unit area of rice K i (CNY•hm−2), labor employment per unit area of rice L i (Day•hm−2) and annual proportion of crop disasters Z i (%) from 1978 to 2020 in China. The time is year denoted by t i , i = 1978 , , 2020 . From the perspective of cost-benefit statistics, rice production costs include material and service costs and labor costs. However, China’s rice production is still based on small-scale decentralized production and management, and it is difficult to accurately measure the labor costs. In China Agricultural Production Cost-Benefit Compilation, the labor costs are mainly composed of the discount of family labor. To avoid double counting the discount of family labor into the model, the material and service cost per unit area and labor employment per unit area in rice represent capital investment and labor input, respectively. The natural disaster is represented by the annual proportion of crop disasters. The value of material and service costs per unit area are adjusted by the price index of means of agricultural production, in which 1977 is the base period.

3.2. Results of Estimation

In the varying-coefficient production function model, the optimal kernel function the Epanechnikov kernel function H ( z ) = 0.75 ( 1 z 2 ) + is selected. The optimal bandwidths are h 1 , o p t = 0.5383 , h 2 , o p t = 0.1495 , h 3 , o p t = 0.1497 , respectively, determined by the cross-validation method. Based on the optimal bandwidths, the annual estimated values of capital elasticity and labor elasticity from 1978 to 2020 are obtained (See Appendix A Table A1). The results show that from 1978 to 2020, the value of capital elasticity is between 0.3209 and 0.3589, with mean 0.3437, and the value of labor elasticity is between −0.1759 and −0.1640, with mean −0.1730.
In this study, the key problem is the estimation of capital elasticity and labor elasticity. To test the reliability of our estimation, we use the same data sets to estimate and compare the capital elasticity and labor elasticity with the varying-coefficient production function model and the C-D production function model, respectively. The data sets are 1978–2016, 1978–2017, 1978–2018, 1978–2019 and 1978–2020, respectively. In the C-D production function model, natural disaster data is missing. The results are shown in Table 1.
It is observed that the average estimation results of capital elasticity and labor elasticity of the two models are different because of the different construction mechanisms of the models, but the change trends of capital elasticity and labor elasticity are basically the same. Such trends can be explained in other words that, in the past five years, the capital elasticity has an increasing trend, whereas the absolute value of labor elasticity has a slightly decreasing trend. This shows that the application of the varying-coefficient production function model to measure capital elasticity and labor elasticity is worthy of reference. Here, the labor elasticity is negative, indicating that labor input has a negative impact on the growth of rice yield per unit area, and China’s rice production is in the stage of declining marginal return from the labor force.
In fact, using the C-D production function to measure capital elasticity and labor elasticity severely depends on the data set and model variable structure. Different data sets give different elastic estimation. For example, using the C-D production function, Tian et al. [36] obtained that the capital elasticity was 0.1072 and the labor elasticity was −0.0631 based on China’s rice production data from 1991 to 2008. Li et al. [37], based on China’s rice production data from 1978 to 2012, concluded that the capital elasticity was 0.4824 and the labor elasticity was −0.0826. Based on the rice production data of the Sichuan Province from 1980 to 2015, Chen et al. [38] attained that the capital elasticity was 0.2814 and the labor elasticity was −0.0655. As we all know, the C-D production function model is a parametric model, and the varying-coefficient model is a nonparametric model. In statistics, nonparametric models are generally more robust than parametric models. Moreover, using the mean square error (MSE) to evaluate the fitting of the model, we find that based on the rice production data in China from 1978 to 2020, the MSE of the estimation in the varying-coefficient production model is 2.23 × 10−5, less than the MSE of the estimation in the C-D production function model, which is 1.54 × 10−3. To sum up, the use of the varying-coefficient production function model in measuring capital elasticity and labor elasticity depends on the data set, but also the fitting effect and robustness of the model; the varying-coefficient production function model is better than the C-D production function model.
To show the estimated effect, the fitting results of the model are showed in Figure 1, where LnY represents the natural logarithm of Y. The dynamic characteristics of output elasticity of capital and labor are presented visually in Figure 2 and Figure 3, respectively.
In Figure 1, the real values and the fitting values of LnY are shown with asterisks and circles, respectively. The difference between the real values and the fitting values is quite small, thus the estimation of the model is satisfactory. Figure 2 and Figure 3 show that capital elasticity and labor elasticity have obvious fluctuation characteristics in different years. From 1978 to 2004, both capital elasticity and labor elasticity shows an obvious downward trend. From 2005 to 2020, capital elasticity shows an increasing trend, whereas labor elasticity is relatively stable. When the proportion of crop disasters is significantly larger, such as in the years 1980, 1991, 1994, 1997, 2000, 2001 and 2003, the elasticity of capital and labor decrease significantly. It indicates that natural disasters greatly effect rice yield per unit area, as well as the capital elasticity and labor elasticity. As a result, hypothesis (1) in Section 1 is reasonable.
To explore the impact of natural disasters on the output elasticity of capital and labor, scatter plots are drawn with the annual proportion of crop disasters as abscissa and the corresponding elasticity value as ordinate. The results are given in Figure 4 and Figure 5.
It can be seen from Figure 4 that there is a significant negative linear correlation between capital elasticity and annual proportion of crop disasters, with a correlation coefficient of −0.9817. Figure 5 shows that there is also a negative correlation between the absolute value of labor elasticity and annual proportion of crop disasters, with a correlation coefficient of −0.8752. As the annual proportion of crop disasters increases, the decreasing speed of capital elasticity and labor elasticity tends to increase. When the annual proportion of crop disasters is more than about 18%, the elasticity of capital and labor show a significant decline curve. It indicates that capital elasticity and labor elasticity are greatly affected by natural disasters. When natural disasters reach a certain proportion, such as 18%, the capital elasticity and the labor elasticity will decrease rapidly, and at this time, the ability of existing rice production technology to resist natural disasters will decrease rapidly.

3.3. Contribution Rate of Rice Input Factors

To investigate the dynamic characteristics of the contribution rate of input factors to the rice yield growth in China, taking the production level in 1978 as the common comparison object, the contribution rate of each input factor in rice yield growth from 1979 to 2020 are calculated by using the measurement methods described in Section 2.3. The results are shown in Appendix A Table A2. Figure 6 shows the change trend of the contribution rate of the three main production factors to the rice yield growth. The names of the three main contribution rates are total capital contribution rate (abbreviated to TKR), total labor contribution rate (abbreviated to TLR) and generalized technical progress contribution rate (abbreviated to GTPR). We emphasize that this study does not calculate the contribution rate in rice yield growth each year relative to the previous year, because technological progress is a gradual process, and there is little difference in the technical level of rice production between two adjacent years. It is difficult and meaningless to calculate the technological progress of each year, compared with the previous year in rice production.
Figure 6 shows that, compared with 1978, the GTPR of the rice yield growth in China from 1979 to 2020 shows a declining trend in fluctuations, whereas the TKR shows a rising trend in fluctuations and the TLR is relatively stable in the same period. Specifically, from 2000 to 2020, the GTPR decreased significantly. Its value decreased from 52.58% in 2000 to 15.30% in 2020. The value of the TKR increased from 25.80% in 2000 to 65.52% in 2020, an increase of 153.95%. Capital investment became the main factor in rice yield increase. It shows that China’s rice yield growth in recent years has mainly depended on the increase of material and service cost, indicating that it is an extensive and unsustainable growth mode.

3.4. Decomposition of Capital Contribution Rate

It can be seen from Section 3.3 that the increase in rice yield per unit area in China has mainly depended on the increase of capital investment since 2000. To identify which inputs are playing the leading role, this section will decompose the total capital contribution rate for the second time.
In this study, capital investment is measured by the material and service cost per unit area of rice. The material and service cost are composed of seed cost, chemical fertilizer cost, farm fertilizer cost, pesticide cost, agricultural film cost, lease cost (including mechanical operation cost, irrigation and drainage cost, animal power cost, etc.), fuel power cost, technical service cost, tools and materials cost, repair and maintenance cost, other direct costs and indirect costs. We define the weight of each input is the proportion of the cost of each input factor to the material and service cost. The capital contribution rate will be decomposed according to the weight of the input factors. The data is extracted from China Agricultural Product Cost-Benefit Compilation. The results are listed in Appendix A Table A3, and the trend for visual perception is shown in Figure 7.
Figure 7 shows that the top four input factors are machinery, fertilizer, seed and pesticide in order of contribution rate. The contribution rate of machinery increases fastest, and machinery becomes the factor with the largest contribution rate among capital investment. The contribution rate of machinery increased from 2.76% in 2000 to 24.24% in 2020. It proves that the mechanization of rice production in China has made great progress. The second most important input factor is chemical fertilizer. The contribution rate of chemical fertilizer increased from 7.38% in 2000 to 18.25% in 2012, and slightly decreased from 17.56% in 2013 to 16.47% in 2020 ranking it second. The third important input factor is seed. The contribution rate of seed increased from 1.98% in 2000 to 8.18% in 2020. The fourth input factor is pesticide. The contribution rate of pesticide increased from 1.93% in 2000 to 7.35% in 2020. Thus, it can be seen that, firstly, since the 21st century, with the rapid growth of agricultural labor costs, machinery plays an increasingly important role in China’s rice production, and the development of rice production mechanization is of great significance to improve the rice yield. Secondly, chemical fertilizer and pesticide play an important role in the rice yield growth, indicating that the current rice production mode in China is not sustainable. In the future, people expect to drastically reduce the use of fertilizers and pesticides in rice production. Finally, the most fundamental way to improve rice yield per unit area is to cultivate good varieties. The application of superior rice varieties has played an increasingly important role in China’s rice yield growth. Based on the improvement of rice varieties and the integration of agro-machinery and agronomy, China’s rice production mode will quickly enter the modernization stage.

4. Discussion

4.1. Policy Implications

The contribution rate of rice input factors is an important basis for the formulation of rice industrial policy. The empirical results presented in this study provide important insights into policy design in the rice industry.
One of the main contributions of this study is to provide a method for measuring the dynamic elasticity and contribution rate of input factors. Using such method, the characteristics of the contribution rate of rice input factors in China from 1979 to 2020 were investigated, and the conclusions from which can be used as a reference for the adjustment of rice industrial policy.
It can be seen from the conclusion in Section 3.2 that, compared with results calculated by our proposed model, the capital elasticity is overestimated, whereas the labor elasticity is underestimated, when using the traditional C-D production function to measure capital elasticity and labor elasticity. As a result, the capital contribution rate is overestimated, the labor contribution rate is underestimated, and the contribution rate of generalized technological progress is distorted. The results in this study show that the contribution of generalized technological progress to China’s rice production presents a downward trend in fluctuations from 1979 to 2012, which is basically consistent with the conclusion by Li et al. (2016) [37]. However, from 2012 to 2020, the downward trend of the contribution rate of generalized technological progress has been significantly alleviated, the labor contribution rate has become more stable, and the capital contribution rate has a slight increasing trend, but the growth rate has slowed down significantly. The results show that since 2012, China’s rice industrial policy has achieved certain success in changing the extensive mode of rice production. In addition, it can be seen from Figure 2 and Figure 3 that the watershed between the dynamic changes of capital elasticity and labor elasticity was 2004. In the same year, the Law of the People’s Republic of China on Promotion of Agricultural Mechanization was promulgated and implemented. Perhaps, such a discovery can be used as one of the bases to test the impact of the Law of the People’s Republic of China on Promotion of Agricultural Mechanization on China’s rice industry.
In the analysis of economic growth, capital is a complex concept. The secondary decomposition of capital contribution rate can measure the contribution rate of various material inputs, thus providing a more targeted basis for the formulation of industrial policies [39]. The results in this study show that, since 2000, although China has made great progress in rice production technology, including the success of super hybrid rice breeding and the gradual maturity of cultivation technology and mechanization technology, the blind expansion of capital investment has reduced the contribution of technological progress to the growth of rice yield per unit area. Among the numerous capital investments in rice production, machinery, chemical fertilizer, seed and pesticide are the four most important input factors. The increase in machinery and seed input at this stage is in line with China’s needs to develop rice mechanization and strengthen variety cultivation. However, the increase of chemical fertilizers and pesticides has restricted the sustainable development of China’s rice industry. With the development and progress of human society, not starving is no longer the only goal of human life, and grain yield is also no longer the only goal of agricultural production. Healthy and delicious food is gradually getting people’s favor. The environmental pollution and food toxicity caused by the residues of chemical fertilizers and pesticides are also being realized and studied [40,41]. In future rice production, we cannot unilaterally pursue the improvement of yield at the cost of environmental pollution and food poisoning. The extensive mode of rice production in China has not been fundamentally changed. Technology and resource constraints remain the “key” issues that perplex the development of China’s rice industry [42]. Strengthening scientific and technological support as well as optimizing the policy system are the inevitable choice for the development of China’s rice industry.

4.2. Advantage of the Proposed Model and Applications in Future Work

Based on the idea of varying-coefficient modeling, the C-D production function is extended to a varying-coefficient production function, which can be used to estimate the dynamic elasticity of input factors and calculate the contribution rate of input factors, providing a new method to study the dynamic characteristics of the contribution rate of input factors in economic growth.
In particular, the analytical ability of the model is expanded due to the introduction of the natural disaster factor. In this model, the contribution rate of an input factor can be decomposed into the net factor contribution rate and the contribution rate of interaction between input factor and natural disaster. For example, it can be seen from Appendix A Table A2 that the natural disaster has a certain impact on the rice yield growth. The value of KDR is between −31.25% and 41.13%, with mean 1.11%. The value of LDR is between −11.85% to 28.98%, with mean 0.61%. The law is concluded as follows: when the annual proportion of crop disaster increases, the KDR is negative, whereas the LDR value is positive and vice versa. This is consistent with intuitive perception. Therefore, when the traditional C-D production function is extended by varying-coefficient production function, the new model can not only be used to measure the dynamic output elasticity and contribution rate of input factors in the study period, but also can be used to explore the interaction between input factors and natural disasters in studying rice yield growth. Such a model can provide a new way for in-depth analysis of the factors in studying economic growth.
Of course, the model can be further improved. When the model is applied, its adaptability needs to be considered. Firstly, in essence, the varying-coefficient production function model belongs to the nonparametric statistical model. To get good results, large samples are generally required. Secondly, the bandwidth has a great influence on the estimation of the coefficient functions. How to determine the optimal bandwidth is still a problem to be further studied [31,32]. Thirdly, the variable, Z, in the elastic functions α ( Z ) and β ( Z ) is not unique. In this study, Z means the annual proportion of crop disasters. In other applications, Z can be represented by other variables such as temperature, concentration, humidity and so on.

5. Conclusions

From the results presented in this study, the following conclusions can be drawn:
(1)
From 1978 to 2020, the value of capital elasticity of rice yield growth in China is between 0.3209 and 0.3589, with mean 0.3437, and the value of labor elasticity is between −0.1759 to −0.1640 with mean −0.1730, indicating capital elasticity and labor elasticity are not constant in different years.
(2)
The correlation coefficient between capital elasticity and the annual proportion of crop disasters is −0.9817, and correlation coefficient between labor elasticity (absolute value) and the annual proportion of crop disasters is −0.8752, presenting a negative relationship for both. With an increase in the annual proportion of crop disasters, the decreasing speed of capital elasticity and labor elasticity tends to increase. When the annual proportion of crop disasters is more than about 18%, the capital elasticity and the labor elasticity show a significant decline curve, proving that natural disasters have a great impact on capital elasticity and labor elasticity. Therefore, we should focus on the long-term strengthening of disaster prevention and resistance. In view of the possibility of the more frequent occurrence of extreme weather events, we should systematically consider various links such as variety cultivation, water conservancy projects, farmland construction, agricultural machinery and agronomy, explore scientific methods, and promote the improvement of the agricultural disaster prevention and relief work system and mechanism.
(3)
Compared with 1978, the GTPR of China’s rice yield growth from 1979 to 2020 shows a declining trend in fluctuations, whereas the TKR shows a rising trend in fluctuations and the TLR is relatively stable in the same period. The value of GTPR is between 15.30% and 58.71%, the TKR is between 9.57% and 65.52%, and the TLR is between 0.92% and 69.373%. Since 2000, the increase in rice yield per unit area in China mainly depended on the increase of capital investment. The top four input factors are machinery, chemical fertilizer, seed and pesticide in the order of contribution rate. To ensure the sustainable and healthy development of China’s rice industry, it is suggested to reasonably limit the use of chemical fertilizers and pesticides, vigorously develop agricultural mechanization technology, and strengthen rice variety cultivation.

Author Contributions

Z.L. and X.M. designed the study. Z.L., Conceptualization, methodology, formal analysis, investigation, resources, writing-original draft preparation, funding acquisition; X.W. (Xiaola Wu), software, validation, data curation, data editing; X.W. (Xicheng Wang), Writing-review and editing, supervision; H.Z. and J.C., Resources, data editing; X.M., Writing-review and editing, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China [grant No. 52175226], Laboratory of Lingnan Model Agriculture Project [grant No. NT2021009], Department of Science and Technology of Guangdong Province [grant No. KTP20210196], and the Earmarked Fund for Modern Agro-industry Technology Research System [grant No. CARS-01-47].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting the reported results can be found at China Statistical Yearbook and China Agricultural Product Cost-Benefit Compilation.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The values of output elasticity of capital and labor from 1978 to 2020.
Table A1. The values of output elasticity of capital and labor from 1978 to 2020.
Yearαβ
19780.3426−0.1735
19790.3504−0.1757
19800.3280−0.1671
19810.3466−0.1748
19820.3497−0.1756
19830.3495−0.1756
19840.3502−0.1757
19850.3408−0.1728
19860.3394−0.1722
19870.3444−0.1741
19880.3382−0.1717
19890.3388−0.1719
19900.3482−0.1753
19910.3337−0.1697
19920.3370−0.1712
19930.3411−0.1729
19940.3246−0.1656
19950.3428−0.1735
19960.3446−0.1742
19970.3302−0.1681
19980.3399−0.1724
19990.3378−0.1715
20000.3209−0.1640
20010.3277−0.1670
20020.3366−0.1710
20030.3239−0.1653
20040.3505−0.1757
20050.3467−0.1749
20060.3399−0.1724
20070.3395−0.1722
20080.3439−0.1739
20090.3455−0.1745
20100.3488−0.1754
20110.3546−0.1758
20120.3555−0.1757
20130.3532−0.1759
20140.3547−0.1758
20150.3550−0.1758
20160.3539−0.1759
20170.3564−0.1755
20180.3566−0.1755
20190.3589−0.1751
20200.3589−0.1751
Table A2. Input factor contribution rates in rice yield growth from 1979 to 2020.
Table A2. Input factor contribution rates in rice yield growth from 1979 to 2020.
YearGTPR/%NKR/%KDR/%TKR/%NLR/%LDR/%TLR/%
197942.8915.0641.1356.1912.77−11.850.92
198020.7040.82−31.259.5740.7528.9869.73
198132.3823.8516.1540.0032.93−5.3127.63
198258.0812.6811.1123.7921.20−3.0718.13
198356.7717.928.9726.8918.78−2.4416.34
198458.7117.517.9825.4917.82−2.0115.80
198558.4322.51−2.1120.4120.490.6821.17
198654.1128.58−3.5125.0719.681.1320.82
198740.6638.441.9340.3719.53−0.5618.98
198832.6649.63−5.2044.4321.271.6422.91
198938.2646.19−3.9242.2818.251.2119.46
199036.5643.085.0748.1416.58−1.2815.30
199141.6746.05−8.5237.5318.052.7520.80
199237.6148.73−4.8743.8717.041.4818.52
199329.9454.21−1.2852.9416.770.3617.13
199434.7258.95−15.8043.1517.204.9322.13
199531.5552.790.1752.9615.54−0.0515.49
199636.3748.231.5049.7314.29−0.3913.89
199746.7944.78−8.5836.2014.402.6117.01
199842.5143.33−1.7941.5415.480.4815.96
199939.6046.19−3.3342.8616.630.9017.53
200052.5841.04−15.2425.8017.194.4321.62
200145.0344.16−10.9933.1718.683.1221.81
200236.0647.88−4.4343.4519.321.1620.49
200341.5348.63−14.5234.1120.334.0324.36
200431.9643.755.4649.2119.79−0.9618.83
200533.5843.482.9246.4020.63−0.6120.02
200631.6747.98−1.9546.0321.860.4422.30
200737.2043.97−2.0641.9020.440.4520.89
200834.6844.710.8245.5319.95−0.1619.79
200933.6444.581.8246.4020.30−0.3319.96
201025.8749.773.9853.7421.02−0.6320.39
201123.0750.147.3657.5020.14−0.7119.43
201218.6854.077.8861.9620.01−0.6519.36
201317.3256.386.5362.9120.44−0.6819.76
201417.5255.827.2363.0520.06−0.6219.43
201517.8055.867.2763.1319.65−0.5919.06
201616.3257.436.7164.1520.12−0.5919.53
201715.7456.808.0064.8019.95−0.4819.47
201818.2155.027.8362.8419.38−0.4418.94
201916.7655.309.0464.3419.25−0.3418.91
202015.3056.429.1065.5219.51−0.3319.18
Table A3. Capital contribution rate and its decomposition from 2000 to 2020.
Table A3. Capital contribution rate and its decomposition from 2000 to 2020.
YearSeed Cost/%Chemical Fertilizer Cost/%Farm Fertilizer Cost/%Pesticide Cost/%Agricultural Film Cost/%Mechanical Operation
Cost /%
Irrigation and Drainage Cost/%Animal Power Cost/%Fuel Power Cost/%Other Cost/%
20001.987.380.911.930.362.762.332.380.015.77
20012.269.391.172.600.433.643.212.990.007.47
20023.2012.001.453.260.544.903.883.510.0110.70
20032.489.651.002.870.384.043.162.730.087.71
20043.6215.831.794.830.667.113.993.860.217.29
20053.8916.291.575.030.707.953.513.260.104.10
20063.9315.181.366.060.679.283.572.990.122.88
20073.5213.590.995.670.599.403.032.750.072.30
20083.5016.541.065.650.5610.912.622.560.082.04
20094.1015.041.105.660.4912.192.812.360.322.34
20105.4215.881.346.470.5815.722.972.250.222.89
20115.9717.441.116.250.5817.563.152.090.333.01
20126.6018.251.156.690.6020.103.041.850.343.33
20136.9217.561.106.630.5921.463.201.440.393.60
20147.2816.221.046.740.6022.892.771.360.443.72
20157.3016.071.136.750.5923.172.731.180.453.77
20167.6115.881.136.790.6223.932.740.990.513.95
20177.9616.041.166.900.5924.032.640.750.654.07
20187.7415.991.176.540.5423.302.530.530.573.93
20197.8816.611.146.860.5123.732.740.410.653.81
20208.1816.471.177.350.5124.242.600.320.723.97

References

  1. Mohidem, N.A.; Hashim, N.; Shamsudin, R.; Che Man, H. Rice for food security: Revisiting its production, diversity, rice milling process and nutrient content. Agriculture 2022, 12, 741. [Google Scholar] [CrossRef]
  2. FAOSTAT. Production/Yield Quantities of Rice, Paddy in World + (Total) 1978–2020. Food and Agriculture Organization of the United Nations. 2022. Available online: https://www.fao.org/faostat/en/#data/QCL/visualize (accessed on 2 September 2022).
  3. Huang, M.; Shan, S.L.; Xie, X.B.; Cao, F.B.; Zou, Y.B. Why high grain yield can be achieved in single seedling machine-transplanted hybrid rice under dense planting conditions? J. Integr. Agric. 2018, 17, 1299–1306. [Google Scholar] [CrossRef]
  4. Zhang, X.F.; Wang, D.Y.; Fang, F.P.; Zhen, Y.K.; Liao, X.Y. Food safety and rice production in China. Res. Agric. Mod. 2005, 26, 85–88. [Google Scholar]
  5. Ma, G.H.; Yuan, L.P. Hybrid rice achievements, development and prospect in China. J. Integr. Agric. 2015, 14, 197–205. [Google Scholar] [CrossRef]
  6. Li, J.J.; Li, J.P. International comparison of rice production cost and benefit and China’s development prospects. China Rice 2021, 27, 22–30. [Google Scholar]
  7. Huang, M.; Tang, Q.Y.; Ao, H.J.; Zou, Y.B. Yield potential and stability in super hybrid rice and its production strategies. J. Integr. Agric. 2017, 16, 1009–1017. [Google Scholar] [CrossRef]
  8. Xu, L.; Yuan, S.; Man, J.G. Changes in rice yield and yield stability in China during the past six decades. J. Sci. Food Agric. 2020, 100, 3560–3569. [Google Scholar] [CrossRef] [PubMed]
  9. Li, Z.H.; Ma, X.; Chen, L.T.; Li, H.W.; Huang, Y.Q.; Li, J.R.; Luo, G.W.; Yao, J.H. Effects of coupling of nursing seedling densities and seedling fetching area on transplanting quality and yield of hybrid rice. Trans. Chin. Soc. Agric. Eng. (Trans. CSAE) 2019, 35, 20–30. [Google Scholar]
  10. Espe, M.B.; Hill, J.E.; Leinfelder-Miles, M.; Espino, L.A.; Mutters, R.; Mackill, D.; Kessel, C.V.; Linquist, B.A. Rice yield improvements through plant breeding are offset by inherent yield declines over time. Field Crops Res. 2018, 222, 59–65. [Google Scholar] [CrossRef]
  11. Umetsu, C.; Lekprichakul, T.; Chakravorty, U. Efficiency and technical change in the Philippine rice sector: A Malmquist total factor productivity analysis. Am. J. Agric. Econ. 2003, 85, 943–963. [Google Scholar] [CrossRef]
  12. Elahi, E.; Cui, W.J.; Zhang, H.M.; Nazeer, M. Agricultural intensification and damages to human health in relation to agrochemicals: Application of artificial intelligence. Land Use Policy 2019, 83, 461–474. [Google Scholar] [CrossRef]
  13. Liu, S.W.; Zhang, P.Y.; He, X.L.; Wang, Z.Y.; Tan, J.T. Productivity and efficiency change in China’s grain production during the new farm subsidies years: Evidence from the rice production. Custos E Agronegocio 2015, 11, 106–123. [Google Scholar]
  14. Mariano, M.J.; Villano, R.; Fleming, E. Are irrigated farming ecosystems more productive than rainfed farming systems in rice production in the Pilippines? Agric. Ecosyst. Environ. 2010, 139, 603–610. [Google Scholar] [CrossRef]
  15. Mariyono, J. Decomposition total factor productivity of Indonesian rice production. Econ. J. Emerg. Mark. 2018, 10, 121–127. [Google Scholar] [CrossRef]
  16. Wang, J.Q.; Song, H.H.; Tian, Z.W.; Bei, J.L.; Zhang, H.Y.; Ye, B.; Ni, J. A method for estimating output elasticity of input factors in Cobb-Douglas production function and measuring agricultural technological progress. IEEE Access 2021, 9, 26234–26250. [Google Scholar] [CrossRef]
  17. Ball, E.V.; Wang, S.L.; Nehring, R.; Mosheim, R. Productivity and economic growth in U.S. agriculture: A new look. Appl. Econ. Perspect. Policy 2016, 38, 30–49. [Google Scholar] [CrossRef]
  18. Li, Z.; Zhang, H.P. Productivity growth in China’s agriculture during 1985–2010. J. Integr. Agric. 2013, 12, 1896–1904. [Google Scholar] [CrossRef]
  19. Zhang, S.F.; Xu, B. Production functions with time-varying elasticities and under the catch-up consensus: Total factor productivity. China Econ. Q. 2009, 8, 551–568. [Google Scholar]
  20. Solow, R.W. Technical change and the aggregate production function. Rev. Econ. Stat. 1957, 39, 312–320. [Google Scholar] [CrossRef]
  21. Shumway, R.H. Applied Statistical Time Series Analysis; Prentice Hall: Englewood Cliffs, NJ, USA, 1988. [Google Scholar]
  22. Cleveland, W.S.; Grosse, E.; Shyu, W.M. Local regression models. In Statistical Models in S; Chambers, J.M., Hastie, T.J., Eds.; Wadsworth & Brooks: Pacific Grove, CA, USA, 1991; pp. 309–376. [Google Scholar]
  23. Hastie, T.J.; Tibshirani, R.J. Varying-coefficient models. J. R. Stat. Soc. Ser. B 1993, 55, 757–796. [Google Scholar] [CrossRef]
  24. Cai, Z.W.; Fan, J.Q.; Li, R.Z. Efficient estimation and inferences for varying-coefficient models. J. Am. Stat. Assoc. 2000, 95, 888–902. [Google Scholar] [CrossRef]
  25. Ahmad, I.; Leelahanon, S.; Li, Q. Efficient estimation of a semiparametric partially linear varying coefficient model. Ann. Stat. 2005, 33, 258–283. [Google Scholar] [CrossRef]
  26. Luo, X.H.; Yang, Z.H.; Zhou, Y. Nonparametric estimation of the production function with time-varying elasticity coefficients. Syst. Eng. Theory Pract. 2009, 29, 144–149. [Google Scholar]
  27. Zhang, S.F.; Xu, B.; Gu, W.T. Statistical and economic tests for the time-varying elasticity production function model. Stat. Res. 2011, 28, 91–96. [Google Scholar]
  28. Zhang, S.F.; Zhu, J.N.; Fang, Q.; Liu, Y.X.; Xu, S.W.; Tsai, M.H. Solving the time-varying Cobb-Douglas production function using a varying-coefficient quantile regression model. J. Adv. Comput. Intell. Intell. Inform. 2019, 23, 831–837. [Google Scholar] [CrossRef]
  29. Yang, J.; Xiong, W.; Yang, X.G.; Cao, Y.; Feng, L.Z. Geographic variation of rice yield response to past climate change in China. J. Integr. Agric. 2014, 13, 1586–1598. [Google Scholar] [CrossRef]
  30. Song, C.X.; Liu, R.F.; Oxley, L.; Ma, H.Y. Do farmers care about climate change? Evidence from five major grain producing areas of China. J. Integr. Agric. 2019, 18, 1402–1414. [Google Scholar] [CrossRef]
  31. Zhang, R.Q.; Li, G.Y. Averaged estimation of functional-coefficient regression models with different smoothing variables. Stat. Probab. Lett. 2007, 77, 455–461. [Google Scholar] [CrossRef]
  32. Zhang, R.Q.; Feng, J.Y. Varying-coefficient regression models with different smoothing variables. Acta Math. Appl. Sin. 2007, 30, 444–451. [Google Scholar]
  33. Fan, G.L.; Xu, H.X. Averaged and integrated estimations of varying-coefficient regression models with dependent observations. Discret. Dyn. Nat. Soc. Hindawi 2019, 2019, 7146793. [Google Scholar] [CrossRef]
  34. Lee, Y.K.; Mammen, E.; Park, B.U. Backfitting and smooth backfitting in varying coefficient quantile regression. Econom 2014, 17, 20–38. [Google Scholar] [CrossRef]
  35. Yang, S.J. Estimation for semiparametric varying coefficient models with different smoothing variables under random right censoring. J. Korean Stat. Soc. 2018, 47, 161–171. [Google Scholar] [CrossRef]
  36. Tian, Y.; Li, B.; Zhang, J.B. Measurement of contribution rate of rice technological progress in China. Stat. Decis. 2012, 350, 93–95. [Google Scholar]
  37. Li, Z.H.; Ma, X. Measurement and dynamic characteristics of contribution rate of input factors of rice yield growth in China. Stat. Decis. 2016, 447, 94–97. [Google Scholar]
  38. Chen, Y.; Fu, Y.; Zhou, W.; Sun, H.; Wang, T.; Ren, W.J. The change of rice input-output characteristic and influence of technical progress on rice yield in Sichuan province. Chin. J. Agric. Resour. Reg. Plan. 2018, 39, 24–31. [Google Scholar]
  39. Wan, Z.; Zhang, C.; Fang, W. Calculation and secondary separation on contribution rate of scientific and technological progress of rice industry in Guangdong province. Technol. Econ. 2012, 31, 63–67. [Google Scholar]
  40. Carvalho, F.P. Pesticides, environment, and food safety. Food Energy Secur. 2017, 6, 48–60. [Google Scholar] [CrossRef]
  41. Arunrat, N.; Sereenonchai, S.; Chaowiwat, W.; Wang, C.; Hatanno, R. Carbon, nitrogen and water footprints of organic rice and conventional rice over 4 years of cultivation: A case study in the lower North of Thailand. Agronomy 2022, 12, 380. [Google Scholar] [CrossRef]
  42. Xu, C.C.; Ji, L.; Li, F.B.; Feng, J.F.; Fang, F.P. Situation and strategies of rice industry development in China. J. Huazhong Agric. Univ. 2022, 41, 21–27. [Google Scholar]
Figure 1. The fitting results of the model.
Figure 1. The fitting results of the model.
Agriculture 12 01431 g001
Figure 2. The output elasticity of capital from 1978 to 2020.
Figure 2. The output elasticity of capital from 1978 to 2020.
Agriculture 12 01431 g002
Figure 3. The output elasticity of labor from 1978 to 2020.
Figure 3. The output elasticity of labor from 1978 to 2020.
Agriculture 12 01431 g003
Figure 4. The relationship between the output elasticity of capital and the annual proportion of crop disasters.
Figure 4. The relationship between the output elasticity of capital and the annual proportion of crop disasters.
Agriculture 12 01431 g004
Figure 5. The relationship between the output elasticity of labor and the annual proportion of crop disasters.
Figure 5. The relationship between the output elasticity of labor and the annual proportion of crop disasters.
Agriculture 12 01431 g005
Figure 6. The contribution rates of rice input factors from 1979 to 2020.
Figure 6. The contribution rates of rice input factors from 1979 to 2020.
Agriculture 12 01431 g006
Figure 7. The trend of capital contribution rate and its decomposition from 2000 to 2020.
Figure 7. The trend of capital contribution rate and its decomposition from 2000 to 2020.
Agriculture 12 01431 g007
Table 1. Comparison of estimation results of two models.
Table 1. Comparison of estimation results of two models.
Data SetsC-D Production
Function Model
Varying-Coefficient Production Function Model
Capital ElasticityLabor ElasticityAverage Value of
Capital Elasticity
Average Value of
Labor Elasticity
1978–20160.5046−0.05180.3171−0.1747
1978–20170.5100−0.04620.3315−0.1738
1978–20180.5133−0.04310.3314−0.1737
1978–20190.5159−0.04060.3367−0.1732
1978–20200.5195−0.03740.3437−0.1730
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Z.; Wu, X.; Wang, X.; Zhong, H.; Chen, J.; Ma, X. Measurement and Analysis of Contribution Rate for China Rice Input Factors via a Varying-Coefficient Production Function Model. Agriculture 2022, 12, 1431. https://doi.org/10.3390/agriculture12091431

AMA Style

Li Z, Wu X, Wang X, Zhong H, Chen J, Ma X. Measurement and Analysis of Contribution Rate for China Rice Input Factors via a Varying-Coefficient Production Function Model. Agriculture. 2022; 12(9):1431. https://doi.org/10.3390/agriculture12091431

Chicago/Turabian Style

Li, Zehua, Xiaola Wu, Xicheng Wang, Haimin Zhong, Jiongtao Chen, and Xu Ma. 2022. "Measurement and Analysis of Contribution Rate for China Rice Input Factors via a Varying-Coefficient Production Function Model" Agriculture 12, no. 9: 1431. https://doi.org/10.3390/agriculture12091431

APA Style

Li, Z., Wu, X., Wang, X., Zhong, H., Chen, J., & Ma, X. (2022). Measurement and Analysis of Contribution Rate for China Rice Input Factors via a Varying-Coefficient Production Function Model. Agriculture, 12(9), 1431. https://doi.org/10.3390/agriculture12091431

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop