titolo<-read.csv("c:/dati/intesa.csv",header=T,dec=".",sep=",") attributes(titolo) nrow(titolo) ncol(titolo) names(titolo)[1]<-"Date" row.names(titolo)[1:10] titolot<-t(titolo) length(titolo) p.uff<-sum(titolo[,5])/sum(titolo[,4]) p.uff1 p.uff1<-sum(titolo$CTVL)/sum(titolo$QUANT) p.min<-min(titolo$PREZ) p.min pp.min<-which(titolo$PREZ==p.min) pp.min titolo$ORA[pp.min] p.max<-max(titolo$PREZ) p.max pp.max<-which(titolo$PREZ==p.max) pp.max titolo$ORA[pp.max] K<-p.max-p.min K Kr<-K/p.uff*100 Kr sigma<-sqrt(var(titolo$PREZ)) sigma CV<-sigma/p.uff*100 CV rif<-floor(0.10*nrow(titolo)) titolo.rif<-titolo[(nrow(titolo)-rif):nrow(titolo),] p.rif<-sum(titolo.rif$CTVL)/sum(titolo.rif$QUANT) p.rif p.uff infra.r<-function(nome) { titolo<-read.csv(paste("c:/dati/",nome,".csv",sep=""),header=T,dec=".",sep=",") p.uff<-sum(titolo[,5])/sum(titolo[,4]) p.min<-min(titolo$PREZ) p.max<-max(titolo$PREZ) rif<-floor(0.10*nrow(titolo)) titolo.rif<-titolo[(nrow(titolo)-rif):nrow(titolo),] p.rif<-sum(titolo.rif$CTVL)/sum(titolo.rif$QUANT) K<-p.max-p.min Kr<-K/p.uff*100 sigma<-sqrt(var(titolo$PREZ)) CV<-sigma/p.uff*100 out<-list(p.uff=p.uff,p.max=p.max,p.min=p.min, p.rif=p.rif,Kr=Kr,CV=CV) } out.intesa<-infra.r("intesa") out.intesa out.medio<-infra.r("medio") out.medio