r - バリオグラムとバリオグ関数の比較

Question

私は (おそらく間違って) 最も簡単なケースでvariogは、geoRパッケージ内とパッケージvariogram内の出力はsp同じであると想定しました。

私はこのデータセットを持っています:

head(final)
lat     lon  elev seadist  tradist samples rssi
1 60.1577 24.9111 2.392     125 15.21606     200  -58
2 60.1557 24.9214 3.195     116 15.81549     200  -55
3 60.1653 24.9221 4.604     387 15.72119     200  -70
4 60.1667 24.9165 7.355     205 15.39796     200  -62
5 60.1637 24.9166 3.648     252 15.43457     200  -73
6 60.1530 24.9258 2.733      65 16.10631     200  -57

それは (私が推測する) 投影されていないデータでできているので、それらを投影します

#data projection
#convert to sp object:
coordinates(final) <- ~ lon + lat #longitude first
library(rgdal)
proj4string(final) =  "+proj=longlat +datum=WGS84"
UTM <- spTransform(final, CRS=CRS("+proj=utm +zone=35V+north+ellps=WGS84+datum=WGS84"))

gstatライブラリに従って傾向のないバリオグラムを作成します

var.notrend.sp<-variogram(rssi~1, UTM)
plot(var.notrend.sp)

geoR私が行くのと同じ出力を得ようとしています

UTM1<-as.data.frame(UTM)
UTM1<-cbind(UTM1[,6:7], UTM1[,1:5])
UTM1
coords<-UTM1[,1:2]
coords
var.notrend.geoR <- variog(coords=coords, data=rssi,estimator.type='classical')
plot(var.notrend.geoR)