Esther Quiñones Luna
MQECONOMICS
PAQUETERIA
library(tidyverse)
library(dplyr)
library(stringr)
library(lubridate)
library(ggplot2)
library(plotly)
library(sp)
library(RCurl)
CÓDIGO
#Indice de precios al productor
inpp <- read.csv("C:\\INPP.csv", encoding = "UTF-8")
inpp <- rename(inpp, Fecha=X.U.FEFF.Fecha)
inpp1 <-ts(inpp,freq=12,start=c(2010,6))
dinpp3 <- diff(inpp1, lag=12)
dinpp3 <- as.data.frame(dinpp3)
linpp2 <- lag(inpp,12)
linpp2 <-linpp2[-c(1:12),]
dinpp3 <- select (dinpp3, -Fecha)
linpp2 <- select(linpp2, -Fecha)
dinpp3 <- as.matrix (dinpp3)
linpp2 <- as.numeric(linpp2)
linpp2 <- as.matrix(linpp2)
dlinpp4 <- (dinpp3/linpp2)*100
dlinpp4 <- as.data.frame(dlinpp4)
n <- nrow(inpp)
Fecha <- inpp$Fecha
Fecha <- Fecha [13:n]
dlinpp4$Fecha <-Fecha
d1 <- "01/06/2011"
d2 <- "01/02/2022"
dia1=as.Date(d1,format="%d/%m/%Y")
dia2=as.Date(d2,format="%d/%m/%Y")
fecha<-seq(dia1,dia2,by="month")
dlinpp4$fecha <-fecha
#Cambiar de vuelta a la configuración inicial
n <-nrow(dlinpp4)
dlinpp4$N <- 1:n
dlinpp4<- data.frame(dlinpp4)
ginpp1 <- select(dlinpp4, Fecha, fecha, N, X23.Construcción,
X313.Fabricación.de.insumos.textiles.y.acabado.de.textiles,
X314.Fabricación.de.productos.textiles..excepto.prendas.de.vestir,
X315.Fabricación.de.prendas.de.vestir,
X316.Curtido.y.acabado.de.cuero.y.piel,
X336.Fabricación.de.equipo.de.transporte,
X71.Serv.esparcimiento.culturales.y.deportivos..y.otros,
X72.Serv.de.alojamiento.temporal.y.de.preparación.de.alimentos.y.bebidas)
ginpp1 <- rename(ginpp1, Construcción=X23.Construcción,
Fabricación.de.insumos.textiles.y.acabado.de.textiles=X313.Fabricación.de.insumos.textiles.y.acabado.de.textiles,
Fabricación.de.productos.textiles..excepto.prendas.de.vestir=X314.Fabricación.de.productos.textiles..excepto.prendas.de.vestir,
Fabricación.de.prendas.de.vestir=X315.Fabricación.de.prendas.de.vestir,
Curtido.y.acabado.de.cuero.y.piel=X316.Curtido.y.acabado.de.cuero.y.piel,
Fabricación.de.equipo.de.transporte= X336.Fabricación.de.equipo.de.transporte,
Serv.esparcimiento.culturales.y.deportivos..y.otros=X71.Serv.esparcimiento.culturales.y.deportivos..y.otros,
Serv.de.alojamiento.temporal.y.de.preparación.de.alimentos.y.bebidas= X72.Serv.de.alojamiento.temporal.y.de.preparación.de.alimentos.y.bebidas)
#Formato xts
library(xts)
library(quantmod)
ginpp1 <-as.data.frame(ginpp1)
rownames(ginpp1) <- as.Date(ginpp1$fecha)
ginpp1 <-as.matrix(ginpp1)
class(ginpp1)
ginpp1.xts <- xts(x=ginpp1, order.by=fecha, frequency = 1)
#Gráfico 1
library(tidyquant)
library(plotly)
inppg2 <- data.frame(ginpp1.xts$Construcción,
ginpp1.xts$Fabricación.de.insumos.textiles.y.acabado.de.textiles,
ginpp1.xts$Fabricación.de.productos.textiles..excepto.prendas.de.vestir,
ginpp1.xts$Fabricación.de.prendas.de.vestir,
ginpp1.xts$Curtido.y.acabado.de.cuero.y.piel,
ginpp1.xts$Fabricación.de.equipo.de.transporte,
ginpp1.xts$Serv.esparcimiento.culturales.y.deportivos..y.otros,
ginpp1.xts$Serv.de.alojamiento.temporal.y.de.preparación.de.alimentos.y.bebidas
)
inppg2 <- data.frame(inppg2,rownames(inppg2))
fecha <- inppg2$rownames.inppg2.
fig <- plot_ly(inppg2,type = 'scatter', mode = 'lines')%>%
add_trace(x = ~fecha, y = ~Construcción, name = 'Construcción')%>%
add_trace(x = ~fecha, y = ~Fabricación.de.insumos.textiles.y.acabado.de.textiles, name = 'Fabricación de insumos textiles y acabado de textiles')%>%
add_trace(x = ~fecha, y = ~Fabricación.de.productos.textiles..excepto.prendas.de.vestir, name = 'Fabricación de productos textiles excepto prendas de vestir') %>%
add_trace(x = ~fecha, y = ~Fabricación.de.prendas.de.vestir, name = 'Fabricación de prendas de vestir') %>%
add_trace(x = ~fecha, y = ~Construcción, name = 'Construcción')%>%
add_trace(x = ~fecha, y = ~Curtido.y.acabado.de.cuero.y.piel, name = 'Curtido y acabado de cuero y piel')%>%
add_trace(x = ~fecha, y = ~Fabricación.de.equipo.de.transporte, name = 'Fabricación de equipo de transporte') %>%
add_trace(x = ~fecha, y = ~Serv.esparcimiento.culturales.y.deportivos..y.otros, name = 'Servicios de esparcimiento culturales, deportivos y otros') %>%
add_trace(x = ~fecha, y = ~Serv.de.alojamiento.temporal.y.de.preparación.de.alimentos.y.bebidas, name = 'Servicios de alojamiento temporal y de preparación de alimentos y bebidas')%>%
layout(title = 'Indice Nacional de Precios Productor',legend=list(title=list(text='INPP')),
xaxis = list(rangeslider = list(visible = T),
rangeselector=list(
buttons=list(
list(count=6, label="6m", step="month", stepmode="backward"),
list(count=1, label="YTD", step="year", stepmode="todate"),
list(count=1, label="1y", step="year", stepmode="backward"),
list(step="all")
))))
fig <- fig %>%
layout(
xaxis = list(title= 'Fecha',
zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
yaxis = list(title= 'Variación anual',
zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
plot_bgcolor='#e5ecf6')
fig
RESULTADO
Datos utilizados
Comments