ブロック クリギングに問題があります。まず、グリッド (5000*5000 m) を作成します。お互いに非常に近い駅を除外します。次のステップでは、バリオグラムと fit.variogram を作成します。投影と座標を割り当てます。その後、通常のクリギングに対してクリゲ関数が機能します。
x_2 <-as.integer(range(data$x))
y_2 <-as.integer(range(data$y))
grid_1 <-expand.grid(x=seq(from=x_2[1],to=x_2[2], by=5000),y=seq(y_2[1],to=y_2[2],by=5000))
coordinates(grid_1) <- ~x+y #Koordinaten zuweisen, Variablennamen festlegen, auf die man sich bezieht)
gridded(grid_1) =TRUE
grid.stats <- data.frame(data)[,1:2]
coordinates(grid.stats) <- ~x+y
ind.zd <- zerodist(grid.stats)
data.new <- data[-ind.zd,]
p.sill <- var(log(data.new$variable+0.1)) #0.4005127
nugget <- var(log(data.new$variable+0.1))/2 #0.2002564
vgm_data <- variogram(log(variable+0.1)~1, ~x+y, data.new)
m.fit <- fit.variogram(vgm_data, vgm(0.4, "Sph", 1000000, 0.2))
coordinates(data.new) <- ~x+y
krig_1 <- krige(log(variable+0.1)~x+y, data.new, grid_1, nmax=7, nmin=4, m=m.fit)
krig_block <- krige(log(variable+0.1)~1, data.new, grid_1, m=m.fit, block =c(5000/5000))
ブロック クリギングに「ブロック」を追加しようとしましたが、結果が得られません。私の問題を知っている人はいますか？事前に感謝します。敬具