Partage
  • Partager sur Facebook
  • Partager sur Twitter

package R

    3 juillet 2021 à 22:51:39

    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 


    • Partager sur Facebook
    • Partager sur Twitter

    package R

    × Après avoir cliqué sur "Répondre" vous serez invité à vous connecter pour que votre message soit publié.
    • Editeur
    • Markdown