Bonsoir svp je rencontre un problème avec le package gstat du logiciel ,j'aimerai que vous m'aidiez ...merci
je joins à ceci le code et l'erreur:
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
package R
× Après avoir cliqué sur "Répondre" vous serez invité à vous connecter pour que votre message soit publié.
× Attention, ce sujet est très ancien. Le déterrer n'est pas forcément approprié. Nous te conseillons de créer un nouveau sujet pour poser ta question.