Chapitre 5 Analyses longitudinales

library(plm)

library(survival)
library(survminer)

5.1 Pseudo-panels

5.2 Modèles linéaires

Modèle de base: \[y_{i,t} = \alpha_i + \beta_i k_{i,t} + \gamma_i n_{i,t} + \epsilon_{i,t}\]

lineaire <- lm(y ~ x1 + X2, data=d)

Modèle pooled13: \[y_{i,t} = \alpha + \beta k_{i,t} + \gamma n_{i,t} + \epsilon_{i,t}\]

pooled <- plm(y ~ x1 + X2,
              data = d,
              index = c("var_individu","var_temporelle"),
              model = "pooling")

Modèle à effet fixe: \[y_{i,t} = \alpha_i + \beta k_{i,t} + \gamma n_{i,t} + \epsilon_{i,t}\]

effet_fixe <- plm(y ~ x1 + x2,
                  data,
                  index,
                  model = "within")

# twoways: estimer l'effet fixe individuelle et l'effet fixe temporel
model_sete <- plm(y ~ x1 + x2,
                  data,
                  index,
                  model = "within",
                  effect = "twoways")

Tests

## Test de Fischer ##
# Test sur la pertinence de l'hypothèse de l'existence d'un effet fixe individuel (s'applique aussi au modèle pooled). Seule la p-value nous intéresse.
coeftest(effet_fixe, vcov. = vcovHC, type = "HC1")

## pFtest: F-test pour comparer le pouvoir explicatif du modèle A par rapport au modèle B ##
pFtest(effet_fixe, pooled)

## Test de Hausman ##
phtest(effet_fixe, pooled)

5.3 Modèles dynamiques

plm(y ~ lag(y) + lag(X1) + x2,
    data,
    index,
    model = "within", # ou autre
    effect = "individual") # ou autre

5.4 Modèles non-linéaires

5.5 Modèles de durée

5.6 Optimal matching


  1. Moindres Carrés Ordinaires askip???↩︎