serie chrono et extraction

Bonjour

Je dispose d'une série de données sur plusieurs années où je veux mesurer INDI:

X              INDI
1994         4.36
1998         3.26
1999          5.63
....            ...

J'utilise le test de Bai et Perron afin qu'il m estime les breakpoints (estimations des points où il y a une tendance inverse).
J'obtiens ce résultat :

Corresponding to breakdates:
  2.5 % breakpoints 97.5 %
1  1994       1995   1996
2  2003        2004   2005

Je souhaite créer une colonne dummy_variable dans mon fichier qui prenne 1 si l’année est un breakpoint ou 0 si elle ne l'est pas (donc 1 uniquement pour 1994 et 2003).

Je ne vois pas comment extraire le résultat de cet output pour la création de la nouvelle variable bien que je connaisse a peu près la syntaxe attendue :
data$dummy <- ifelse(data$year est breakpoint, 1, 0) quelque chose comme ça.


Merci pour votre aide mes chers!

Code:

---

1 2 3

data <- data.frame(Ma_Serie_Chronologique) data$dummy_variable[data$X == 1994 || data$X == 2003] <- 1 data$dummy_variable[data$X != 1994 && data$X != 2003] <- 0

---

Cordialement.

Citation:

---

Envoyé par dev_ggy

Code:

---

1 2 3

data <- data.frame(Ma_Serie_Chronologique) data$dummy_variable[data$X == 1994 || data$X == 2003] <- 1 data$dummy_variable[data$X != 1994 && data$X != 2003] <- 0

---

Cordialement.

---

Bonjour dev et merci pour ta reponse.

Le truc c est que si j ai envie de faire ca pour chaque tableau de x tableaux. je ne connais pas a l avance la valeur des annees resultante. Le but est que je fasse une boucle qui me genere cette dummy en prenant comme variable le resultat des breakpoints générés pour chaque sous tableaux (ici pour chaque pays), contenu dans le meme jeu de donnes (analyse par matrice). Ces annees bp changent selon le pays. Comment pourrai je associer les bp de l output a la variable de test? merci

Il faut que tu fasses un str() de l'objet qui contient tes breakpoints, ce qui te montrera où est stockée l'information que tu veux (ie les deux années des breakpoints).

Ensuite tu n'auras plus qu'à remplacer tes dates par le bon argument de l'objet:

Code:

---

ta$dummy_variable[data$X == objet$arg1 || data$X == objet$arg2] <- 1

---

8-) merci sengou

Voilà que je fais ma boucle et je n'obtiens que des 0 dans ma variable alors qu'il existe des bp donc 1.

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

data<-ts(tabo, start=1988)   row.names="year"<-data ## mettre en label l annee   Fonction.bp<-function(b) {   bp.inf <- breakpoints(year ~ tradevaluein1000usd , data = data [b,], h = 8)   t<-breakdates(confint(bp.inf,breaks=2)) bb<-t[1,1] vv<-t[2,1]     tabo$res <- ifelse(tabo$year == bb | tabo$year == vv , 1, 0)   return( tabo$res)   }   data <- data.frame(data) numero<-1:nrow(tabo) tapply(numero, tabo$code_o, Fonction.bp)

---

J'obtiens :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

$ABW (1) [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [61] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0   $AFG (2) [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [61] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0   $AGO (3) [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 [61] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

---

Merci de votre aide je suis désespéree.

Rappel input :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

code_o year varqt *Dummy a construire 1 1988 375.059 1 1989 3458.656 1 1990 2924.484 etc 2 1988 121861.2 2 1989 125059.8 2 1996 124654.5 etc 3 1989 3120863 3 1990 4082248 3 1991 3290223 etc

---

svp de l'aide sur l optimisation de mon code, meme si je sais que je suis une quiche dans ce domaine sincere merciiiii