I. Introduction

The world has gone through a significant demographic upheaval in the last several decades. A growing working-age population is a crucial determinant of economic prosperity. However, since the world’s working-age population is shrinking while the non-working population is growing, there is the threat of an increase in the cost of care, thereby putting pressure on government-sponsored programmes like pension systems, medicare and other social security programs. According to United Nations demographics estimates, by the end of 2050, the old-age dependency ratio will be doubled. Although rising life expectancy reflects the success of the economy in medical advancements and in the public health sector, it also raises concerns about future economic prospects, the functioning of social security programs, and pension systems. Demographics do not directly dictate the course of economic growth, but act as a significant channel in determining an economy’s development potential as demographic shifts affect the economy’s demand and supply in both qualitative and quantitative ways.

The Life Cycle Hypothesis of consumption and saving (Ando & Modigliani, 1963) has impacted research into monetary policy mechanisms that incorporate the impact of age composition. Regarding the proponents of the Life Cycle Hypothesis posit that income varies over people’s lifetime, and saving allows them to move income from those times in life when income is high to the times when income is low. Thus, people tend to accumulate their savings towards their retirement when working, and then spend whatever has been saved when they retire and grow old.

In the same vein, the conclusion could be drawn that an economy’s national savings and age structure could be highly related to this. Thus, changes in demographics would have an impact on an economy by altering the balance between savings and investment, resulting in monetary variables such as inflation rate being adjusted (Canon et al., 2015), exchange rate (Kohli et al., 2020; Lee & Park, 2015; Rose et al., 2009) and interest rate (Kara & von Thadden, 2016; Padha et al., 2021; Park & Kim, 2012). The age distribution of a country’s population is shown to be a crucial factor in the connection between monetary expansion and price increases, according to the research conducted by Kopecky (2022). In addition, recent studies conducted by Cravino et al. (2020), Cravino et al. (2020) and Auclert et al. (2021) establish an important relationship between monetary policy instruments and demographics.

India is one the fastest growing economies in this post-modernism era of global flux. Given its population size and considering life cycle patterns in economic activities, we are keenly motivated to understand how changes in the age-composition of India affect the performance of its monetary policy, particularly as no investigation has been done in the Indian context in this regard. The main thrust of this paper is to understand the impact of monetary policy instruments on India’s demography. This research considers the repo rate as a monetary policy variable and uses a policy-restricted structural vector autoregressive model (SVAR) to determine how the working-age population influences the response of inflation and unemployment to a policy rate shock over time.

II. Materials and Methodology

This study considers the repo rate (Repo) as the monetary policy variable. Inflation rate (Inf), unemployment rate (Unra) and exchange rate (Exra) are taken as responsive channels in order to trace the impact of policy rate shock. The population variable (Pop) being considered is the population of working age cohorts representing the demographic structure of the economy collectively, as the working-age cohorts are the major stakeholders in economic activities of an economy. In this study, the repo rate is considered a shock in the economy to evaluate its real impact on the demographic. We use annual data on the age composition of India, given by working age population and aggregate unemployment rate.

Firstly, in order to check the oddity of the variables used in the analysis, this study runs a unit root test as time series data is used. Augmented Dickey Fuller (ADF) test is used to determine variables’ time series characteristics. Following that, to inspect how age composition affects the response of the inflation rate and unemployment rate to the monetary policy shocks, this study considers the impulse response function with one standard deviation shock in a SVAR environment.

A. Structural Vector-Autoregression

In order to understand the relationship between monetary policy instruments and India’s demographic structure, the following SVAR is formulated.

RYt=α0+α1Yt1+ut

The generalized form can be represented as

RYt=nj=1α1Yt1+ut,utN(0,K)

In the above generalized form, Yt(=Y1t, Y2t,Y3t,Ykt) is a k×1 vector of endogenous variables, α1αk is a k×k matrix of lag coefficients to be estimated on jth lag, ut is a k×1 vector of white noise innovation processes representing the structural shocks. R is a matrix that reflects the contemporaneous relationship among the endogenous variables.

Now considering the study , here vector Yt has five variables which are (Repo), (Inf), (Unra), (Exra) and (Pop) and the system is represented as shown below.

a11Repot+a12Inft+a13Unrat+a14Exrat+a15Popt=α10+α11Repot1+α12Inft1+α13Unrat1+a14Exrat1+α15Popt1+uRepot

a21Repot+a22Inft+a23Unrat+a24Exrat+a25Popt=α20+α21Repot1+α22Inft1+α23Unrat1+a24Exrat1+α25Popt1+uInft

a31Repot+a32Inft+a33Unrat+a34Exrat+a35Popt=α30+α31Repot1+α32Inft1+α33Unrat1+a34Exrat1+α35Popt1+uUnrat

a41Repot+a42Inft+a43Unrat+a44Exrat+a45Popt=α40+α41Repot1+α42Inft1+α43Unrat1+α44Exrat1+α45Popt1+uExrat

a51Repot+a52Inft+a53Unrat+a54Exrat+a55Popt=α50+α51Repot1+α52Inft1+α53Unrat1+α54Exrat1+α55Popt1+uPopt

The matrix form of the above system is given as

[1a12a13a14a15a211a23a24a25a31a321a34a35a41a42a431a45a51a52a53a541][RepotInftUnratExratPopt]=[α10α20α30α40α50][α11α12 α13α14α15α21α22α23α24α25α31α32α33α34α35α41α42α43α44α45α51α52α53α54α55]+[Repot1Inft1Unrat1Exrat1Popt1]+[uRepotuInftuUnratuExrat uPopt]

Now in order to derive reduced form VAR, multiply both sides of Equation (1) by R1 ; we get

R1RYt=R1α0+R1α1Yt1+R1ut

We know that multiplying a matrix with its inverse gives an identity matrix i.e., R1R=I, where I is the identity matrix. Thus, the above equation can be written as

Yt=Z0+Z1Yt1+et

where vector Yt depends on lag of itself and the forecast error et. The matrix R here also indicates the forecast errors of VAR, et and the structural shocks ut i.e., et=R1ut. In reduced form VAR, et is the linear makeup of ut.

Structural VAR is basically a conceptually constructed framework and is non-observable. Thus, it cannot be computed directly. As we have time series or history of variables, we use these to estimate the VAR. Following that, the regression of each variable is run against its lag and the lag of other variables included in the VAR. Hence, we get the coefficient Z and the forecast error et.

Starting from the reduced form VAR estimation, our goal is to get the structural model that isolates the purely exogenous shocks and responds to the variables of interest. To do this, we would need to get matrix R, or more precisely, we need to restrict matrix R.

Get R and multiply the reduced form VAR by R to get the structural model, shocks and contemporaneous channelization among the variables.

RYt=RZ0+RZ1Yt1+Ret

RYt=α0+α1Yt1+ut(

After getting the structural model, we need to impose restrictions on contemporaneous relationships among endogenous variables of the structural model. This what identification is about. Identification means imposing restrictions on matrix R based on economic intuition.

So, SVAR with imposing short-run restriction is given as

R = \begin{bmatrix} \mbox{1} & \mbox{0} & C(7) & \mbox{0} & C(13) \\ \mbox{0} & \mbox{1} & C(8) & C(11) & C(14) \\ C(1) & C(4) & \mbox{1} & \mbox{0} & C(15) \\ C(2) & C(5) & C(9) & \mbox{1} & C(16) \\ C(3) & C(6) & C(10) & C(12) & \mbox{1} \\ \end{bmatrix}

And SVAR with imposing long-run restriction is given as

R = \begin{bmatrix} \mbox{1} & \mbox{0} & \mbox{0} & C(8) & C(11) \\ \mbox{0} & \mbox{1} & C(6) & C(9) & C(12) \\ C(1) & C(4) & \mbox{1} & \mbox{0} & C(13) \\ C(2) & \mbox{0} & \mbox{0} & \mbox{1} & C(14) \\ C(3) & C(5) & C(7) & C(10) & \mbox{1} \\ \end{bmatrix}

III. Econometrics Results and Discussion

This study uses the Akaike Information Criterion (AIC), Schwarz’s Information Criterion (SIC), and Hannan-Quinn Information Criterion(HIC) to determine the optimal lag length for the VAR model. Considering these criteria the optimal lag length turns out to be 2. To examine the time series traits of variables considered in the analysis, we conduct unit root tests The ADF test shows that all variables, Repo, Inf, Unra, Exra, and Pop are stationary at I(1). Repo, Inf, Exra, and Pop are significant at the 5% level whereas Unra is significant at the 10% level.

Table 1.Results of unit root test
ADF (Level) ADF (First Difference)
Variables Constant
(P — Value)
Constant & Trend
(P — Value)
Constant
(P — Value)
Constant & Trend
(P — Value)
Repo 0.0087 0.0905 0.0004*** 0.0013***
Inf 0.1338 0.0110 0.000*** 0.0001***
Unmp 0.5990 0.2472 0.0255*** 0.0835**
InEx 0.8390 0.7461 0.0010*** 0.0062***
Popw 0.9570 0.1781 0.0003*** 0.0024***

Note: This table reports ADF unit root test results. *** and ** indicates statistical significance at the 1% and 5% levels, respectively.

Table 2.VAR lag length selection criterions.
Lags AIC SC HQ
0 0.432020 0.678866L 0.494101
1 0.275734 1.756814 0.648222
2 -1.121162*** 1.594150 -0.438269***

This table reports lag length criterion based results. *** indicates statistical significance at the 1% level.

A. Estimation of SVAR for Short Run Restrictions

We have discussed the matrix with short run restriction earlier in Section II. The matrix is represented as

\small{S = \left\lbrack \begin{matrix} . \\ Repo \\ Inf \\ Unra \\ Exra \\ Pop \\ \end{matrix}\ \begin{matrix} Repo & Inf & Unra & Exra & Pop \\ \mbox{1} & \mbox{0} & C(7) & \mbox{0} & C(13) \\ \mbox{0} & \mbox{1} & C(8) & C(11) & C(14) \\ C(1) & C(4) & \mbox{1} & \mbox{0} & C(15) \\ C(2) & C(5) & C(9) & \mbox{1} & C(16) \\ C(3) & C(6) & C(10) & C(12) & \mbox{1} \\ \end{matrix} \right\rbrack}

Now the estimated matrix with coefficient values is given as

\small{S = \ \begin{bmatrix} \mbox{1} & \mbox{0} & -6.339896 & \mbox{0} & -4.037263 \\ \mbox{0} & \mbox{1} & \mbox{16.24683} & \mbox{1.469397} & 10.56553 \\ -\mbox{0.101010} & 0.056936 & \mbox{1} & \mbox{0} & 0.539517 \\ -\mbox{0.698480} & \mbox{0.410654} & 7.148152 & \mbox{1} & 3.965317 \\ -\mbox{0.185562} & \mbox{0.126049} & \mbox{1.723454} & 0.259399 & \mbox{1} \\ \end{bmatrix}}

The coefficient values obtained are represented in Table 3.

Table 3.Structural VAR estimates
Coefficient Std. Error z- Statistic Probability
C(1) -0.101010 0.093654 -1.078545 0.2808
C(2) -0.698480 0.751231 -0.929780 0.3525
C(3) -0.185562 0.160352 -1.157214 0.2472
C(4) 0.056936 0.155809 0.365419 0.7148
C(5) 0.410654 1.093718 0.375466 0.7073
C(6) 0.126049 0.200473 0.628754 0.5295
C(7) -6.339896 2.451572 -2.586054 0.0097***
C(8) 16.24683 10.45665 1.553731 0.1202
C(9) 7.148152 1.949300 3.667036 0.0002***
C(10) 1.723454 0.518550 3.323604 0.0009***
C(11) 1.469397 1.796597 0.817878 0.4134
C(12) 0.259399 0.081132 3.197253 0.0014***
C(13) -4.037263 1.707367 -2.364614 0.0180***
C(14) 10.56553 6.366176 1.659635 0.0970**
C(15) 0.539517 0.161728 3.335950 0.0009***
C(16) 3.965317 1.212054 3.271567 0.0011***

Note: This table reports results obtained using the structural VAR model. *** and * indicates statistical significance at the 1% and 5% levels, respectively.

Considering the output from matrix S, we fInd that the coefficients’ values for the shock’s impact on Pop from Repo, Inf, Unra, and Exra are significant. With a shock in the repo rate (monetary policy control variable), its impact is significantly transferred to the working-age groups, which means that the working class is responsive to the shocks in the policy rate. The impact of the shock is highly negative on Pop in the short-run as shown by the coefficient value (-4.037263). This finding is similar to Cravino et al. (2020), Cravino et al. (2020) and Kopecky (2022).

B. Estimation of SVAR for Long-run Restrictions

The matrix F imposed with long-run restriction has also been discussed earlier in Section II , and the matrix is represented as

\small{F = \left\lbrack \begin{matrix} . \\ Repo \\ Inf \\ Unra \\ Exra \\ Pop \\ \end{matrix}\ \begin{matrix} Repo & Inf & Unra & Exra & Pop \\ \mbox{1} & \mbox{0} & \mbox{0} & C(8) & C(11) \\ \mbox{0} & \mbox{1} & C(6) & C(9) & C(12) \\ C(1) & C(4) & \mbox{1} & \mbox{0} & C(13) \\ C(2) & \mbox{0} & \mbox{0} & \mbox{1} & C(14) \\ C(3) & C(5) & C(7) & C(10) & \mbox{1} \\ \end{matrix} \right\rbrack}

The estimated matrix with coefficient values is then given as

\small{F = \begin{bmatrix} \mbox{1} & \mbox{0} & \mbox{0} & 2.937215 & -1.923777 \\ \mbox{0} & \mbox{1} & -8.111631 & - \mbox{13.64615} & 0.439337 \\ -0.300726 & -0.011206 & \mbox{1} & \mbox{0} & 0.865437 \\ 0.167861 & \mbox{0} & \mbox{0} & \mbox{1} & -0.524561 \\ -0.327495 & - \mbox{0.027185} & 0.772721 & -0.483892 & \mbox{1} \\ \end{bmatrix}}

The coefficient values obtained are represented in the Table 4.

Table 4.Structural VAR estimates
Coefficient Std. Error z — Statistic Probability
C(1) -0.300726 0.143442 -2.096497 0.0360***
C(2) 0.167861 0.072017 2.330836 0.0198***
C(3) -0.327495 0.118907 -2.754201 0.0059***
C(4) -0.011206 0.032607 -0.343683 0.7311
C(5) -0.027185 0.096274 -0.282372 0.7777
C(6) -8.111631 2.464636 -3.291209 0.0010***
C(7) 0.772721 0.125526 6.155858 0.0000***
C(8) 2.937215 0.976310 3.008487 0.0026***
C(9) -13.64615 2.302084 -5.927739 0.0000***
C(10) -0.483892 0.126812 -3.815813 0.0001***
C(11) -1.923777 1.110844 -1.731815 0.0833**
C(12) 0.439337 5.793302 0.075835 0.9396
C(13) 0.865437 0.215622 4.013672 0.0001***
C(14) -0.524561 0.268953 -1.950377 0.0511**

Note: *** and * indicates statistical significance at the 1% and 5% levels, respectively.

Considering the output of matrix F, we find that the impact of the shock on Pop from Repo, Unra, and Exra turn out to be significant, whereas in the case of Inf it turns out to be insignificant. Also, the impact of demographic shocks on Repo turn out to be significant. So, based on that, it can be said that matrix F has long-run restriction. Over time, the impact of shock in the repo rate has significant effects on the working-age population. As witnessed from the matrix, the impact of shock on Pop is negative, represented by the coefficient value (-1.923777). Even in this case, Pop is negatively affected by shock, although to a lesser extent than that caused by short-run constraints. These findings are consistent with the findings of Auclert et al. (2021) and Stephen et al. (2020).

IV. Conclusion

Our study revolves around tracing the direction of the interaction between demographic changes and monetary policy tools. We conclude that the shock in policy rates has a significant impact on the working-age population variable (which, for the purposes of this study, represents the demographic structure), and that this impact has a negative effect on the working-age population variable. This conclusion can be drawn because the shocks in policy rate are from short to medium term time horizon.

In placing the long-run restrictions, the impact of shocks on the working-age population is also negative. This study is only a step towards understanding the interaction between monetary policy and the demographics of the economy, in this case, for India.

Based on the relationship between monetary policy and demographic structure (represented by the working-age population), which was drawn through empirical analysis, the government should keep rolling out new incentive schemes from time to time in order to keep working cohorts population of different sectors motivated and above a certain benchmark. This study suggests that the working-age population is the key economic growth and prosperity driver. Thus, our study demands rich demographic data for deeper insights, which can lead us to better understand the macroeconomics of demography.