r - ggplot2 - resaltar los puntos seleccionados y el comportamiento extraño

Question

Quiero resaltar los puntos seleccionados y encontré un comportamiento extraño. Primero algunos datos ficticios:

a <- 1:50 
b <- rnorm(50) 
mydata <- data.frame(a=a,b=b) 
ggplot(mydata,aes(x=a,y=b)) + geom_point() 

Esto funciona correctamente. Ahora, para poner de relieve algunos puntos, añado otra capa geom_point:

ggplot(mydata[20:40,],aes(x=a,y=b)) + 
    geom_point() + 
    geom_point(aes(x=a[c(10,12,13)],y=b[c(10,12,13)]),colour="red") 

Tenga en cuenta que estoy mostrando sólo una gama limitada de los datos ([20:40]). Ahora viene el comportamiento extraño:

ggplot(mydata[10:40,],aes(x=a,y=b)) + 
    geom_point() + 
    geom_point(aes(x=a[c(10,12,13)],y=b[c(10,12,13)]),colour="red") 

Cambiar el tamaño del rango seleccionado, me sale un error, más o menos traducido del alemán: Error...: Arguments implying different number of rows. Extrañamente, esto varía con el rango seleccionado. [23:40] funcionará, [22:40] no lo hará.

---

El error en Inglés es:

Error in data.frame(x = c(19L, 21L, 22L), y = c(0.28198, -0.6215, : 
    arguments imply differing number of rows: 3, 31 

Answer 1

Si los datos es diferente entre las diferentes capas, a continuación, es necesario especificar los nuevos datos para cada capa.

Esto se hace con el argumento data=... para cada geom que necesita de datos diferente:

set.seed(1) 
mydata <- data.frame(a=1:50, b=rnorm(50)) 
ggplot(mydata,aes(x=a,y=b)) + 
    geom_point(colour="blue") + 
    geom_point(data=mydata[10:13, ], aes(x=a, y=b), colour="red", size=5) 

Answer 2

Otra opción de añadir las condiciones para ambos atributos, color y tamaño, dentro geom_point. Luego controlamos manualmente aquellos que usan scale_colour_manual y scale_size_manual respectivamente.

set.seed(1) 
mydata <- data.frame(a = 1:50, b = rnorm(50)) 
ggplot(mydata) + 
    geom_point(aes(x = a, y = b, colour = a > 9 & a < 14, size = a > 9 & a < 14)) + 
    scale_colour_manual(values = c("blue", "red")) + 
    scale_size_manual(values =c(1, 4))+ 
    theme(legend.position = "none")