require(e1071)#pour coef
require(sp)
require(geoR)
require(gstat)
#importer les données
test_mar <- read.table("data.dat",sep=",",header = TRUE )

test_mar$X<-as.numeric(test_mar$X)
test_mar$Y<-as.numeric(test_mar$Y)
#transformer les donnes en classe de donnees spatiale
coordinates(test_mar)<- ~X+Y
class(test_mar)
test_mar$marbre<-as.numeric(test_mar$marbre)

bubble(test_mar,"marbre",fill=TRUE,maxsize = 2.5,scales = list(draw=TRUE),
       main="Repartition des points échantillonnés")

vr <- variogram(bre ~ coordinates(marbre.all), loc = marbre.all,cutoff = 700)
vr
plot(vr,pch=11,type="p")
#ajustement du vario
vm <- vgm(0.0475, "Sph", 290, 0.0002)
p<-plot(vr, plot.numbers = T, pch = 11, col = "darkblue", model = vm)


kr <- krige(test_mar$marbre~1, test_mar ,datas.Grid, vm)

> kr <- krige(test_mar$marbre~1, test_mar ,gr, vm)
Error in cbind(newdata@bbox[cbind(c(1, 1, 1, 1), c(1, 2, 1, 1))], newdata@bbox[cbind(c(2,  : 
  trying to get slot "bbox" from an object (class "data.frame") that is not an S4 object