Exemple avec R: Indice stratifié

Le calcul d’un indice stratifié est simple en pratique. Une bonne façon de le faire est d’utiliser une simple régression linéaire.

Considérez les données suivantes avec des transactions sur deux périodes pour deux groupes de marchandises.

# Faire quelques données
df <- data.frame(period = c(0, 0, 0, 1, 1, 1, 1, 1), 
                 group = letters[1:2],
                 price = 1:8)
df

##   period group price
## 1      0     a     1
## 2      0     b     2
## 3      0     a     3
## 4      1     b     4
## 5      1     a     5
## 6      1     b     6
## 7      1     a     7
## 8      1     b     8

Les indices pour chaque groupe peuvent être calculés avec une seule régression linéaire, puis agrégés avec une moyenne.

# Apportez dans la bibliothèque gpindex
library(gpindex)

# Calculer les indices au niveau des strates avec une régression linéaire
mdl <- lm(log(price) ~ group + group:period - 1, df)

# Transformer les coefficients de régression en indices
index <- exp(coef(mdl)[-seq_len(length(unique(mdl$model$group)))])

# Agrégat, en supposant que les deux strates ont un poids égal
geometric_mean(index) * 100

## [1] 313.886