今回使うパッケージ:geoR
https://cran.r-project.org/web/packages/geoR/geoR.pdf
In [ ]:
In [3]:
library(geoR)
In [4]:
data(s100)
In [5]:
class(s100)
Out[5]:
In [14]:
plot.geodata(s100)
In [15]:
points.geodata(s100)
In [7]:
cloud1 <- variog(s100, option = "cloud", max.dist = 1)
cloud2 <- variog(s100, option = "cloud", estimator.type = "modulus", max.dist = 1)
bin1 <- variog(s100, uvec = seq(0, 1, l = 11))
bin2 <- variog(s100, uvec = seq(0, 1, l = 11), estimator.type = "modulus")
In [9]:
plot(cloud1, main = "classical estimator")
In [11]:
plot(bin1, main = "classical estimator")
In [10]:
plot(cloud2, main = "modulus estimator")
In [16]:
plot(bin2, main = "modulus estimator")
In [19]:
bin1 <- variog(s100, uvec = seq(0, 1, l = 11), bin.cloud = TRUE)
bin2 <- variog(s100, uvec = seq(0, 1, l = 11), estimator.type = "modulus", bin.cloud = TRUE)
par(mfrow = c(1, 2))
plot(bin1, bin.cloud = TRUE, main = "classical estimator")
plot(bin2, bin.cloud = TRUE, main = "modulus estimator")
In [ ]: