Générer une variable fictive

Une autre option qui peut mieux fonctionner si vous avez beaucoup de variables est factor et model.matrix.

> year.f = factor(year)
> dummies = model.matrix(~year.f)

Cela comprendra une colonne d'interception (tous les uns) et une colonne pour chacune des années de votre ensemble de données, sauf une, qui sera la valeur par défaut ou d'interception.

Vous pouvez changer la façon dont la valeur par défaut est choisie en modifiant contrasts.arg dans model.matrix.

De plus, si vous voulez omettre l'intercept, vous pouvez simplement laisser tomber la première colonne ou ajouter +0 à la fin de la formule.

J'espère que cela vous sera utile.

La manière la plus simple de produire ces variables fictives est la suivante :

> print(year)
[1] 1956 1957 1957 1958 1958 1959
> dummy <- as.numeric(year == 1957)
> print(dummy)
[1] 0 1 1 0 0 0
> dummy2 <- as.numeric(year >= 1957)
> print(dummy2)
[1] 0 1 1 1 1 1

Plus généralement, vous pouvez utiliser ifelse pour choisir entre deux valeurs en fonction d'une condition. Ainsi, si au lieu d'une variable fictive 0-1, vous vouliez, pour une raison quelconque, utiliser 4 et 7, vous pourriez utiliser ifelse(année == 1957, 4, 7).

Ce que je fais normalement pour travailler avec ce genre de variables fictives est :

(1) comment puis-je générer une variable muette pour l'observation n°10, c'est-à-dire pour l'année 1957 (valeur = 1 en 1957 et zéro sinon).

data$factor_year_1 <- factor ( with ( data, ifelse ( ( year == 1957 ), 1 , 0 ) ) )

(2) comment générer une variable muette qui est nulle avant 1957 et prend la valeur 1 à partir de 1957 et jusqu'en 2009 ?

data$factor_year_2 <- factor ( with ( data, ifelse ( ( year < 1957 ), 0 , 1 ) ) )

Je peux ensuite introduire ce facteur comme variable muette dans mes modèles. Par exemple, pour voir s'il existe une tendance à long terme dans une variable y :

summary ( lm ( y ~ t,  data = data ) )

J'espère que cela vous aidera !