从经度和模拟中模拟地理数据R的纬度

Question

我正在尝试对我正在测试的模型进行一些模拟，这些模型具有系统发育和地理网络。

我希望独立测试这些效果，以比较每个维度的方法估算能力。但是，我无法生成显示重要地理空间自动=相关性的连续变量。到目前为止，我已经使用了年平均温度（这是国家级数据），但这不允许我控制空间自相关的大小。另外它的缺点是不知道空间自相关的真实值。

我有一个经度和纬度点列表，我想知道是否有一种简单或常用的方法来生成一个连续响应变量，该变量将显示空间自相关，并能够控制该相关的大小？然而，优选在R内，只要它们可以转换成R-相容形式，其他方法也是可以接受的。

任何建议都将不胜感激，请询问是否有任何有用的信息。

Answer 1

遵循Dormann等人的方法。 （2007年），你可以这样做：

N <- 3000
p <- 1/N

# generate some points
set.seed(1234)
x.coord <- runif(N,0,100)
y.coord <- runif(N,0,100)
points <- cbind(x.coord,y.coord)

# distance matrix between points
Dd <- as.matrix(dist(points))

# weights matrix
w <- exp(-p * Dd)
Ww <- chol(w)

# errors
z <- t(Ww) %*% rnorm(N,0,1) 

# plot
df <- data.frame(x = x.coord, y = y.coord, z = z)
require(ggplot2)
ggplot(df, aes(x = x, y = y, col = z)) +
  geom_point() +
  scale_colour_gradient(low="red", high="white")

其中变量p控制自相关的大小（这里我将其设置为1/3000 = 0.000333）。 p = 0将没有相关性。

参考文献：Dormann，C.F.，McPherson，J.M.，Araujo，M.B.，Bivand，R.，Bolliger，J.，Carl，G.，...... Wilson，R。（2007）。在物种分布数据分析中考虑空间自相关的方法：综述。 Ecography，30（5），609-628。