<p>Unfortunately the paper is behind a paywall, so the answer is based entirely on your question.</p>
<p>What is done under the hood to get the result is a rather straight forward log-linearisation and, hence, approximation. If your outline is correct, the paper uses the &quot;tricks&quot; of taken continuous time derivative of discrete functions and discrete integration
derivation and discrete integration, but ultimately it comes down to the same thing.</p>
<p>I'll provide an alternative derivation of the result, making the approximation explicit as well as identifying the weights <span class="math-container">$w$</span>. So even if it's not the same, you can follow and confirm the result.</p>
<p>We use that we can approximate a variable <span class="math-container">$X_t$</span> as
<span class="math-container">$$
X_t \approx X\biggl(1+\log\frac{X_t}{X}\biggr)
$$</span>
where <span class="math-container">$X$</span> is a constant. The closer <span class="math-container">$X$</span> is to <span class="math-container">$X_t$</span> the better the approximation. Note that <span class="math-container">$V(T)$</span> may be far from <span class="math-container">$V(0)$</span> in which case the approximation is bad, say if the rate of change of <span class="math-container">$V$</span> from <span class="math-container">$0$</span> to <span class="math-container">$T$</span> is more than 10%.</p>
<p>Let's start by approximating <span class="math-container">$V_j(T)$</span> around <span class="math-container">$V_j(0)$</span>. Apart from the definitions, we use that <span class="math-container">$\log \frac{x}{y}=\log x - \log y$</span>. We get
<span class="math-container">\begin{equation*}
  \begin{split}
     V_j (T) &amp; \approx V_j(0) \biggl( 1+\log \frac{V_j(T)}{V_j(0)}\biggr)= V_j(0)\biggl(1+\log V_j(T) - \log V_j(0)\biggr) \\
       &amp; = V_j(0)\biggl(1+\log \prod_{i=1}^n x_{ij}(T) - \log \prod_{i=1}^n x_{ij}(0)\biggr) \\
       &amp; = V_j(0)\biggl(1+ \sum_{i=1}^n \log x_{ij}(T) - \sum_{i=1}^n \log x_{ij}(0)\biggr) \\
       &amp; = V_j(0)\biggl(1+ \sum_{i=1}^n [\log x_{ij}(T) - \log x_{ij}(0) ] \biggr) \\
       &amp; = V_j(0)\biggl(1+ \sum_{i=1}^n \log \frac{x_{ij}(T)}{x_{ij}(0)} \biggr)
   \end{split}
\end{equation*}</span></p>
<p>Subtracting <span class="math-container">$V_j(0)$</span> from both sides, we get
<span class="math-container">\begin{equation*}
  \begin{split}
V_j (T) -  V_j(0) &amp;\approx V_j(0)\biggl(1+ \sum_{i=1}^n \log \frac{x_{ij}(T)}{x_{ij}(0)} \biggr) - V_j(0) \\
    &amp; = V_j(0)\sum_{i=1}^n \log \frac{x_{ij}(T)}{x_{ij}(0)}
   \end{split}
\end{equation*}</span>
At this point you can already recognize similarities in the structure. Now, all that remains to do is to sum up, simplify, and identify the weights.
<span class="math-container">\begin{equation*}
  \begin{split}
V(T)-V(0) &amp; = \sum_{j=1}^m [V_j (T) -  V_j(0)] \\
    &amp; \approx \sum_{j=1}^m V_j(0)\sum_{i=1}^n \log \frac{x_{ij}(T)}{x_{ij}(0)}\\
    &amp; = \sum_{j=1}^m \biggl(V_j(0) \log \frac{x_{1j}(T)}{x_{ij}(0)}+ \cdots + V_j(0)\log \frac{x_{ij}(T)}{x_{nj}(0)} \biggr)
   \end{split}
\end{equation*}</span></p>
<p>Note that the last line correspond exactly to the result in your question, with the weights identified as <span class="math-container">$w_j=V_j(0)$</span>. So both derivations are equivalent.</p>

CC BY-SA 4.0
2022-12-09T14:18:14.960Z
Looks legit. The only further thing I can contribute is the paragraph from the paper following the one written down in the question. This might be of further help: "In the literature, different sets of weights have been proposed by researchers. This has led to several Divisia decomposition methods: Parametric Divisia method 1 (PDM1) (Ang, 1994, Ang and Lee, 1994), AMDI (Boyd et al., 1987), LMDI-I (Ang and Liu, 2001), and LMDI-II (Ang and Choi, 1997). PDM1 and AMDI are parametric, while LMDI-I and LMDI-II are non-parametric and do not rely on any specification for the weights."
VARulle
CC BY-SA 4.0
2022-12-09T14:18:50.283Z
"A summary of the decomposition formulae for AMDI, LMDI-I, and LMDI-II can be found in Ang et al. (2003) and Ang (2004). Eqs. (9), (10) imply that the observed arithmetic and ratio changes of the aggregate indicator may not be exactly the same as the sum and product of the estimated effects, respectively. A residual term may appear in the decomposition results. In IDA, a perfect decomposition method is one that does not give a residual term. Among the above Divisia methods, LMDI-I and LMDI-II are the only ones that give perfect decomposition results at the aggregate level (Ang et al., 2003)."
VARulle
CC BY-SA 4.0
2022-12-09T15:57:41.367Z
@VARulle: appreciate that. For me the paper is behind a paywall.
BrsG
CC BY-SA 4.0
2022-12-07T13:55:15.323Z
2022-12-08T14:44:43.547Z
2022-12-08T14:44:43.547Z
BrsG
<blockquote>
<p>Can anyone point me in the right direction to understand this better?</p>
</blockquote>
<p>I give you a partial answer, that I hope can help you to have a guidance to understand the matter.</p>
<p>Unfortunately, the paper you linked is not readable, so I read only the abstract.  And we can't know a priori the particular method used there to calculate the integral. So, without reading the paper, a complete, specific, answer is difficult to give.</p>
<blockquote>
<p>the paper states that as the data we want to evaluate the
integral over <strong>is not continuous (which it's not)</strong> we need to use a
discrete integration</p>
<p>[...] I've no idea <strong>how we went from the integral to the sum</strong>, using
some discrete integration, nor what these weights are that have been
introduced.</p>
</blockquote>
<p>As I suppose you know, besides analytical methods to solve integrals, there are <em>numerical</em> or <em>discrete</em> methods to calculate integrals.</p>
<p>These methods are important in two particular cases:</p>
<ol>
<li><p>When the integral cannot be calculated via analytical methods, or, as it is said, in closed form: many integrals cannot be calculated or doesn't exist a closed form (a formula) for the integral, or better, for the primitive.</p>
</li>
<li><p>We have data. Data are <strong>by nature discontinuos</strong>, you have specific numerical values indexed, for instance, by time. That is, you <strong>don't know the function</strong> to be integrated, you don't have it, but only some specific, discrete values of it.
Without knowing the function, of course, we cannot make the integral, so we have to resort to numerical or discrete methods.</p>
</li>
</ol>
<p>Also discontinuos functions can be integrable, and in most cases they are, but the fact is that in your case <strong>we have no function, but only points of the function.</strong></p>
<p>That's the reason why you pass from the integral to a sum, because a sum is on discrete values, integration requires a function (and that function must be integrable).</p>
<p>The subject of numerical or discrete integration is part of <em>numerical analysis.</em></p>
<p>There are many numerical methods to solve integrals.</p>
<p>We should know the specific method used in your paper, so that we can understand also what these 'weights' <span class="math-container">$w$</span> are in that particular case:</p>
<blockquote>
<p>Numerical integration methods can generally be described as combining
evaluations of the integrand to get an approximation to the integral.
The integrand is evaluated at a finite set of points called
integration points and a weighted sum of these values is used to
approximate the integral. The integration points and weights depend on
the specific method used and the accuracy required from the
approximation.</p>
</blockquote>
<p>This quote is from an article of Wikipedia on numerical integration:</p>
<p><a href="https://en.wikipedia.org/wiki/Numerical_integration" rel="nofollow noreferrer">https://en.wikipedia.org/wiki/Numerical_integration</a></p>
<p>From a mathematical point of view, if one could read the original paper, we could  see if , in numerical integration, some interpolatory method has been used, as Lagrange polynomials to find the weights <span class="math-container">$w_j$</span>. In this case the interpolation points are given by the available observations.</p>
<p>Without reading the paper, it is just a guess.
But, as I understand from the comments of Varulle, who quotes parts of the original article, this is not the case, other methods are used.</p>
<p>To continue the guessing, it is enough clear that a numerical quadrature has been used, that is a formula as</p>
<p><span class="math-container">$\int_0^T f(t) = \sum_ {j=1,..,n}w_jf(x_j)$</span> ,</p>
<p>what is not clear is how the derivatives <span class="math-container">$\frac{dln(x_{1j}(t))}{dt}$</span> disappear.</p>

CC BY-SA 4.0
2022-12-06T13:39:54.420Z
2022-12-09T19:58:55.857Z
2022-12-09T19:58:55.857Z
BakerStreet
<p>Sorry for this down-to-earth answer, but I think that the underlying assumption simply states that <span class="math-container">$V_j(t)$</span> is constant (up to an uncorrelated random term), or in other words:
<span class="math-container">$$ V_j(t)=w_j,$$</span>
which implies that
<span class="math-container">$$ \int_0^T V_j(t) \frac{dln(x_{1j}(t))}{dt} dt = w_j\int_0^T \frac{dln(x_{1j}(t))}{dt} dt= w_jln\left(\frac{x_{1j}(T)}{x_{1j}(0)}\right),$$</span>
and yields the expression for <span class="math-container">$V(T)-V(0)$</span>. I do not see in the paper any effort for using numerical integration methods.</p>
<p><strong>EDIT</strong>: If we try to go beyond constant values for <span class="math-container">$V_j$</span>, the result will be different from the one given by the authors, in general. For instance, in the linear case, <span class="math-container">$V(t)=a_j+b_jt$</span>, with <span class="math-container">$V'(t)=b_j \neq 0,$</span> we find that</p>
<p><span class="math-container">\begin{align*}
\int_0^T \left( a_j+b_js \right) \frac{d\ln(x_{1j})}{dt}(s) ds &amp;= a_j\int_0^T \frac{d\ln(x_{1j})}{dt}(s) ds + b_j\int_0^T s \frac{d\ln(x_{1j})}{dt}(s) ds \\
&amp;= a_j\ln\left(\frac{x_{1j}(T)}{x_{1j}(0)}\right) 
+ b_j G_j(T,0), 
\end{align*}</span>
which is generally different from the equation given by the authors. If in top of that, we either assume that <span class="math-container">$\ln(x_{1j}(t))' = c_{1j},$</span> or that <span class="math-container">$x_{1j}'(t) = c_{1j},$</span>  we still end up with an expression with <span class="math-container">$G_j \neq 0$</span> and NOT linear in <span class="math-container">$ \ln \left(x_{1j}(T)/x_{1j}(0)\right) $</span>.</p>
<p>The authors explicitly say that their equation is an approximation, and its quality will depend on the empirical properties of the data. A possibility would be to determine the weights by a regression, provided sufficient observations are available.</p>

CC BY-SA 4.0
2022-12-07T10:20:12.717Z
To tell the truth, I don’t think that the $V_j$ are constant. In the question it is written that 
$\frac{dV(t)}{d(t)} = \sum_j^m \frac{dV_j(t)}{d(t)}$
Therefore, if $V_j(t)$ would be constant, the above derivatives should be $0$, so end of the story. 
I see that that $w_{ij}$ come from the discrete integration, are not an initial hypothesis abot $V_j$ being constant.
BakerStreet
CC BY-SA 4.0
2022-12-07T12:06:49.113Z
Of course, $V_j(t)$ always gives the same value for a specific $t$. But the whole purpose is to see what happens if $t$ changes. The authors abuse notation with respect to derivatives and integrals, which can be confusing when the function is really discrete.
BrsG
CC BY-SA 4.0
2022-12-07T13:00:22.707Z
It sems to me that the author thinks of $()$ as a function of time defined on an interval, say $[0,∞)$, so the integral makes sense. But we have data, discrete, not the function itself, that is only punctual values of the function $()$ . Therefore, the author must resort to discrete integration methods.
BakerStreet
CC BY-SA 4.0
2022-12-07T13:48:59.417Z
In other words, if one wants to speak of a discrete function, we can say that we have a function $V(t)$ defined on an interval $[0, \infty)$, but as we  have discrete data, we know only a function,  let’s call it $V’(t)$, defined on a discrete set, with the same values of $V(t)$ on this set. The $V’(t)$ can be thought of as a restriction to these discrete values of $V(t)$. $V(t)$ is integrable, $V’(t)$ of course not.
BakerStreet
CC BY-SA 4.0
2022-12-07T14:16:13.967Z
I agree with your comments and have added a paragraph to address them.
Bertrand
CC BY-SA 4.0
2022-12-07T08:02:18.000Z
2022-12-07T14:15:04.267Z
2022-12-07T14:15:04.267Z
Bertrand
<p>I’m currently reading a paper on index decomposition. The paper is here for reference : <a href="https://www.sciencedirect.com/science/article/pii/S0140988315001772" rel="nofollow noreferrer">https://www.sciencedirect.com/science/article/pii/S0140988315001772</a></p>
<p>The paper is setting out how it has gone about deriving the index decomposition, and I’m getting a bit lost at the final step. The aim is to decompose the changes in an aggregate variable between time <span class="math-container">$T$</span> and time <span class="math-container">$0$</span>.</p>
<p>We start from the position that we want to break an aggregate category <span class="math-container">$V$</span> down into a series <span class="math-container">$j$</span> of subcategories, <span class="math-container">$V_j$</span>, each of which is impacted by a range of <span class="math-container">$i$</span> factors over time. So we have :</p>
<p><span class="math-container">$$ V(t) = \sum_j^m V_j (t) =  \sum_j^m \big ( \prod_1^i x_{ij}(t)  \big )$$</span></p>
<p>So what we are interested in is:
<span class="math-container">$$V(T) - V(0) = \sum_j^m \big ( \prod_1^i x_{ij}(T) \big ) - \sum_j^m \big ( \prod_1^i x_{ij}(0) \big ) $$</span></p>
<p>We can also consider <span class="math-container">$V(T) - V(0)$</span> as:
<span class="math-container">$$V(T) - V(0) = \int_0^T \frac{dV(t)}{d(t)}$$</span></p>
<p>Which then leads us to consider <span class="math-container">$\frac{dV(t)}{d(t)}$</span>, from the definition above we have:</p>
<p><span class="math-container">$$\frac{dV(t)}{d(t)} = \sum_j^m \frac{dV_j(t)}{d(t)}$$</span></p>
<p><span class="math-container">$$=\sum_j^m \left( \frac{dx_{1j}(t)}{dt}x_{2j}(t)x_{3j}(t)...x_{nj}(t)+\frac{dx_{2j}(t)}{dt}x_{1j}x_{3j}(t)...x_{nj}(t)+... +\frac{dx_{nj}(t)}{dt}x_{1j}(t)x_{2j}(t)...x_{(n-1)j} \right)$$</span></p>
<p>We can then re-write this as:</p>
<p><span class="math-container">$$\frac{dV(t)}{dt} =\sum_j^m \left(  \frac{V_j(t)}{x_{1j}(t)}\frac{dx_{1j}(t)}{dt}+\frac{V_j(t)}{x_{2j}(t)}\frac{dx_{2j}(t)}{dt} + ... + \frac{V_j(t)}{x_{nj}(t)}\frac{dx_{nj}(t)}{dt} \right) $$</span></p>
<p>Which, as  <span class="math-container">$\frac{dln(x_{ij}(t))}{dt}) = \frac{1}{x_{ij}(t)}\frac{dx_{ij}(t)}{dt}$</span> we can re-write as:</p>
<p><span class="math-container">$$ \frac{dV(t)}{dt} = \sum_j^m V_j(t) \left( \frac{dln(x_{1j}(t))}{dt}+\frac{dlnx_{2j}(t))}{dt} +\frac{dln(x_{3j}(t))}{dt}+ ... + \frac{dln(x_{nj}(t))}{dt} \right) $$</span></p>
<p>Which gives us finally that:
<span class="math-container">$$V(T)-V(0) = \int_0^T \frac{dV(t)}{dt} = \int_0^T \left(\sum_j^m V_j(t) \left( \frac{dln(x_{1j}(t))}{dt}+\frac{dlnx_{2j}(t))}{dt} +\frac{dln(x_{3j}(t))}{dt}+ ... + \frac{dln(x_{nj}(t))}{dt} \right) \right)$$</span></p>
<p>It is at this point that the paper state that as the data we want to evaluate the integral over is not continuous (which it's not) we need to use a discrete integration, which gives the solution:
<span class="math-container">$$V(T)-V(0) = \sum_j^m \left( w_j \times ln \left( \frac{x_{1j}(T)}{x_{1j}(0)} \right) + w_j \times ln \left( \frac{x_{2j}(T)}{x_{2j}(0)} \right) + ... + w_j \times ln \left( \frac{x_{nj}(T)}{x_{nj}(0)} \right)  \right)$$</span></p>
<p>For some weights <span class="math-container">$w_j$</span>.</p>
<p>At this point I'm lost. I've no idea how we went from the integral to the sum, using some discrete integration, nor what these weights are that have been introduced.</p>
<p>Can anyone point me in the right direction to understand this better?</p>
<p>Thanks for any help,</p>
<p>Hmmm16</p>

<mathematical-economics><applied-econometrics>
Derivation of index decomposition analysis
hmmmm
# Derivation of index decomposition analysis

(Asked by: username_0 on 2022-12-03)

I’m currently reading a paper on index decomposition. The paper is here for reference : <https://www.sciencedirect.com/science/article/pii/S0140988315001772>

The paper is setting out how it has gone about deriving the index decomposition, and I’m getting a bit lost at the final step. The aim is to decompose the changes in an aggregate variable between time $T$ and time $0$.

We start from the position that we want to break an aggregate category $V$ down into a series $j$ of subcategories, $V\_j$, each of which is impacted by a range of $i$ factors over time. So we have :

$$ V(t) = \sum\_j^m V\_j (t) = \sum\_j^m \big ( \prod\_1^i x\_{ij}(t) \big )$$

So what we are interested in is:
$$V(T) - V(0) = \sum\_j^m \big ( \prod\_1^i x\_{ij}(T) \big ) - \sum\_j^m \big ( \prod\_1^i x\_{ij}(0) \big ) $$

We can also consider $V(T) - V(0)$ as:
$$V(T) - V(0) = \int\_0^T \frac{dV(t)}{d(t)}$$

Which then leads us to consider $\frac{dV(t)}{d(t)}$, from the definition above we have:

$$\frac{dV(t)}{d(t)} = \sum\_j^m \frac{dV\_j(t)}{d(t)}$$

$$=\sum\_j^m \left( \frac{dx\_{1j}(t)}{dt}x\_{2j}(t)x\_{3j}(t)...x\_{nj}(t)+\frac{dx\_{2j}(t)}{dt}x\_{1j}x\_{3j}(t)...x\_{nj}(t)+... +\frac{dx\_{nj}(t)}{dt}x\_{1j}(t)x\_{2j}(t)...x\_{(n-1)j} \right)$$

We can then re-write this as:

$$\frac{dV(t)}{dt} =\sum\_j^m \left( \frac{V\_j(t)}{x\_{1j}(t)}\frac{dx\_{1j}(t)}{dt}+\frac{V\_j(t)}{x\_{2j}(t)}\frac{dx\_{2j}(t)}{dt} + ... + \frac{V\_j(t)}{x\_{nj}(t)}\frac{dx\_{nj}(t)}{dt} \right) $$

Which, as $\frac{dln(x\_{ij}(t))}{dt}) = \frac{1}{x\_{ij}(t)}\frac{dx\_{ij}(t)}{dt}$ we can re-write as:

$$ \frac{dV(t)}{dt} = \sum\_j^m V\_j(t) \left( \frac{dln(x\_{1j}(t))}{dt}+\frac{dlnx\_{2j}(t))}{dt} +\frac{dln(x\_{3j}(t))}{dt}+ ... + \frac{dln(x\_{nj}(t))}{dt} \right) $$

Which gives us finally that:
$$V(T)-V(0) = \int\_0^T \frac{dV(t)}{dt} = \int\_0^T \left(\sum\_j^m V\_j(t) \left( \frac{dln(x\_{1j}(t))}{dt}+\frac{dlnx\_{2j}(t))}{dt} +\frac{dln(x\_{3j}(t))}{dt}+ ... + \frac{dln(x\_{nj}(t))}{dt} \right) \right)$$

It is at this point that the paper state that as the data we want to evaluate the integral over is not continuous (which it's not) we need to use a discrete integration, which gives the solution:
$$V(T)-V(0) = \sum\_j^m \left( w\_j \times ln \left( \frac{x\_{1j}(T)}{x\_{1j}(0)} \right) + w\_j \times ln \left( \frac{x\_{2j}(T)}{x\_{2j}(0)} \right) + ... + w\_j \times ln \left( \frac{x\_{nj}(T)}{x\_{nj}(0)} \right) \right)$$

For some weights $w\_j$.

At this point I'm lost. I've no idea how we went from the integral to the sum, using some discrete integration, nor what these weights are that have been introduced.

Can anyone point me in the right direction to understand this better?

Thanks for any help,

Hmmm16

**Tags:** mathematical-economics, applied-econometrics

---

## 3 answers

### Answer by username_1 on 2022-12-07. Score: 6

Unfortunately the paper is behind a paywall, so the answer is based entirely on your question.

What is done under the hood to get the result is a rather straight forward log-linearisation and, hence, approximation. If your outline is correct, the paper uses the "tricks" of taken continuous time derivative of discrete functions and discrete integration
derivation and discrete integration, but ultimately it comes down to the same thing.

I'll provide an alternative derivation of the result, making the approximation explicit as well as identifying the weights $w$. So even if it's not the same, you can follow and confirm the result.

We use that we can approximate a variable $X\_t$ as
$$
X\_t \approx X\biggl(1+\log\frac{X\_t}{X}\biggr)
$$
where $X$ is a constant. The closer $X$ is to $X\_t$ the better the approximation. Note that $V(T)$ may be far from $V(0)$ in which case the approximation is bad, say if the rate of change of $V$ from $0$ to $T$ is more than 10%.

Let's start by approximating $V\_j(T)$ around $V\_j(0)$. Apart from the definitions, we use that $\log \frac{x}{y}=\log x - \log y$. We get
\begin{equation\*}
\begin{split}
V\_j (T) & \approx V\_j(0) \biggl( 1+\log \frac{V\_j(T)}{V\_j(0)}\biggr)= V\_j(0)\biggl(1+\log V\_j(T) - \log V\_j(0)\biggr) \\
& = V\_j(0)\biggl(1+\log \prod\_{i=1}^n x\_{ij}(T) - \log \prod\_{i=1}^n x\_{ij}(0)\biggr) \\
& = V\_j(0)\biggl(1+ \sum\_{i=1}^n \log x\_{ij}(T) - \sum\_{i=1}^n \log x\_{ij}(0)\biggr) \\
& = V\_j(0)\biggl(1+ \sum\_{i=1}^n [\log x\_{ij}(T) - \log x\_{ij}(0) ] \biggr) \\
& = V\_j(0)\biggl(1+ \sum\_{i=1}^n \log \frac{x\_{ij}(T)}{x\_{ij}(0)} \biggr)
\end{split}
\end{equation\*}

Subtracting $V\_j(0)$ from both sides, we get
\begin{equation\*}
\begin{split}
V\_j (T) - V\_j(0) &\approx V\_j(0)\biggl(1+ \sum\_{i=1}^n \log \frac{x\_{ij}(T)}{x\_{ij}(0)} \biggr) - V\_j(0) \\
& = V\_j(0)\sum\_{i=1}^n \log \frac{x\_{ij}(T)}{x\_{ij}(0)}
\end{split}
\end{equation\*}
At this point you can already recognize similarities in the structure. Now, all that remains to do is to sum up, simplify, and identify the weights.
\begin{equation\*}
\begin{split}
V(T)-V(0) & = \sum\_{j=1}^m [V\_j (T) - V\_j(0)] \\
& \approx \sum\_{j=1}^m V\_j(0)\sum\_{i=1}^n \log \frac{x\_{ij}(T)}{x\_{ij}(0)}\\
& = \sum\_{j=1}^m \biggl(V\_j(0) \log \frac{x\_{1j}(T)}{x\_{ij}(0)}+ \cdots + V\_j(0)\log \frac{x\_{ij}(T)}{x\_{nj}(0)} \biggr)
\end{split}
\end{equation\*}

Note that the last line correspond exactly to the result in your question, with the weights identified as $w\_j=V\_j(0)$. So both derivations are equivalent.

> **Comments:**
> **username_2**: Looks legit. The only further thing I can contribute is the paragraph from the paper following the one written down in the question. This might be of further help: "In the literature, different sets of weights have been proposed by researchers. This has led to several Divisia decomposition methods: Parametric Divisia method 1 (PDM1) (Ang, 1994, Ang and Lee, 1994), AMDI (Boyd et al., 1987), LMDI-I (Ang and Liu, 2001), and LMDI-II (Ang and Choi, 1997). PDM1 and AMDI are parametric, while LMDI-I and LMDI-II are non-parametric and do not rely on any specification for the weights."  
> **username_2**: "A summary of the decomposition formulae for AMDI, LMDI-I, and LMDI-II can be found in Ang et al. (2003) and Ang (2004). Eqs. (9), (10) imply that the observed arithmetic and ratio changes of the aggregate indicator may not be exactly the same as the sum and product of the estimated effects, respectively. A residual term may appear in the decomposition results. In IDA, a perfect decomposition method is one that does not give a residual term. Among the above Divisia methods, LMDI-I and LMDI-II are the only ones that give perfect decomposition results at the aggregate level (Ang et al., 2003)."  
> **username_1**: @VARulle: appreciate that. For me the paper is behind a paywall.  

---

### Answer by username_3 on 2022-12-06. Score: 4

> Can anyone point me in the right direction to understand this better?

I give you a partial answer, that I hope can help you to have a guidance to understand the matter.

Unfortunately, the paper you linked is not readable, so I read only the abstract. And we can't know a priori the particular method used there to calculate the integral. So, without reading the paper, a complete, specific, answer is difficult to give.

> the paper states that as the data we want to evaluate the
> integral over **is not continuous (which it's not)** we need to use a
> discrete integration
>
> [...] I've no idea **how we went from the integral to the sum**, using
> some discrete integration, nor what these weights are that have been
> introduced.

As I suppose you know, besides analytical methods to solve integrals, there are *numerical* or *discrete* methods to calculate integrals.

These methods are important in two particular cases:

1. When the integral cannot be calculated via analytical methods, or, as it is said, in closed form: many integrals cannot be calculated or doesn't exist a closed form (a formula) for the integral, or better, for the primitive.
2. We have data. Data are **by nature discontinuos**, you have specific numerical values indexed, for instance, by time. That is, you **don't know the function** to be integrated, you don't have it, but only some specific, discrete values of it.
   Without knowing the function, of course, we cannot make the integral, so we have to resort to numerical or discrete methods.

Also discontinuos functions can be integrable, and in most cases they are, but the fact is that in your case **we have no function, but only points of the function.**

That's the reason why you pass from the integral to a sum, because a sum is on discrete values, integration requires a function (and that function must be integrable).

The subject of numerical or discrete integration is part of *numerical analysis.*

There are many numerical methods to solve integrals.

We should know the specific method used in your paper, so that we can understand also what these 'weights' $w$ are in that particular case:

> Numerical integration methods can generally be described as combining
> evaluations of the integrand to get an approximation to the integral.
> The integrand is evaluated at a finite set of points called
> integration points and a weighted sum of these values is used to
> approximate the integral. The integration points and weights depend on
> the specific method used and the accuracy required from the
> approximation.

This quote is from an article of Wikipedia on numerical integration:

<https://en.wikipedia.org/wiki/Numerical_integration>

From a mathematical point of view, if one could read the original paper, we could see if , in numerical integration, some interpolatory method has been used, as Lagrange polynomials to find the weights $w\_j$. In this case the interpolation points are given by the available observations.

Without reading the paper, it is just a guess.
But, as I understand from the comments of Varulle, who quotes parts of the original article, this is not the case, other methods are used.

To continue the guessing, it is enough clear that a numerical quadrature has been used, that is a formula as

$\int\_0^T f(t) = \sum\_ {j=1,..,n}w\_jf(x\_j)$ ,

what is not clear is how the derivatives $\frac{dln(x\_{1j}(t))}{dt}$ disappear.


---

### Answer by username_4 on 2022-12-07. Score: 4

Sorry for this down-to-earth answer, but I think that the underlying assumption simply states that $V\_j(t)$ is constant (up to an uncorrelated random term), or in other words:
$$ V\_j(t)=w\_j,$$
which implies that
$$ \int\_0^T V\_j(t) \frac{dln(x\_{1j}(t))}{dt} dt = w\_j\int\_0^T \frac{dln(x\_{1j}(t))}{dt} dt= w\_jln\left(\frac{x\_{1j}(T)}{x\_{1j}(0)}\right),$$
and yields the expression for $V(T)-V(0)$. I do not see in the paper any effort for using numerical integration methods.

**EDIT**: If we try to go beyond constant values for $V\_j$, the result will be different from the one given by the authors, in general. For instance, in the linear case, $V(t)=a\_j+b\_jt$, with $V'(t)=b\_j \neq 0,$ we find that

\begin{align\*}
\int\_0^T \left( a\_j+b\_js \right) \frac{d\ln(x\_{1j})}{dt}(s) ds &= a\_j\int\_0^T \frac{d\ln(x\_{1j})}{dt}(s) ds + b\_j\int\_0^T s \frac{d\ln(x\_{1j})}{dt}(s) ds \\
&= a\_j\ln\left(\frac{x\_{1j}(T)}{x\_{1j}(0)}\right)
+ b\_j G\_j(T,0),
\end{align\*}
which is generally different from the equation given by the authors. If in top of that, we either assume that $\ln(x\_{1j}(t))' = c\_{1j},$ or that $x\_{1j}'(t) = c\_{1j},$ we still end up with an expression with $G\_j \neq 0$ and NOT linear in $ \ln \left(x\_{1j}(T)/x\_{1j}(0)\right) $.

The authors explicitly say that their equation is an approximation, and its quality will depend on the empirical properties of the data. A possibility would be to determine the weights by a regression, provided sufficient observations are available.

> **Comments:**
> **username_3**: To tell the truth, I don’t think that the $V_j$ are constant. In the question it is written that  $\frac{dV(t)}{d(t)} = \sum_j^m \frac{dV_j(t)}{d(t)}$ Therefore, if $V_j(t)$ would be constant, the above derivatives should be $0$, so end of the story.  I see that that $w_{ij}$ come from the discrete integration, are not an initial hypothesis abot $V_j$ being constant.  
> **username_1**: Of course, $V_j(t)$ always gives the same value for a specific $t$. But the whole purpose is to see what happens if $t$ changes. The authors abuse notation with respect to derivatives and integrals, which can be confusing when the function is really discrete.  
> **username_3**: It sems to me that the author thinks of $()$ as a function of time defined on an interval, say $[0,∞)$, so the integral makes sense. But we have data, discrete, not the function itself, that is only punctual values of the function $()$ . Therefore, the author must resort to discrete integration methods.  
> **username_3**: In other words, if one wants to speak of a discrete function, we can say that we have a function $V(t)$ defined on an interval $[0, \infty)$, but as we  have discrete data, we know only a function,  let’s call it $V’(t)$, defined on a discrete set, with the same values of $V(t)$ on this set. The $V’(t)$ can be thought of as a restriction to these discrete values of $V(t)$. $V(t)$ is integrable, $V’(t)$ of course not.  
> **username_4**: I agree with your comments and have added a paragraph to address them.  

---

