Discussion :

[sarah_ss]

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!

[dev_ggy]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[sarah_ss]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Sengar]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[sarah_ss]

merci sengou

[sarah_ss]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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