Discussion :

[ANOVA]

Bonjour,

je cherche à faire une image pdf à insérer dans un raport LaTeX. Je souhaite créer cette image avec R. Voici ce que j'ai fait :

Code :

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

1
2
3
4
5
plot(c(0, 100), c(0, 100),type="n", xlab="", ylab="", axes=False)
text(10,97,expression(H[0]~':'~delta~'=0'))
y=c(100,100)
x=c(2,10)
lines(x, y)

Je souhaiterais savoir s'il est possible d'entourer le texte, et si oui comment ? Il y a peut être un paramètre de la fonction text qui le permet, mais je ne l'ai pas trouvé. Quant à utiliser les fonctions lines ou ablines, cela pose problème car dès qu'on agrandit la taille de la fenêtre graphique, les coordonnées qu'on a rentrées ne sont plus valides comme vous pouvez le constater en lançant mon code (les dernières instructions) sur votre ordi.

Je vous remercie.

PS : par ailleurs, je n'ai rien trouvé de mieux que d'utiliser la fonction plot en mettant axes = False pour créer un graphique vierge sur lequel je peut dessiner en entrant les coordonnées. N'y a t-il pas un moyen plus adapté pour cela ?

[A. D.]

Bonjour,

La première idée qui me vient est la suivante :

Code :

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

1
2
3
4
5
 
par(mar=c(16.5,15.5,16.5,15.5))
plot(0,0,type="n",xlab="",ylab="",axes=FALSE)
text(0,0,expression(H[0]~':'~delta~'= 0'))
box(lty=1)

Ce n'est sûrement pas la meilleure façon de procéder, mais je pense que cela donne ce que vous recherchez.
Cela dit, je n'ai pas bien compris l'intérêt : vous souhaitez créer une image pdf du texte " Ho : delta = 0 " encadré pour l'insérer dans un document LaTeX? Pourquoi ne pas écrire directement la formule en LaTeX et utiliser la commande /boxed() (ou quelque chose du genre) pour l'encadrer?

Cordialement,

A.D.

[ANOVA]

Je me connecte rapidement. Je vous remercie pour votre réponse. Je vais regarder ça ce soir.

En fait, j'ai simplifié le problème pour ne pas vous embrouiller avec des détails inutils. Mais ce que je souhaite faire est un peu plus élaboré que ça, sinon, effectivement LaTeX le ferait bien mieux.

En fait, j'ai fait un test de Dickey Fuller, et pour expliquer la séquence des tests, je me suis dit que ce serait pas mal de faire un graphique. Pour l'instant mon graphique donne ceci :

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
plot(c(-10, 90), c(0, 100),type="n", xlab="", ylab="", axes=F)
 
title("Séquence de test ADF")
 
text(10,97,expression(H[0]~':'~delta~'=0'))
text(40,97,expression(H[0]~':'~Phi~'=0'), col = "grey50")
text(70,97,expression(X[t]~'~I(0)'), col = "grey50")
text(40,82,expression(X[t]~'~I(1)'), col = "grey50")
 
text(10,67,expression(H[0]~':'~gamma~'=0'))
text(40,67,expression(H[0]~':'~Phi~'=0'), col = "grey50")
text(70,67,expression(X[t]~'~I(0)'), col = "grey50")
text(40,52,expression(X[t]~'~I(1)'), col = "grey50")
 
text(10,37,expression(H[0]~':'~Phi~'=0'))
text(40,37,expression(X[t]~'~I(0)'), col = "grey50")
text(10,22,expression(X[t]~'~I(1)'))
 
y=c(94,70)
x=c(11,11)
lines(x, y)
text(14,82,expression(H[0]))
 
y=c(64,40)
x=c(11,11)
lines(x, y)
text(14,52,expression(H[0]))
 
y=c(34,25)
x=c(11,11)
lines(x, y)
text(14,30,expression(H[0]))
 
y=c(97,97)
x=c(20,30)
lines(x, y, col = "grey50")
text(25,93,expression(H[1]), col = "grey50")
 
y=c(67,67)
x=c(20,30)
lines(x, y, col = "grey50")
text(25,63,expression(H[1]), col = "grey50")
 
y=c(37,37)
x=c(20,30)
lines(x, y, col = "grey50")
text(25,34,expression(H[1]), col = "grey50")
 
y=c(97,97)
x=c(50,60)
lines(x, y, col = "grey50")
text(55,93,expression(H[1]), col = "grey50")
 
y=c(67,67)
x=c(50,60)
lines(x, y, col = "grey50")
text(55,63,expression(H[1]), col = "grey50")
 
y=c(94,85)
x=c(40,40)
lines(x, y, col = "grey50")
text(43,90,expression(H[0]), col = "grey50")
 
y=c(64,55)
x=c(40,40)
lines(x, y, col = "grey50")
text(43,60,expression(H[0]), col = "grey50")

Le graphique représente tous les tests à réaliser pour faire un test ADF. En noir, c'est le chemin que je suis (sur mes données, mais ça aurait pu être un tout autre chemin sur d'autres données).

[A. D.]

Bonjour,

Je ne suis pas sûre d'avoir bien saisi ce que vous souhaitez faire, mais en ajoutant les commandes suivantes après les votres, on obtient peut-être quelque chose qui s'en rapproche (après ça reste du bidouillage ) :

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
x=c(0,0)
y=c(15,103)
lines(x,y)
 
x=c(80,80)
y=c(15,103)
lines(x,y)
 
x=c(0,80)
y=c(15,15)
lines(x,y)
 
x=c(0,80)
y=c(103,103)
lines(x,y)

Cordialement,

A.D.