Al ver la pregunta, si se omitiera las escalas y los valores de los ejes Y -no sus nombres-, pareciera que se quiere encontrar una especie de relación o correlación entre las dos series.
Si bien puede ser válido, hay que tener presente que correlación no siempre va de la mano con la causa. Esto le pasó al proyecto google flu, incluso hay un artículo en la revista science que habla al respecto The Parable of Google Flu: Traps in Big Data Analysis
Sobre lo anterior, también es posible acceder a muchos ejemplos reales en la página spurious correlations
Lo anterior quiere decir que, una variable de interés sólo puede ser modela a partir de otras variables subyacentes, siempre y cuando estas últimas sean independientes y la primera dependiente de las mismas.
Expuesto lo anterior existe una biblioteca latticeExtra
, se generan los datos:
dias <-100
capinst<-rnorm(dias, mean=80, sd=10)
infla<-rweibull(dias, shape=1, scale=0.75)
index <- 1:dias
capinfl<- data.frame(
capinst=capinst,
infla=infla,
index = index
)
Se genera el gráfico con las dos series
library(latticeExtra)
# se puede mostrar las dos variables en el mismo tiempo
xyplot(capinst + infla ~ index, capinfl, type="l")
# se agrega dos ejes Y
# se construye cada serie por separado
obj1 <- xyplot(capinst ~ index, capinfl, type = "l" , lwd=2, ylab="capinst", xlab="xDías")
obj2 <- xyplot(infla ~ index, capinfl, type = "l", lwd=2, ylab="inflacion")
# --> se realiza la grafica con el segundo eje Y
doubleYScale(obj1, obj2, add.ylab2 = TRUE)
Como gráfico interno de trabajo lo veo adecuado pero si va en un informe final o digamos mediante una implementación con Shiny para visualizarlo de forma dinámica en una Web App tomaría las palabras @mpaladino, podría desorientar al usuario que lo va leer.
No obstante hay un caso que encontré dónde si es válido aplicar un gráfico como el que buscan, -humedad relativa para termodinámica- está incluso el código para lograrlo con ggplot2
, sólo tendrían que hallar el factor de escala (en el ejemplo está con 5) que se adecue a sus necesidades.
ggplot2 version 2.2.0 - Demonstration of dual y-axes using sec.axis
Aquí una sugerencia de como hallar ese factor de escala
library(ggplot2)
# Mostrando la serie capinst
p <- ggplot(capinfl, aes(x = index))
p <- p + geom_line(aes(y = capinst), colour="blue")
# como capinst tiene una escala mayor que infla se
# requiere ajustar los datos
sf<-max(capinst)/max(infla)
# Se agrega infla a Y multiplicada por el factor
p <- p + geom_line(aes(y = infla*sf), colour="red")
p <- p + scale_y_continuous(sec.axis = sec_axis(~.*sf, name = "Inflación Escalada"))
p <- p + labs(y = "Capinst",
x = "xFechas",
colour = "Parameter")
p <- p + theme(legend.position = c(0.8, 0.9))
p
Pero si observan el Eje Y secundario (Inflación Escalada) apareceré la situación que se reflejaría en el Shiny
si nos situamos en cada punto de la serie roja (inflación)