Chapitre 5 Analyses longitudinales
library(plm)
library(survival)
library(survminer)5.2 Modèles linéaires
Modèle de base: yi,t=αi+βiki,t+γini,t+ϵi,t
lineaire <- lm(y ~ x1 + X2, data=d)Modèle pooled13: yi,t=α+βki,t+γni,t+ϵi,t
pooled <- plm(y ~ x1 + X2,
data = d,
index = c("var_individu","var_temporelle"),
model = "pooling")Modèle à effet fixe: yi,t=αi+βki,t+γni,t+ϵ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)Moindres Carrés Ordinaires askip???↩︎