2

データセットから適応したこのコードを使用しようとしていますmeuse

data<-list(var1,var2,x,y)

coordinates(data)=~x+y

grid = spsample(data, type = "regular", cellsize = c(0.05,0.05))

vt <- variogram(var1 ~ var2,data=data)

vt.fit <- fit.variogram(vt, vgm(0.2, "Sph", 800, 0.05))

gstatobj <- gstat(id = 'var1', formula = var1 ~ var2, model=vt.fit, set = list(gls=1))

私の目標は、のようなグリッドを作成することですmeuse.grid。しかし、うまくいきcoordinatesません... list は適切なコマンドではありません。

何を使おうかな？グリッドの作成に使用している方法は正しいですか?

4

1 に答える 1

1

次の再現可能な例は、jlhoward のコメントが正しく、Darko の応答が間違っていることを示しています。

library(gstat)
var1 = 1:3; var2 = 1:3; x = 1:3; y = 1:3
data<-list(var1,var2,x,y)
coordinates(data) = ~x+y
Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘coordinates<-’ for signature ‘"list"’
data<-data.frame(var1,var2,x,y)
coordinates(data) = ~x+y
class(data)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

これをもう一度行うと混乱するかもしれません。

coordinates(data) = ~x+y
Error in `coordinates<-`(`*tmp*`, value = ~x + y) : 
  setting coordinates cannot be done on Spatial objects, where they have already been set

ただし、既存の（そして正しい）ものはそのままにしますdata。

于 2015-02-28T20:07:31.397 に答える