Combinar los dos gráficos es relativamente sencillo en este caso ya que estamos trabajando básicamente con los mismos datos y los mismos ejes. Podemos hacer lo siguiente:
S %>%
ggplot() +
geom_point(data=dataNU, aes(x=NU, y=Rend)) +
geom_smooth(data=dataNU,
formula = y ~ x + I(x^2),
method = "lm",
se=FALSE,
aes(colour="", x=NU, y=Rend) ) +
scale_colour_manual(name="NU", values=c("blue")) +
labs(x="NU",
y="Kg") +
geom_bar(data=S,
aes(x=NCP, y=Rend),
stat = "summary",
width=30,
alpha = c(0,0.3,0.3),
fill = "red",
fun.y = "mean") +
theme_minimal()
Resutado:

Detalle:
- Como habrás visto, tenemos que usar el
dataframe
completo ya que queremos graficar una barra para NP==360
y usando el conjunto subset(S, NCP == "0")
excluimos estos valores. Pero en el resto de los componentes usamos el subset
que ya armaste dataNU
- Agregamos un
geom_bar()
a la gráfica, la misma va a ser del tipo summary
y la función sumaria va a ser la media (fun.y = "mean"
)
- En tu ejemplo pareciera que no quieres una barra para
NP == 0
, si usamos el dataframe
completo, nos generará una barra para este grupo, pero hay un truco que es configurar la transparencia para que no se vea: alpha = c(0,0.3,0.3)
Lo siguiente es ver como dibujamos las rectas horizontales y verticales. Esto no es del todo trivial, al menos para mí, ya que en primer lugar tenemos que encontrar un punto dónde la recta horizontal cuyo y
(que se corresponde a la media de cada NU
) intersecciona con la curva de regresión, el tema es que justamente la regresión nos sirve para determinar un punto y
dado cualquier x
pero no a la inversa.
Seguramente hay formas más directas, pero me falta teoría, así que como no has mencionado que necesites precisión en este punto, vamos a usar una forma aproximativa para llegar a estos puntos. La idea es usar predict()
obtener todos los y
(Rend
) de cada punto del eje x
(NU
) de 0 a 360 en pasos de 0.1 (se pude ajustar este número para lograr mayor precisión), con estos números simplemente buscaremos aquel que más se aproxime a la media de cada barra:
# Calculamos la media de los valores de NCP 180 y 360
bar_means <- c(mean(S[S$NCP==180,3]), mean(S[S$NCP==360,3]))
# A lo cavernicola, generamos los y de la curva, luego
# buscamos el x mas cercano a la media de cada barra
model<-lm(Rend ~ (NU+ I(NU^2)), data = dataNU)
newx <-data.frame(NU=seq(0,360,0.1))
fitline = predict(model, newdata=newx)
xnew <- c(min(newx[fitline>bar_means[1],]),
min(newx[fitline>bar_means[2],])
)
Con esto puntos ya podemos definir los segmentos desde cada barra a la curva y de esta al eje x
.
# Generamos los datos para dibujar los segmentos verticales y horizontales
lineas_h <- data.frame(
pos = bar_means,
start = c(180, 360),
end = xnew
)
lineas_v <- data.frame(
pos = xnew,
start = bar_means,
end = c(0, 0)
)
Con estos datos, podemos usar geom_segment()
para dibujar estas líneas. Además podremos agregar las leyendas, pero tendremos una para cada barra (incluida la del 0) con colores distintos, por que efectivamente son grupos distintos, el código completo sería:
library(ggplot2)
library(magrittr)
NU <- c(0,180,90,0,0,0,90,180,0,0,180,90,0,0,0)
NCP <- c (180,0,0,0,360,360,0,0,180,0,0,0,180,0,360)
Rend<- c(14083,20856,19273,9655,15171,14021,14794,17578,11196,7201,13869,14114,8281,6041,12251)
S <- data.frame(NU=NU,NCP=NCP,Rend=Rend)
S
dataNU <- subset(S, NCP == "0")
# Calculamos la media de los valores de NCP 180 y 360
bar_means <- c(mean(S[S$NCP==180,3]), mean(S[S$NCP==360,3]))
# A lo cavernicola, generamos los y de la curva y
# buscamos el x mas cercano a la media de cada barra
model<-lm(Rend ~ (NU+ I(NU^2)), data = dataNU)
newx <-data.frame(NU=seq(0,360,0.1))
fitline = predict(model, newdata=newx)
xnew <- c(min(newx[fitline>bar_means[1],]),
min(newx[fitline>bar_means[2],])
)
# Generamos los datos para dibujar los segmentos verticales y horizontales
lineas_h <- data.frame(
pos = bar_means,
start = c(180, 360),
end = xnew
)
lineas_v <- data.frame(
pos = xnew,
start = bar_means,
end = c(0, 0)
)
# Plot sobre el data.frame completo
S %>%
ggplot() +
geom_point(data=dataNU, aes(x=NU, y=Rend)) +
geom_smooth(data=dataNU,
formula = y ~ x + I(x^2),
method = "lm",
se=FALSE,
aes(colour="", x=NU, y=Rend)
) +
scale_colour_manual(name="NU", values=c("blue")) +
labs(x="NU",
y="Kg") +
geom_bar(data=S,
aes(x=NCP, y=Rend, fill = factor(NCP)),
stat = "summary",
width=30,
alpha = c(0,0.3,0.3),
fun.y = "mean") +
geom_segment(data=lineas_h,aes(x=start, y=pos, xend=end, yend=pos), linetype=2) +
geom_segment(data=lineas_v,aes(y=start, x=pos, yend=end, xend=pos), linetype=2) +
labs(fill="NCP") +
theme_minimal()
Y el resultado final

El resultado: