# INIZIALIZZAZIONE
library(tseries)
options(digits=6) # APPROSSIMAZIONE ALLA 6° CIFRA DECIMALE

# INPUT E PREPARAZIONE DATI
dati<-read.table("http://www.climatemonitor.it/Download/R/UAH.csv",header=FALSE,sep="", dec=".",col.names=c("Anno","Mese","Anomalia"))
UAH_mensile<-data.frame(dati)
attach(UAH_mensile)
c<-nrow(UAH_mensile)
ultimo_mese<-UAH_mensile$Mese[c]
ultimo_anno<-UAH_mensile$Anno[c]
ultimo_valore<-Anomalia[c]
UAH<-ts(Anomalia,start=c(1978,12),frequency=12)

# GRAFICO PRINCIPALE
plot(UAH,type="s",col="darkgrey",xlab="",ylab = "°C - Anomalia termica",xlim=c(1979, 2010), ylim=c(-0.6, 0.80))
lines(UAH,type="h",col="lightgrey") # GRAFICO DI FONDO
abline(h=0,col="darkgrey")

# TITOLI
mtext("Climate Monitor - www.climatemonitor.it",side=1,line=2,adj=1,cex=0.7)
mtext("Dati: http://vortex.nsstc.uah.edu/public/msu/t2lt/tltglhmam_5.2",side=1,line=3,adj=1,cex=0.7)
title(main="Anomalie termiche - Dati UAH",cex=0.9)

# REGRESSIONE
fit<-lm(UAH~Anno)
a<-coef(fit)[1]
b<-coef(fit)[2]
Anno1<-min(Anno)
Anno2<-max(Anno)
y1<-a+b*Anno1
y2<-a+b*Anno2
x_coord<-c(Anno1,Anno2)
y_coord<-c(y1,y2)
lines(x_coord,y_coord,type="l",col="red")

# LOWESS
lines(lowess(UAH,f=.1),col="blue",lwd=2)