# w = 1/(D[ind]^2)
# w = w/sum(w)
# W = w %o% rep(1,48)
+
+#print("Voisins + ij")
+#print(sort(D)[1:nbvois])
+#print(dateJPrev)
+#print(rownames(data)[sort(ind)])
+#print(data[sort(ind), 1:48])
+
JourMoy = apply(data[ind, 1:48], 2, mean)
#JourMoy = apply(W*data[ind, 1:48], 2, sum)
NomFile = paste("Voisins_Epandage_PMjour_Hc_",Hc,".png",sep="")