minimo=c(2.135829,0,0.9998784,1.575533,0,0,7.065609) massimo=c(72.22222,29.32653,22.1049,34.73541,35.87713,1.477646,75.19857) sd_training_set=c(18.2467950,6.3022608,5.7965420,8.3413784,6.0488286,0.3525393,17.5918777) random=matrix(ncol=7,nrow=34) perf_random=matrix(ncol=7,nrow=100) for (z in 1:100) { for (i in 1:34) { for (j in 1:7) { random[i,j]=runif(1,min=minimo[j],max=massimo[j]) } } totale=rowSums(random) random1=random/totale*100 perf_random[z,] <- colMeans(abs(random1-sales))/sd_training_set } colMeans(perf_random) rowMeans(perf_random) mean(rowMeans(perf_random))