setwd("C:/MASTER-Profesores/Germán/08-09/Apuntes")


###########################################################################
#######     TEMA 1: ANÁLISIS DESCRIPTIVO DE UNA SERIE DE TIEMPO    ########
###########################################################################

###################### ALGUNAS SERIES  #######################      
pib <- scan("pib.txt")
t <- 1960:2006
plot(t, pib, type="o", xlab="Año", ylab="PIB (en miles de euros)")

temp <- scan("temperatura.txt")
t <- 1:180
plot(t, temp, type="o", xlab="Mes", ylab="Temperatura (en ºC)")

tur <- scan("turismo.txt")
t <- 1:116
plot(t, tur, type="o", xlab="Mes", ylab="Turismo (en miles)")

cons <- scan("consumo_electricidad.txt")
t <- 1:396
plot(t, cons, type="o", xlab="Mes", ylab="Consumo de electricidad")


###################### MEDIAS MÓVILES ######################
pib <- scan("pib.txt")
t <- 1960:2006

ajuste5 <- filter(pib, sides=2, rep(1/5,5))
ajuste11 <- filter(pib, sides=2, rep(1/11,11))
residuos5 <- pib-ajuste5
residuos11 <- pib-ajuste11

plot(t, pib, type="o", xlab="Año", ylab="PIB")
lines(t, ajuste5, col="blue")
lines(t, ajuste11, col="red", lty="dashed")

plot(t, residuos11, type="o", col="red", lty="dashed", xlab="Año", ylab="Residuos")
lines(t,residuos5, type="o", col="blue")
abline(0,0)

par(mfrow=c(3,1))
plot(t, pib, type="l", xlab="Año", ylab="PIB")
plot(t, ajuste11, type="l", xlab="Año", ylab="Tendencia")
plot(t, residuos11, type="l", xlab="Año", ylab="Residuos")

plot(t, pib, type="o", xlab="Año", ylab="PIB")
lines(t,ajuste11, col="red", lty="dashed")

##########################################################
tur <- scan("turismo.txt")
tur.ts <- ts(tur, start=c(1995, 1), frequency=12)

descomposicion <- decompose(tur.ts, type="additive")

plot(descomposicion)

plot(tur.ts, type="o", xlab="Año", ylab="Turismo")
lines(descomposicion$trend + descomposicion$seasonal, col="red", lty="dashed")

##########################################################
cons <- scan("consumo_electricidad.txt")
cons.ts <- ts(cons, start=c(1972, 1), frequency=12)

descomposicion <- decompose(cons.ts, type="multiplicative")

plot(descomposicion)

plot(cons.ts, type="o", xlab="Año", ylab="Consumo de electricidad")
lines(descomposicion$trend*descomposicion$seasonal, col="red", lty="dashed")

##########################################################
temp <- scan("temperatura.txt")
temp.ts <- ts(temp, start=c(1988, 1), frequency=12)

descomposicion <- decompose(temp.ts, type="additive")

plot(descomposicion)

plot(temp.ts, type="o", xlab="Mes", ylab="Temperatura")
lines(descomposicion$trend + descomposicion$seasonal, col="red", lty="dashed")