# Lettura dei dati e ordinamento
titolo<-read.csv("c:/dati/tenaris.csv",dec=".",sep=",",header=T)
durata<-nrow(titolo)
ordine<-as.matrix(seq(durata,1,-1))
titolo<-cbind(ordine,titolo)
titolo<-titolo[sort.list(titolo[,1]),]
close.p<-as.matrix(titolo$Close)
close.ts<-ts(close.p,start=c(2002,12,17),frequency=250,names="Tenaris")
volumi<-(titolo$Volume)

# Finestra grafica
par(mfrow=c(2,1))
ts.plot(close.ts)
barplot(volumi)
par(mfrow=c(1,1))

# Calcolo rendimenti giornalieri (RA_t, R_t, r_t)
close.ra<-diff(close.ts)
close.Rt<-((close.ts/lag(close.ts,-1))-1)*100
close.rt<-(diff(log(close.ts)))*100

# Grafico prezzi e rendimenti
par(mfrow=c(2,1))
ts.plot(close.ts,ylab="Prezzi",xlab="Tempo",main="Tenaris")
ts.plot(close.rt,ylab="Rendimenti log",xlab="Tempo")
par(mfrow=c(1,1))
names(titolo)
#Sistemazione dell'output in tabella
tabella<-cbind(titolo[(durata-9):durata,2],
close.ts[(durata-9):durata],
close.ra[(durata-10):(durata-1)],
close.Rt[(durata-10):(durata-1)],
close.rt[(durata-10):(durata-1)])
tabella<-as.data.frame(tabella)
names(tabella)<-c("Data","P_t","RA_t","R_t","r_t")
write.table(tabella,"c:/dati/output.csv",sep=",",append=F)