EconBase
← Back to paper

An econometric analysis of the Italian cultural supply

Extracted main text — title through conclusion, appendix excluded. This is what our citation measures are computed over, published so the extraction can be checked by eye.

56,840 characters · 6 sections · 37 citation commands

Rendered from LaTeX for readability, not typeset faithfully. Citation keys are highlighted; maths is left as source; figures, tables and equation environments are summarised rather than reproduced; unrecognised commands are greyed out so nothing is silently dropped. Email addresses are removed.

\doublespacing

center[center omitted — 99 chars of source]
center[center omitted — 368 chars of source]
quote\begin{center} {\sl Abstract} \end{center} Price indexes in time and space is a most relevant topic in statistical analysis from both the methodological and the application side. In this paper a price index providing a novel and effective solution to price indexes over several periods and among several countries, that is in both a multi-period and a multilateral framework, is devised. The reference basket of the devised index is the union of the intersections of the baskets of all periods/countries in pairs. As such, it provides a broader coverage than usual indexes. Index closed-form expressions and updating formulas are provided and properties investigated. Last, applications with real and simulated data provide evidence of the performance of the index at stake. {\sl Keywords:} multi-period price index, multilateral price index, stochastic approach, updating process, OLS. {\sl JEL code: C43; E31; C01.}

Introduction

Multi-period and multilateral price indexes, used to compare sets of commodities over time and across countries, respectively, are of prominent interest for statisticians Biggeri2010. Several approaches to the problem have been carried out in the literature.

One of them is the axiomatic approach balk1995axiomatic, which rests on the availability of both quantities and prices, dealt with as independent random variables, and aims at obtaining price indexes, to enjoy suitable properties fisher1921best, fisher1922making.

A second approach hinges on the economic theory\footnote{This approach is also known as preference field approach or functional approach divisia1926indice.} diewert1979economic, caves1982economic and rests on the idea that consumption choices come from the optimization of a utility function under budget constraints. Here, prices play the role of independent variables, while quantities arise as solutions to an optimization problem in accordance with the decision maker's preference scheme.

A third approach is the stochastic one clements2006stochastic, diewert2010stochastic, which can be traced back to the works of jevons1863serious, jevons1869depreciation and edgeworth1887report, edgeworth1925memorandum. Thanks to balk1980method, clements1987measurement, this approach has recently been reappraised, and its role in inflation measurements duly acknowledged zahid2010measuring. In this framework, prices are assumed to be affected by measurement errors whose bias effect must be duly minimized. The stochastic approach (hereafter, SA) turns out to be somewhat different from other approaches, insofar as it is closely related to regression theory theil1960best, clements1987measurement. In fact, the SA enables the construction of tests and confidence intervals for price indexes, which provide useful pieces of information clements2006stochastic. Furthermore, the SA has less limits than other approaches\footnote{The SA, differently from the index number theory originating from theil1967economics does not need to account for the economic importance of single prices.} and clears the way to further extensions, as shown in diewert2004stochastic, diewert2005weighted, silver2009hedonic, rao2004country.

In this paper, we devise a multi-period/multilateral price index, MPL index henceforth, within the stochastic framework. The derivation of the MPL index, which is the solution to an optimization problem, calls for quantities and values of the commodities (not prices), like walsh1901measurement. The reference basket, namely the set of commodities for all periods/countries, is made up of the union of the intersections of all the couples of year/country baskets in pairs. Namely, the price index of a commodity can be always computed once the latter is present in at least two periods/countries. Thus, the reference basket turns out to be more representative than the ones commonly used by the majority of statistical agencies, which either align the reference basket to that of the first period, or make it tally with the intersection of the commodity sets of all periods/countries. Eventually, such a reference basket is likely to be scarcely representative of the commodities present in each period/country. In this sense, just like hedonic pakes2003reconsideration, GESKS balk2012price and country/time-product-dummy (CPD/TPD) approaches with incomplete price tableaus rao2004country, the MPL index does not drop any observation on the account of having no counterpart in the reference basket. The lack of a commodity in a period/country $t$ requires putting both its quantity and value equal to zero in that period/country. However, unlike the aforesaid approaches, the MPL index is built on quantities and values, not on prices. Neither any preliminary computation of binary price indexes, as in the GESKS approach ivancic2011scanner, nor the use of any type of weighting matrix for dealing with missing values or quantities, as in the case of CPD/TPD indexes, are needed.

The updating of the MPL index is easy to accomplish and suitable formulas, tailored to the multi-period or multilateral nature of the data, are provided. In fact, while the inclusion of fresh values and quantities, of a set of commodities corresponding to an extra period, does not affect the previous values of the MPL index, the inclusion of a new country affects all former MPL indexes. Hence, two updating formulas have been proposed for the MPL index: one for the multi-period case and another for the multilateral case. A comparison of the said index to CPD/TPD index -- a multilateral/multi-period index that, like the MPL one, can be read as a solution to an optimization problem -- provides evidence of an easier implementation and greater accuracy of the former index.

To sum up, a threefold novelty characterizes the paper. First, a price index, which proves effective either for the multi-period or the multilateral case, is devised. Second, updating formulas tailored to the multilateral and the multi-period version of the index are provided. Third, the grater simplicity of use and accuracy of the said index is highlighted in comparison with well-known standard multilateral/multi-period indexes.

The paper is organized as follows. In Section (ref), we briefly go over the SA to price indexes and point out how several indexes are solutions to optimization problems. In Section (ref), within the SA, we devise the MPL index according to a minimum-norm criterion as well as its updating formulas for the multi-period and the multilateral cases, respectively. Section (ref) is devoted to the properties of the MPL index. Section (ref) provides an application of the MPL index to the Italian cultural supply data to shed light on its potential as both a multi-period and a multilateral index. To gain a better insight into the performance of the MPL index, a comparison is made of the CPD/TPD indexes by using both real and simulated data. Section (ref) completes the paper with some concluding remarks and hints. For the sake of easier readability, an Appendix has been added with proofs and technicalities.

The stochastic approach revisited

In this section, we review the main features of the SA and show how several well-known price indexes can be obtained within this framework.

The SA works out price indexes as solutions to an optimization problem consisting in finding a line (more generally a plane or a hyperplane) which lies as closely as possible to the points whose coordinates are the $N$ commodity prices in the $T$ periods under examination. Following theil1967economics and selvanathan1994index and by assuming for exposition purposes $T=2$, the idea underlying the SA is that, in both the periods taken into account, all prices move almost proportionally. Namely,

equation[equation omitted — 94 chars of source]

where $p_{it}$ is the price of commodity $i$ in period $t$ ($t=\{1,2\}$) and $\lambda$ is a scalar factor acting as price index. Eq. (ref) can be conveniently rewritten as follows

equation[equation omitted — 103 chars of source]

where $\eta_i$ are error terms which, as a rule, are assumed to be non-systematic and uncorrelated between commodities with variances, that can be either constant (with respect to commodities), that is

equation[equation omitted — 89 chars of source]

or not. In the latter case, variances are frequently specified as inversely proportional to the commodity budget share, namely

equation[equation omitted — 104 chars of source]

Here $w_{i1}=\frac{p_{i1}q_{i1}}{\boldsymbol{p}_1'\boldsymbol{q}_1}$ and $\boldsymbol{p}_1'\boldsymbol{q}_1$ is the total expenditure in the base period ($t=1$), being $\boldsymbol{p}_t$ and $\boldsymbol{q}_t$ the vectors of prices and quantities at time $t$. Eq. (ref) can be written in compact notation as follows

equation[equation omitted — 109 chars of source]

where $\boldsymbol{u}$ is an $N$-dimensional unit vector, $\boldsymbol{p}_1^{-H}$ is defined as the vector of the reciprocals of the non-null entries of $\boldsymbol{p}_1$ and zeros otherwise, that is

equation[equation omitted — 218 chars of source]

and $\boldsymbol{\eta}$ is a vector of $N$ non-systematic and spherical random variables, unless otherwise specified. Taking the Hadamard product of both sides of Eq. (ref) by $\boldsymbol{p}_1$ , yields the following

equation[equation omitted — 1,023 chars of source]

The appeal of the SA lays in the possibility of evaluating, besides point estimates, also price index standard errors which increase as the relative price variability increases. The computation of price index standard errors allow to verify the intuitive notion that the less prices move proportionally, the less precise are price index estimators. Further, standard errors prove useful to build confidence intervals for price indexes.

The following theorem shows that several well-known price indexes can be seen as offspring of Eq. (ref).

theoremLet $\boldsymbol{p}_t$ be a vector of the prices of $N$ commodities at time $t$ and $\boldsymbol{q}_t$ be the vector of the corresponding quantities. The Laspeyeres, Paasche, Marschall-Edgeworth and Walsh indexes are solutions to an optimization problem of the form \begin{equation} \min_{\lambda}\,\,\, || \boldsymbol{p}_2\,- \lambda \boldsymbol{p}_1||_N=\min_{\lambda}\,\,\, ||\boldsymbol{e}||_N \end{equation} where $||\,\cdot\,||_{N}=(\boldsymbol{e}'\boldsymbol{A}\boldsymbol{e})^{\frac{1}{2}}$ stands for a (semi)norm of the reference price index and $\boldsymbol{A}$ is a properly chosen non-negative definite matrix.
proofSee Appendix (ref).

The MPL index as solution to an optimization problem

In this section, a multi-period/multilateral price index is derived in the wake of the SA introduced in the previous section. The construction of this kind of index, MPL index hereafter, poses several issues, like the choice of the reference basket and its updating. When the prices of commodity sets in two periods/countries are compared, the reference basket, $K_{\tau}$, is generally set to be a subset of the commodities of the first period/country ($K_{\tau}=K_{1}$), which is also assumed to be representative for the other period/country.

The simplest solution is to take the “intersection” of the two baskets as reference basket, ($K_{\tau}= K_{1}\cap K_{2}$). When the comparison is among more than two periods/countries, statistical agencies generally align the reference basket to that of the first period/country ($K_{\tau}=K_{1}$), which turns out to play the role of base period/country. Of course, this choice is somewhat arbitrary as it leaves open the basket updating problem and it does not take into account links among baskets corresponding to couples, triples, $\dots$, of periods. An alternative approach sets the reference basket as the intersection of all the commodities considered in each single period ($K_{\tau}= \bigcap_{t} K_{t}$). As a result, the reference basket is likely to be partially representative of those commodities which are peculiar to each period.

Taking the SA as the reference frame, we derive a multi-period/multilateral price index, satisfying a minimum-norm criterion, whose reference basket -- over a set of periods or across a set of countries -- is the union of the intersections of the commodity baskets of various periods/countries, taken in pairs.\footnote{Some similarities arise with the chaining rule forsyth1981theory, von2001chain. This is a specific type of temporal aggregation method based on the use of the complete time series when computing a price index from $0$ to $t$. The resulting price index is a measure of the cumulated effect of adjacent periods from 0 to 1, 1 to 2, $\dots$, $t-1$ to $t$. However, chain indexes leave unresolved the reference basket updating and are not applicable in a multilateral perspective, differently from the MPL index. Chain indexes compare the current and the previous periods in order to evaluate the evolution over many periods. For a comparison of this approach with the fixed base one see diewert2001consumer} Such a reference basket proves to be an effective solution for several reasons. First, it is broader and more representative than the ones built on the intersection of commodities which are present in all periods/countries. Second, the price index of a commodity is well defined, provided the latter is present in at least two baskets. The computation of the said price index requires quantities and values of the commodities, not their prices. The lack of a commodity in a given period/country $t$ simply entails setting both its quantity and value equal to zero in that period/country. Figure (ref) shows the reference basket corresponding to the aforementioned approaches for the case of two and three periods/countries.

figure[figure omitted — 368 chars of source]

The MPL index hinges, first of all, on the idea that in each time/country $t$, the $N$ commodity prices move proportionally to a set of $N$ reference prices, that is

equation[equation omitted — 176 chars of source]

or

equation[equation omitted — 202 chars of source]

Here $\boldsymbol{p}_t$ is the actual price vector of the $N$ commodities at period/country $t$, $\tilde{\boldsymbol{p}}$ is the vector of (unknown (time invariant) reference prices, $\lambda_t$ is a scalar factor acting as price index at period/country $t$ and $\boldsymbol{\varsigma}_t$ is a vector of error terms. Eq. (ref) can be viewed as a generalization of Eq. (ref) for two periods/countries. As per Eq. (ref), in each period/country $t$, the $N$ prices ($\boldsymbol{p}_t$) can be represented by a point in a $N$-dimensional space. Accordingly, $N$ prices in $T$ periods/countries, $\boldsymbol{P}=[\boldsymbol{p}_{1},...,\boldsymbol{p}_{t},..,\boldsymbol{p}_{T}]$, can be represented by $T$ points in a $N$-dimensional space. If all prices move proportionally, these points would lie on a hyperplane, $\tilde{\boldsymbol{p}}$, and, in particular, on a straight line crossing the origin for $T=2$.

In general, this is only approximately true and the price “line” crossing the origin enjoys the property of fitting the observed price points, according to a criterion which minimizes the deviations of the data from the “line”. In compact notation, Eq. (ref) can be more conveniently reformulated as follows

equation[equation omitted — 203 chars of source]

where $\boldsymbol{\lambda}$ is the vector of the $T$ price indexes and $\tilde{\boldsymbol{p}}$ is the vector of the $N$ (unknown) reference prices. According to Eq. (ref), the problem of determining a set of $T$ price indexes can be read as the problem of approximating the price matrix, $\boldsymbol{P}$, with a matrix of unit rank, $\boldsymbol{\Pi}$, defined as the product of a vector of price indexes by a (unknown) reference price vector. Moving from prices to values, in light of Eq. (ref), the matrix $\boldsymbol{V}$ of the values of $N$ commodities in $T$ periods/countries has the representation

equation[equation omitted — 275 chars of source]

where $\boldsymbol{Q}$ is the matrix of the quantities of $N$ commodities in $T$ periods/countries. Eq. (ref) can be rewritten as

equation[equation omitted — 293 chars of source]

where $\boldsymbol{v}_{t}$ and $\boldsymbol{q}_{t}$ are the $t^{th}$ columns of $\boldsymbol{V}$ and $\boldsymbol{Q}$ respectively, and $\lambda_{t}$ is the $t^{th}$ element of $\lambda_{t}$. Eq. (ref) entails the following

equation[equation omitted — 145 chars of source]

where ${\boldsymbol{\varepsilon}_t}$ is a non-systematic stochastic term. The above formula, taking into account the identity \[\underset{(N,1)}{\tilde{\boldsymbol{p}}}*\underset{(N,1)}{\boldsymbol{q}_t}=\underset{(N,N)}{\boldsymbol{D}_{\tilde{\boldsymbol{p}}}} \underset{(N,1)}{\boldsymbol{q}_t}, \] can be re-written as

equation[equation omitted — 295 chars of source]

where $\delta_t=(\lambda_t)^{-1}$ takes the role of the deflator, and $\boldsymbol{D}_{\tilde{\boldsymbol{p}}}$ is a diagonal matrix with diagonal entries equal to the elements of $\tilde{\boldsymbol{p}}$ defined as in Eq. (ref).\footnote{The matrix $\boldsymbol{D}_{\tilde{\boldsymbol{p}}}$ is defined as follows: \[\boldsymbol{D}_{\tilde{\boldsymbol{p}}}=

bmatrix[bmatrix omitted — 142 chars of source]

. \]} Eq. (ref) expresses the value, $v_{it}$, of each commodity $i$ at time $t$ (discounted by a factor $\delta_t$) as the product between the (time invariant) reference price, $\tilde{p}_i$, and the corresponding quantity, $q_{it}$, plus an error term, $\varepsilon_{it}$. Note that $\lambda_t=\delta_t^{-1}$ plays the role of the multi-period/multilateral price index in Eq. (ref). Over $T$ periods/countries, the model can be written as

equation[equation omitted — 267 chars of source]

where $\boldsymbol{D}_{\boldsymbol{\delta}}$ is a $T\times T$ diagonal matrix with diagonal entries equal to the elements of $\boldsymbol{\delta}$ (see Footnote (ref)). Without lack of generality, we assume that the first period is the base period (that is $\delta_1$=$\lambda_1=1$), and write the first equation separately from the others $T-1$. Accordingly, the system in Eq. (ref) can be written as

equation[equation omitted — 531 chars of source]

or as

equation[equation omitted — 281 chars of source]

which, after some computations\footnote{Use has been made of the relationships \[vec(\boldsymbol{ABC})=(\boldsymbol{C}'\otimes \boldsymbol{A})\,vec(\boldsymbol{A}) \,\,\,\,\,\,\mbox{and}\,\,\, \,\,\, \underset{(N^2,1)}{vec(\boldsymbol{D}_{\boldsymbol{a}})}=\boldsymbol{R}'_N\,\underset{(N,1)}{\boldsymbol{a}}\] where $\boldsymbol{D}_{\boldsymbol{a}}$ is a diagonal matrix whose diagonal entries are the elements of the vector $\boldsymbol{a}$ and $\boldsymbol{R}_N$ is the transition matrix from the Kronecker to the Hadamard product faliva1996hadamard. }, can be also expressed as

equation[equation omitted — 427 chars of source]

Here, $\boldsymbol{\eta}= vec(\boldsymbol{E}_1)$, $\boldsymbol{\delta}$ is a vector whose elements are the diagonal entries of $\tilde{\boldsymbol{D}}_{\delta}$ as specified in Eq. (ref) and $\boldsymbol{R}_j$ denotes the transition matrix from the Kronecker to the Hadamard product.\footnote{The matrix $\boldsymbol{R}_{j}'$ is defined as follows \[\underset{(j^2\times j)}{\boldsymbol{R}_j'} =

bmatrix[bmatrix omitted — 266 chars of source]

, \] where $\boldsymbol{e}_i$ represents the $N$ dimensional $i$-th elementary vector. } The vector $\boldsymbol{\delta}$, whose (non-null) entries are the reciprocals of the elements of the price index $\boldsymbol{\lambda}$\footnote{The reciprocal of the deflator vector, $\boldsymbol{\delta}$, is defined in Eq. (ref).}, can be obtained by estimating Eq. (ref) with the ordinary least squares (OLS), under the assumption that

equation[equation omitted — 245 chars of source]

In this connection, we state the following result.

theoremThe MPL index can be obtained as reciprocal of the OLS estimate of the deflator vector, $\boldsymbol{\delta}$, in the system in Eq. (ref). This OLS estimate is given by \begin{equation}\begin{split} \underset{(T-1,1)}{\hat{{\boldsymbol{\delta}}}}=&\left\{ (\boldsymbol{I}_{T-1}\ast\boldsymbol{V}_1'\boldsymbol{V}_1)-(\boldsymbol{Q}_1'\ast\boldsymbol{V}_1')\left[ (\boldsymbol{q}_1\boldsymbol{q}_1'+\boldsymbol{Q}_1\boldsymbol{Q}_1')\ast\boldsymbol{I}_N \right]^{-1}(\boldsymbol{Q}_1\ast\boldsymbol{V}_1) \right\}^{-1}\&(\boldsymbol{Q}_1'\ast\boldsymbol{V}_1')\left[ (\boldsymbol{q}_1\boldsymbol{q}_1'+\boldsymbol{Q}_1\boldsymbol{Q}_1')\ast\boldsymbol{I}_N \right]^{-1}(\boldsymbol{q}_1\ast \boldsymbol{v}_1).\end{split}\end{equation}
proofSee Appendix (ref).

The case where $T=2$ is worth considering because it sheds light on the index structure. In this case, the price index turns out to be a ratio of weighted price averages, with weights depending on the harmonic means of the squared quantities as stated in the following corollary.

corollaryWhen $T=2$, the MPL index, $\widehat{\lambda}$, becomes \begin{equation} \hat{{\lambda}}=\frac{\sum_{i=1}^Np_{i2}\pi_i}{\sum_{i=1}^Np_{i1}\pi_i} \end{equation} where $\pi_i= 2 p_{i2}\frac{q_{i1}^2q_{i2}^2}{q_{i1}^2+q_{i2}^2}$ and $p_{it}=\frac{v_{it}}{q_{it}}$. The price index in Eq. (ref) can be also expressed as a convex linear combination of prices, that is \begin{equation} \hat{{\lambda}}=\sum_{i=1}^N\frac{p_{i2}}{p_{i1}}\tilde{\pi}_i \end{equation} with weights $\tilde{\pi}_{i}$ given by \begin{equation} \tilde{\pi}_{i}=\frac{v_{i1}v_{i2}\,\frac{q_{i1}q_{i2}}{q_{i1}^2+q_{i2}^2}}{\sum_{i=1}^N v_{i1}v_{i2}\,\frac{q_{i1}q_{i2}}{q_{i1}^2+q_{i2}^2}}. \end{equation} Furthermore, by setting $\tilde{\boldsymbol{q}}_t$=$\boldsymbol D^{-1/2}\boldsymbol q_t$ with $\boldsymbol D$ specified as follows \begin{equation} \boldsymbol D=\begin{bmatrix} q_{11}^2+q_{12}^2 & 0 & \dots &0 \\ 0 & q_{21}^2+q_{22}^2 & \dots &0 \\ \vdots & \vdots & \ddots &\vdots\\ 0 & 0 & \dots &q_{N1}^2+q_{N2}^2 \end{bmatrix}, \end{equation} the $MPL$ index can be also written in compact form as follows \begin{equation} \hat{\lambda}=\frac{\left(\tilde{\boldsymbol q}_1\ast \boldsymbol v_2\right)'\left(\tilde{\boldsymbol q}_1\ast \boldsymbol v_2\right)}{\left(\tilde{\boldsymbol q}_2\ast \boldsymbol v_2\right)'\left(\tilde{\boldsymbol q}_1\ast \boldsymbol v_1\right)}. \end{equation}
proofSee Appendix (ref).

As a by-product of Theorem (ref), we state the following result.

corollaryThe variance-covariance matrix of the deflator vector, $\boldsymbol{\widehat{\delta}}$, given in Theorem (ref) is \begin{equation} Var(\boldsymbol{\hat{\delta}})=\sigma^2[\boldsymbol{I}_{T-1}*\boldsymbol{V}_1'\boldsymbol{V}_1]^{-1}. \end{equation} The $t^{th}$ diagonal entry of the above matrix provides the variance of the deflator in the $t^{th}$period/country, given by \begin{equation} var({\hat{\delta}_{t}})=\sigma^2[\boldsymbol{v}_t'\boldsymbol{v}_t]^{-1} \end{equation} where $\boldsymbol{v}_t$ denotes the $t^{th}$ column of the matrix $\boldsymbol{V}_1$.
proofSee Appendix (ref).

As for the deflator vector $\boldsymbol{\hat{\delta}}$, its moments and confidence intervals can be easily obtained within the theory of linear regression models, from the result given in Corollary (ref). As the price index vector $\hat{\boldsymbol{\lambda}}$ turns out to be the reciprocal of the said deflator (see Eq. (ref)), its statistical behavior can be derived from the former, following the arguments put forward, for example, in geary1930frequency, curtiss1941distribution and marsaglia1965ratios, merely to quote a few, on ratios (in particular reciprocals) of random variables.

The following corollary provides an approximation of the variance of $\hat{\lambda}_t$, obtained by using the first Taylor expansion of the variance of a ratio of two random variables.

corollaryThe variance of the MPL index ${\hat{\lambda}_{t}}$ is \begin{equation} var({\hat{\lambda}_{t}}) \approx \frac{var({\hat{\delta}_{t}})}{E({\hat{\delta}_{t}}^4)} =\frac{\widehat{\sigma}^2}{[\boldsymbol{v}_t'\boldsymbol{v}_t]E({\hat{\delta}_{t}}^4)} \,\,\,\,\,\,\forall\,t=1,\dots,T \end{equation} In the above equation $\widehat{\sigma}^2=\frac{\widehat{\boldsymbol\mu}'\widehat{\boldsymbol \mu}}{NT-(N+T-1)}$ where $\widehat{\boldsymbol\mu}$ are the OLS residuals of the equation system (ref).
proofSee stuard1994kendall and elandt1980survival.

The following two theorems provide updating formulas for the price index $\boldsymbol{\hat{\lambda}}$. The former proves suitable when the index is used as a multilateral price index, while the latter is appropriate when it is employed as a multi-period index. In the former case, values and quantities of the commodities included in the reference basket are assumed available for an additional $T+1$ country. In the latter case, it is supposed that values and quantities of the commodities included in the reference basket become available at time $T+1$.

theoremShould the values and quantities of $N$ commodities of a reference basket become available for a new additional country, say the $T+1$-th, then, the updated multilateral version of the MPL index, $\boldsymbol{\hat{\lambda}}$, turns out to be the vector of the reciprocals, as defined in Eq. (ref), of the following deflator vector \begin{footnotesize} \begin{equation} \begin{split} \underset{(T,1)}{\hat{\boldsymbol{\delta}}}=&\\ =&\left\{ \begin{bmatrix} \boldsymbol{I}_{T-1}\ast\boldsymbol{V}'_1\boldsymbol{V}_1 & \boldsymbol{0} \\ \boldsymbol{0} & \boldsymbol{v}'_{T+1}\boldsymbol{v}_{T+1} \end{bmatrix}- \begin{bmatrix} \boldsymbol{Q}'_1\ast\boldsymbol{V}'_1 \\ \boldsymbol{v}'_{T+1}\ast\boldsymbol{q}'_{T+1} \end{bmatrix} \begin{bmatrix} \left(\boldsymbol{q}_1\boldsymbol{q}'_1+\boldsymbol{Q}_1\boldsymbol{Q}'_1 +\boldsymbol{q}_{T+1}\boldsymbol{q}'_{T+1}\right) \ast\boldsymbol{I}_N\end{bmatrix}^{-1} \right. \& \left.\begin{bmatrix} \boldsymbol{Q}_1\ast\boldsymbol{V}_1 & \boldsymbol{v}_{T+1}\ast\boldsymbol{q}_{T+1} \end{bmatrix} \right\}^{-1} \begin{bmatrix} \boldsymbol{Q}'_1\ast\boldsymbol{V}'_1 \\ \boldsymbol{v}'_{T+1}\ast\boldsymbol{q}'_{T+1} \end{bmatrix} \begin{bmatrix} \left(\boldsymbol{q}_1\boldsymbol{q}'_1+\boldsymbol{Q}_1\boldsymbol{Q}'_1 +\boldsymbol{q}_{T+1}\boldsymbol{q}'_{T+1}\right) \ast\boldsymbol{I}_N\end{bmatrix}^{-1} \left(\boldsymbol{q}_1\ast\boldsymbol{v}_1\right). \end{split} \end{equation} \end{footnotesize} Here the symbols are defined as in Theorem (ref) and $\boldsymbol{v}_{T+1}$, $\boldsymbol{q}_{T+1}$ denote the vector of values and quantities of $N$ commodities of the new $T+1$-th country, respectively.
proofSee Appendix (ref).
theoremShould the values and quantities of $N$ commodities of a reference basket become available for time $T+1$, then, the updated value $\hat{\lambda}_{T+1}$ of the multi-period version of the MPL index at time $T+1$ turns out to be the reciprocal of the deflator value at time $T+1$ \begin{equation} \begin{split} \underset{(1,1)}{\hat{\delta}_{T+1}}=&\left\{ \boldsymbol{v}_{T+1}'\boldsymbol{v}_{T+1}-\left(\boldsymbol{q}_{T+1}'\ast\boldsymbol{v}_{T+1}'\right)\left[ (\boldsymbol{Q}\boldsymbol{Q}'+\boldsymbol{q}_{T+1}\boldsymbol{q}_{T+1}')\ast\boldsymbol{I}_{N} \right]^{-1}\left(\boldsymbol{q}_{T+1}\ast\boldsymbol{v}_{T+1}\right) \right\}^{-1} \&\left(\boldsymbol{q}_{T+1}'\ast\boldsymbol{v}_{T+1}'\right)\left[ (\boldsymbol{Q}\boldsymbol{Q}'+\boldsymbol{q}_{T+1}\boldsymbol{q}_{T+1}')\ast\boldsymbol{I}_{N} \right]^{-1}\left(\boldsymbol{Q}\ast\boldsymbol{V}\right)\tilde{{\boldsymbol{\delta}}} \end{split} \end{equation} where $\tilde{{\boldsymbol{\delta}}}'=[1, \hat{{\boldsymbol{\delta}}}]'$ and $\hat{{\boldsymbol{\delta}}}$ is defined as in Eq. (ref).
proofSee Appendix (ref).

Figure (ref) highlights the difference between the updating process of the deflator, and thus of the price index, depending on whether it is used in the multilateral or in the multi-period case.

figure[figure omitted — 401 chars of source]

Properties of the MPL index

Let us assume for simplicity $T=2$ and denote with $\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ a generic index number where $\boldsymbol{p}_{t}$ and $\boldsymbol{q}_{t}$ are prices and quantities at time $t$. Without lack of generality, $t=1$ is assumed as the base period. Following predetti2006numeri and martini1992numeri, the main properties of an index number can be summarized as follows:

enumerate• Strong identity: $\hat{\lambda}(\boldsymbol{p}_{2},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=1$. • Commensurability: $\hat{\lambda}(\boldsymbol{\gamma}\ast\boldsymbol{p}_{1},\boldsymbol{\gamma}\ast\boldsymbol{p}_{2},\boldsymbol{\gamma}^{-1}\ast\boldsymbol{q}_{1},\boldsymbol{\gamma}^{-1}\ast\boldsymbol{q}_{2})=\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ with $\underset{(N,1)}{\boldsymbol{\gamma}}>\textbf{0}$ where $\boldsymbol{\gamma}^{-1}\ast \boldsymbol{q}_{t}=[q_{1t}/\gamma_1\,\,\dots\,\,q_{Nt}/\gamma_N]'$. • Proportionality: $\hat{\lambda}(\boldsymbol{p}_{1},\alpha\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\alpha\,\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ with $\alpha>0$. • Dimensionality: $\hat{\lambda}(\alpha\boldsymbol{p}_{1},\alpha\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ with $\alpha>0$. • Monotonicity: $\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{k}\ast\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})>\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ and $\hat{\lambda}(\boldsymbol{k}\ast\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})<\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ with $\underset{(N,1)}{\boldsymbol{k}}>\boldsymbol{u}$ where $\boldsymbol{u}$ is the unit vector.

We now prove the above properties for the MPL index in Eq. (ref), which will be denoted with the extended notation $\hat{{\lambda}}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$.

enumerate• Strong identity: taking the price vectors as equal in the two periods, it can be easily proved that \begin{equation} \hat{{\lambda}}(\boldsymbol{p}_{2},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^Np_{i2}\pi_i}{\sum_{i=1}^Np_{i2}\pi_i}=1\,\,\,where\,\,\,\pi_{i}=2p_{i2}\frac{q_{i1}^2q_{i2}^2}{q_{i1}^2+q_{i2}^2}. \end{equation} • Commensurability: this property is trivially verified upon noting that, when each element $p_{it}$ of the price vectors are multiplied by a positive constant $\gamma_{i}$ and the corresponding quantity $q_{it}$ is divided by $\gamma_i$, then the weights $\tilde{\pi}_i$ become $\tilde{\pi}_i=\frac{{\pi}_i}{\gamma_i}$. Accordingly, \begin{equation} \hat{{\lambda}}(\boldsymbol{\gamma}\ast\boldsymbol{p}_{1},\boldsymbol{\gamma}\ast \boldsymbol{p}_{2},\boldsymbol{\gamma}^{-1}\ast\boldsymbol{q}_{1},\boldsymbol{\gamma}^{-1}\ast\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N \gamma_i p_{i2}(\gamma_i ^{-1}\pi_i)}{\sum_{i=1}^N \gamma_i p_{i1} (\gamma_i^{-1}\pi_i)}=\frac{\sum_{i=1}^Np_{i2}\pi_i}{\sum_{i=1}^Np_{i1}\pi_i}. \end{equation} • Proportionality: this property is satisfied upon noting that $\breve{\pi}_i=\alpha{\pi}_i$. Accordingly, \begin{equation} \hat{{\lambda}}(\boldsymbol{p}_{1},\alpha\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N \alpha p_{i2}\breve{\pi}_i}{\sum_{i=1}^N p_{i1} \breve{\pi}_i}=\alpha\frac{\sum_{i=1}^N p_{i2}\pi_i}{\sum_{i=1}^N p_{i1} \pi_i}. \end{equation} • Dimensionality: the proof of this property follows the same argument used to prove homogeneity. Upon noting that $\breve{\pi_i}=\alpha{\pi}_i$, it follows that \begin{equation} \hat{{\lambda}}(\alpha\boldsymbol{p}_{1},\alpha\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N\alpha p_{i2}\breve{\pi}_i}{\sum_{i=1}^N\alpha p_{i1}\breve{\pi}_i}=\frac{\sum_{i=1}^N p_{i2}\pi_i}{\sum_{i=1}^Np_{i1}\pi_i}. \end{equation} • Monotonicity: if we consider $\boldsymbol{k}>\boldsymbol{u}$, then the associated price index is \[\hat{{\lambda}}(\boldsymbol{p}_{1},\boldsymbol{k}\ast\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N k_i p_{i2}\pi_i}{\sum_{i=1}^N p_{i1}\pi_i}>\frac{\sum_{i=1}^N p_{i2}\pi_i}{\sum_{i=1}^N p_{i1}\pi_i}=\hat{{\lambda}}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})\] if $k_i^2>k_j$ and $k_j^2>k_i$ for every $i\neq j$ given $i,j=1,\dots,N$. While if we consider $\boldsymbol{k}\ast\boldsymbol{p}_{1}$, taking into account that the weight vector $\boldsymbol{\pi}$ is independent from $\boldsymbol{k}$, the following holds \[\hat{{\lambda}}(\boldsymbol{k}\ast\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N p_{i2}\pi_i}{\sum_{i=1}^N k_i p_{i1}\pi_i}<\frac{\sum_{i=1}^N p_{i2}\pi_i}{\sum_{i=1}^N p_{i1}\pi_i}=\hat{{\lambda}}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2}).\]

It is worth noticing that the MPL index enjoys also the following properties:

itemize• Positivity: this property follows straightforward given that the index is the sum of ratios of non-negative quantities ($\pi_i, p_{i2}, p_{i1}\,\,\forall\,i=1,2,\dots, N$). • Inverse proportionality in the base period: let consider a vector of prices in the base period $\alpha \,\boldsymbol{p}_{1}$ with $\alpha>0$. Under this case, $\pi_i$ turns out to be independent from $\alpha$ and the associated index price proves to be proportional to $\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})$ \[\hat{\lambda}(\alpha\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N p_{i2} \pi_i}{\sum_{i=1}^N\alpha p_{i1} \pi_i}=\frac{1}{\alpha} \frac{\sum_{i=1}^N p_{i2} \pi_i}{\sum_{i=1}^Np_{i1} \pi_i}=\frac{1}{\alpha} \hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2}).\] • Commodity reversal property: it follows straightforward that the index price is invariant with respect to any permutation $(i)$: \[\hat{\lambda}(\boldsymbol{p}_{1},\boldsymbol{p}_{2},\boldsymbol{q}_{1},\boldsymbol{q}_{2})= \frac{\sum_{(i)=1}^{(N)} p_{(i)2} \pi_{(i)}}{\sum_{(i)=1}^{(N)}p_{(i)1} \pi_{(i)}}= \frac{\sum_{i=1}^N p_{i2} \pi_{i}}{\sum_{i=1}^Np_{i1} \pi_{i}}.\] • Quantity reversal test: a change in the quantity order only affects $\pi_i$ that remains invariant $\forall\, i=1,\dots,N$. Therefore the index price $\hat{\lambda}$ does not change. • Base reversibility (symmetric treatment of time): for one commodity, that is $N=1$, it is easy to prove that $\hat{\lambda}(p_{1},p_{2}, q_1,q_2)=\frac{p_2}{p_1}$ and, thus, $\hat{\lambda}(p_{2},p_{1}, q_2,q_1)=\frac{p_1}{p_2}=\hat{\lambda}(p_{1},p_{2}, q_1,q_2)^{-1}$. • Transitivity: for $N=1$, $\hat{\lambda}(p_{1},p_{2}, q_1,q_2)=\frac{p_2}{p_1}$, $\hat{\lambda}(p_{2},p_{3}, q_2,q_3)=\frac{p_3}{p_2}$, $\hat{\lambda}(p_{1},p_{3}, q_1,q_3)=\frac{p_3}{p_1}$ therefore, $\hat{\lambda}(p_{2},p_{3}, q_2,q_3)=\hat{\lambda}(p_{1},p_{2}, q_1,q_2) \,\cdot\,\hat{\lambda} (p_{2},p_{3}, q_2,q_3)$. • \textit{Monotonicity}: if $\boldsymbol{p}_{2}=\beta \boldsymbol{p}_{1}$ then \begin{equation} \tilde{\pi}_i=\beta\left(2p_{i1}\frac{q_{i1}^2q_{i2}^2}{q_{i1}^2+q_{i2}^2}\right)=\beta \gamma_i, \end{equation} and the associated index price turns out to be equal to $\beta$ \[\hat{\lambda}(\boldsymbol{p}_{1},\beta\boldsymbol{p}_{1},\boldsymbol{q}_{1},\boldsymbol{q}_{2})=\frac{\sum_{i=1}^N\beta^2\, p_{i1} \gamma_i}{\sum_{i=1}^N \beta p_{i1} \gamma_i}= \beta \frac{\sum_{i=1}^N p_{i1} \gamma_i}{\sum_{i=1}^N p_{i1} \gamma_i}= \beta.\]

An application of the MPL index to the Italian cultural supply data

In this section, we provide an application of the MPL index to Italian cultural supply data, such as revenues and the number of visitors to museums (i.e. monuments, archeological sites, museum circuits, $\dots$). The availability of temporal and geographical data on Italian culture provides a stimulating basis for ascertaining the potential of the MPL price-index methodology set forth in this paper. The flexibility of the MPL index paves the way to moving beyond ISTAT (and similar) analyses, which are confined to price indexes on the supply of data on Italian culture like access to museums and entertainment sectors, aggregated at the national level istatr2018b. In addition, to evaluate the performance of the MPL index, we have made a comparison with the CPD/TPD price indexes diewert2005weighted, rao2004country, using both real and simulated data. Reference has been made to this approach because, under the log-normality assumption of the error term, the maximum likelihood estimator of the said price index tallies with the least square one, likewise with the MPL index.

As for the nature of the data, note that Italian cultural heritage is at the top of various world-class lists and plays a key role in the Italian economy alderighi2018flight.\footnote{The Italian heritage supply chain accounts for 4,976 museums and the like; it generated almost 200 million euro of revenues in 2017 and employs more than 45,000 people istatr2016. } Lately, local cultural supply has evolved significantly. Indeed, most of the Italian museum circuits were founded relatively recently.\footnote{Approximately 2,300 sites (45.5%) of the Italian cultural supply chain were opened between 1960 and 1999, while 2,200 sites (38.6%) were opened in 2000, taking advantage of the investments for economic recovery and infrastructure enhancement made for Italian cultural heritage sites istatr2016.} In the following analysis, we have considered the ranking of the top 30 Italian cultural institutions (museums and the like) according to the highest number of annual visitors since 2004 (data source: \href{www.statistica.beniculturali.it}{www.statistica.beniculturali.it}). Among these, only 20 of the internationally renowned institutions remained ranked in the top 30 on a yearly basis. The other 10 positions were held by museums and institutions which only temporarily experienced an outstanding flow of visitors. These observations led to a twofold issue. First, the need to set-up a price index finalized at registering changes in the period under consideration. Second, a dynamic updating method of the index in order to preserve the information associated with the 10 positions not consistently present in the top 30 ranking (aspects which are not considered in the approaches usually adopted by most statistical agencies). The multi-period version of the MPL price index proves suitable for this scope. Hence, it has been applied to the data set which collects the number of visitors and revenues of the 30 leading Italian cultural institutions which, from 2004 to 2017, were ranked at least two times in the top 30.\footnote{All analyses in this investigation have been made with our own codes, written in R.} Figure (ref) shows the MPL price index together with its annual percentage variations for the period 2004--2017: 2004 being the base year and 2017 the year used for updating the index. We can note that in the early years of our new Millennium, when important investments started being made in the Italian cultural sector, the prices of museums (and the like) tickets grew (Figure (ref)). Thereafter, the price dynamic became more moderate and then tapered in 2009 and 2014 when, the so called “W” recession, namely the international financial and debt crisis in European peripheral countries, hit Italy.

figure[figure omitted — 183 chars of source]

For the sake of further evidence from an empirical standpoint, a comparison of the MPL and the TPD index has been conducted. First TPD price indexes have been computed only for those museums whose prices are available at all times. This has led to a drop in the number of museums/monuments/archeological sites from 36 to 17. Figure (ref) (first panel) shows both the MPL and TPD price indexes together with their 3$\sigma$ confidence bounds. The result is that MPL indexes always fall within the confidence bounds of TPD indexes. This result highlights their alignment with the latter, and provides evidence of their greater accuracy, due to their lower standard errors. In the case when not all items (museums) are priced in all periods, TPD estimates have been obtained by using the time version of the weighted CPD rao2004country. Figure (ref) (second panel) shows both MPL and TPD indexes together with their $3\sigma$ confidence bounds. The same comments of the complete price tableau case apply here.

figure[figure omitted — 399 chars of source]

The availability of data on visitors and revenues in 2017 for museums, monuments, archaeological sites, and museum circuits in the North-West, North-East, Centre and South (which includes the two islands Sicily and Sardinia) has allowed the computation of the multilateral version of the MPL index. Looking at the data, we see that almost half (46.3%) are located in the North, while 28.5% in the Centre, and 25.2% in the South and Islands. The Regions with the highest number of cultural institutions are Tuscany (11%), followed by Emilia-Romagna (9.6%), Piedmont (8.6%) and Lombardy (8.2%) istatr2016. However, alongside the more famous attractions, Italy is home to a wide and rich array of notable locations of cultural interest. A considerable percentage of these places (17.5%) are found in municipalities with less than 2,000 inhabitants, but which can have up to four or five cultural sites in their small area. Almost a third (30.7%) are distributed in 1,027 municipalities with a population varying from 2,000 to 10,000, and a bit more than half (51.8%) are situated in 712 municipalities with a population of 10,000 to 50,000. Italy is, therefore, characterized by a strongly polycentric cultural supply distributed throughout its territory, even in areas considered as marginal from a geographic stance. Table (ref) reports, in the first row, the MPL index computed applying Eq. (ref) to the first three areas (North-West, North-East and Centre) considering the Centre as base area. The second row shows the updated values of the MPL index when the South-Islands are added to the data-set. As for the multi-period case, a comparison of the MPL estimates with the CPD ones is provided. The third row of Table (ref) shows CPD estimates in the case of full price tableau, as all commodities are priced in the four geographic areas. Figure (ref) shows both the MPL and CPD indexes together with their 3$\sigma$ confidence bounds. Once again, the estimates of the MPL index turn out to be more accurate than those provided by the CPD approach, as the former have standard errors lower than the latter. As in the comparison with the TPD index, the confidence bounds of CPD indexes always include MPL estimates, thus suggesting the compatibility of both indexes. It is worth noting that in 2017, access to cultural sites in Southern Italy cost the most: almost twice as much as in the North-Eastern area. While the disparity could be ascribed to several factors, such as different costs of managing museums and similar institutions, tourism flows, etc: that type of analysis goes beyond the scope of the current investigation.

table[table omitted — 511 chars of source]
figure[figure omitted — 167 chars of source]

Finally, in order to investigate more thoroughly the performance of the MPL index as compared to the TPD one, a simulation analysis has been performed. One thousand simulations were carried out by using perturbed values (and prices as a by-product) and assuming fixed quantities (i.e. equal to the original ones). Next, the simulated values (and prices) were used to compute MPL and TPD indexes in different settings: with and without missing values and/or quantities (and accordingly prices). The final MPL and TPD indexes were obtained as averages of all indexes computed on simulated values and prices. Two types of simulations were carried out. First simulated values from the $2^{nd}$ to the $T^{th}$ period (base period values, $\boldsymbol{v}_1$, being kept fixed) were obtained from the original ones ($\boldsymbol{V}_1$) by adding random terms drawn from Normal laws with different means and variances. Plots in Figure (ref) show both the MPL and the TPD indexes, for complete and incomplete price tableau cases, together with the associated $3\sigma$ confidence bounds.

figure[figure omitted — 466 chars of source]

Then, simulated values, from the $2^{nd}$ to the $T^{th}$ period (base period values, $\boldsymbol{v}_1$, being kept fixed), were obtained from simulated values of the previous period with the addition of error terms drawn from Normal laws with given means and variances. Plots in Figure (ref) show both the MPL and the TPD indexes, for the case of complete and incomplete price tableaus, together with the associated $3\sigma$ confidence bounds. In both cases, the MPL estimates are in line with the TPD ones, but are more accurate than the latter as their tighter confidence bounds show.

figure[figure omitted — 526 chars of source]

Conclusion

The paper works out a novel price index that can be used either as a multi-period or as a multilateral index. This index, called MPL index, is obtained as a solution to an “ad hoc” minimum-norm criterion, within the framework of the stochastic approach. The computation of the MPL index does not require the knowledge of commodity prices, but only their quantities and values. The reference basket of the MPL index, over periods or across countries, is more informative and complete than the ones commonly used by statistical agencies, and easy to update. The updating process is twofold depending on the multi-period or the multilateral use of the index. An application of the MPL index to the Italian cultural supply data provides proof of its positive performance. A comparison between the MPL and the CPD/TPD index on both real and simulated data provides evidence of the greater accuracy of the MPL estimates.