-# pdf("Erreur_Epandage_PMjour.pdf")
-ymin = min(na.omit(ErrPrev), Err24)
-ymin = min(ymin, Kvois)
-ymax = max(na.omit(ErrPrev), Err24)
-ymax = max(ymax, Kvois)
-plot(5:24,Err24, type = "l", lwd = 2, cex.axis=1.4, cex.lab = 1.5,
- ylab = "MAE",xlab = "Heures de prévision", ylim= c(ymin,ymax))
-lines(5:24,ErrPrev, lwd=2, lty = 2)
-legend("topright", legend=c("Erreur 24h", "Erreur prévision"), lty = c(1,2), lwd=2, cex=1.5)
-points(5:24, Kvois, cex=1.8, pch = 19)
-# xx = dev.off()
-
-#length(D)
+ list(serie=dataj, prev=JourMoy, err=erreurPrev, line=(nl+ij), neighbs=ind, dates=data[ind,1])
+# erreur24 = mean(abs(dataj[25:48] - JourMoy[25:48]))
+# #png(NomFile)
+# matplot(t(data[ind, 1:48]), type = "l", lwd=1.4, lty=1, col=1:length(ind),
+# cex.axis=1.4, cex.main = 1.7, cex.lab=1.5,
+# xlab="Heures locales", ylab=paste0("PM10 - Erreurs = ",round(erreur24,1),
+# " / ",round(erreurPrev,1)), main=Titre)
+# legend("top",rownames(data)[ind],ncol=2,lwd=1.4, lty=1, col=1:length(ind))
+# lines(1:48, dataj, lwd=2.5)
+# lines(JourMoy, lty = 2, lwd=2)
+# abline(v=c(24.5,H+0.5), lty = 2, lwd=1.2)
+# #xx=dev.off()
+#
+# Err24 = c(Err24, erreur24)
+# ErrPrev = c(ErrPrev, erreurPrev)
+# ResDates = cbind(ResDates, rownames(data)[ind])
+#}
+
+#rownames(ResDates) = 1:10
+#
+#Kvois = NULL
+#for (Col in ncol(ResDates):1)
+#{
+# K = 0
+# for (I in 1:10)
+# {
+# for (J in 1:10)
+# {
+# if (ResDates[I,1] == ResDates[J,Col])
+# K = K +1
+# }
+# }
+# Kvois = c(Kvois, K)
+#}
+#
+## pdf("Erreur_Epandage_PMjour.pdf")
+#ymin = min(na.omit(ErrPrev), Err24)
+#ymin = min(ymin, Kvois)
+#ymax = max(na.omit(ErrPrev), Err24)
+#ymax = max(ymax, Kvois)
+#plot(5:24,Err24, type = "l", lwd = 2, cex.axis=1.4, cex.lab = 1.5,
+# ylab = "MAE",xlab = "Heures de prévision", ylim= c(ymin,ymax))
+#lines(5:24,ErrPrev, lwd=2, lty = 2)
+#legend("topright", legend=c("Erreur 24h", "Erreur prévision"), lty = c(1,2), lwd=2, cex=1.5)
+#points(5:24, Kvois, cex=1.8, pch = 19)
+## xx = dev.off()
+#
+##length(D)
+}