# 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)[ind])
+print(data[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="")