web-dev-qa-db-fra.com

Ajouter une ligne verticale dans la parcelle ggplot

Je trace un graphique en utilisant le code suivant:

library (ggplot2)

png (filename = "graph.png")
stats <- read.table("processed-r.dat", header=T, sep=",")
attach (stats)
stats <- stats[order(best), ]
sp <- stats$A / stats$B
index <- seq (1, sum (sp >= 1.0))
stats <- data.frame (x=index, y=sp[sp>=1.0])
ggplot (data=stats, aes (x=x, y=y, group=1)) + geom_line()
dev.off ()

enter image description here

1 - Comment peut-on ajouter une ligne verticale dans le tracé qui se croise à une valeur particulière de y (par exemple 2)?

2 - Comment faire en sorte que l’axe des y commence à 0,5 au lieu de 1?

rggplot2
43
Shahzad

Vous pouvez ajouter une ligne verticale avec geom_vline(). Dans ton cas:

+ geom_vline(xintercept=2)

Si vous voulez également voir le nombre 0.5 sur votre axe y, ajoutez scale_y_continuous() et définissez limits= Et breaks=

+ scale_y_continuous(breaks=c(0.5,1,2,3,4,5),limits=c(0.5,6))
66
Didzis Elferts

Concernant la première question:

Cette réponse suppose que la valeur de Y que vous désirez est spécifiquement dans votre ensemble de données. Tout d’abord, créons un exemple reproductible car je ne peux pas accéder à votre ensemble de données:

set.seed(9999)
stats <- data.frame(y = sort(rbeta(250, 1, 10)*10 ,decreasing = TRUE), x = 1:250)
ggplot(data=stats, aes (x=x, y=y, group=1)) + geom_line()

initial plot

Ce que vous devez faire est d'utiliser la colonne y de votre cadre de données pour rechercher la valeur spécifique. Essentiellement, vous devrez utiliser

ggplot(data=stats, aes (x=x, y=y, group=1)) + geom_line() + 
    geom_vline(xintercept = stats[stats$y == 2, "x"])

En utilisant les données que j'ai générées ci-dessus, voici un exemple. Comme mon bloc de données ne contient probablement pas la valeur exacte 2, Je vais utiliser la fonction trunc pour le rechercher:

stats[trunc(stats$y) == 2, ]

#           y  x
# 9  2.972736  9
# 10 2.941141 10
# 11 2.865942 11
# 12 2.746600 12
# 13 2.741729 13
# 14 2.693501 14
# 15 2.680031 15
# 16 2.648504 16
# 17 2.417008 17
# 18 2.404882 18
# 19 2.370218 19
# 20 2.336434 20
# 21 2.303528 21
# 22 2.301500 22
# 23 2.272696 23
# 24 2.191114 24
# 25 2.136638 25
# 26 2.067315 26

Maintenant nous savons où sont toutes les valeurs de 2. Comme ce graphique est en baisse, nous l'inverserons, la valeur la plus proche de 2 sera donc au début:

rev(stats[trunc(stats$y) == 2, 1])

#           y  x
# 26 2.067315 26

Et nous pouvons utiliser cette valeur pour spécifier l'endroit où l'interception x devrait être:

ggplot(data=stats, aes (x=x, y=y, group=1)) + geom_line() + 
    geom_vline(xintercept = rev(stats[trunc(stats$y) == 2, "x"])[1])

plot with vertical line

J'espère que ça t'as aidé!

21
ZNK