confr<-function(nome) { titolo<-read.csv(paste("c:/msmf/",nome,".csv",sep=""),header=T,dec=".",sep=",") # Calcolo del prezzo ufficiale della seduta p.uff<-sum(titolo$CTVL)/sum(titolo$QUANT) # Individuazione del prezzo minimo p.min<-min(titolo$PREZ) rp.min<-which(titolo[,3]==p.min) # Individuazione del prezzo massimo p.max<-max(titolo$PREZ) rp.max<-which(titolo[,3]==p.max) # Calcolo del prezzo di riferimento titolo.rif<-titolo[floor(0.9*nrow(titolo)):nrow(titolo),] p.rif<-sum(titolo.rif$CTVL)/sum(titolo.rif$QUANT) # Volatilità c.var<-sqrt(var(titolo$PREZ))/p.uff*100 sqm<-(p.max-p.min)/p.uff*100 # Grafico prezzi ts.plot(titolo$PREZ/titolo$PREZ[1]*100,xlab="Time",ylab="P_t",main=nome) out<-list(p.uff=p.uff,p.max=p.max,p.min=p.min,p.rif=p.rif,c.var=c.var,sqm=sqm) } # Finestra grafica par(mfrow=c(2,1)) intesa<-confr("intinfr") medio<-confr("mediobanca")