なぜ私がこれを手に入れたのか誰か教えてもらえます"Error in (function (classes, fdef, mtable)"か?
R の例で、2 番目の「グリッド」データベースについて何かが見えます。これを作成する必要がありますか? これまでのところ、ターゲット変数、予測子、および座標はすべて同じデータベースにあります。
または、座標系 (CRS) を指定する必要がありますか?
ありがとうございました。
library(gstat)
library(maptools)
clean3145 = read.csv(file="clean3145.csv")
coordinates(clean3145)=~UTMEM+UTMNM
summary(clean3145)
#Variogram of LM residuals
logAs<-lm(logAs1 ~ Elev_m + Basin, clean3145)
residuals(logAs)
variogram(residuals(logAs)~UTMEM+UTMNM, clean3145)
plot(variogram(residuals(logAs)~UTMEM+UTMNM, clean3145), plot.nu=T, pch="+")
#Fit variogram
As.rev=variogram(residuals(logAs)~UTMEM+UTMNM, clean3145)
as.rvgm=fit.variogram(As.rev, vgm(nugget=0.25, model="Exp", range=15000, sill=0.6))
plot(As.rev, as.rvgm, plot.nu=T)
str(as.rvgm)
#Regression Kriging (rk)
library(sp)
AS.rk=krige(logAs, clean3145, as.rvgm)