r中的执行时空模型

我对带有“spdep”包的SAR模型的执行时间有问题

我将同一维度的两个不同数据集传递给同一个函数，这需要非常不同的时间（几秒钟vs几个小时）

我给你写我的代码，如果你有任何想法请告诉我

谢谢 基亚拉

库（spdep）
数据（房屋、包装=“spData”）
hlw我得到：

库（spdep）
数据（房屋、包装=“spData”）
昏暗（房子）
hlw
library(spdep)
data(house, package="spData")



hlw<-nb2listw(LO_nb)
system.time( lagsarlm(log(price) ~ age, data=house, listw=hlw, type="lag", method="Matrix", trs=trMat))

#----------------------------------------------
library(spatstat)
d1<-100
d2<-100
n<-25357 

coord<- runifpoint(n,win=owin(c(0,d1),c(0,d2)))
mat<-cbind(coord$x,coord$y)

X<-rnorm(n,5,2)
Y<-rnorm(n,5,3)
d<-as.data.frame(cbind(X,Y))
cutoff<-  dnearneigh(mat,0,4)
t<-nb2listw(cutoff)
system.time( lagsarlm(Y ~ X, data=d, listw=t, type="lag", method="Matrix"))

library(spdep)
data(house, package = "spData")
dim(house)
hlw <- nb2listw(LO_nb)
# Number of nonzero links: 74874 
system.time(r1 <- lagsarlm(log(price) ~ age, data=house, listw=hlw,
                           type="lag", method="Matrix", trs=trMat))
# user  system elapsed 
# 0.50    0.05    0.55 

#----------------------------------------------
library(spatstat)
set.seed(21)
d1 <- 100
d2 <- 100
coord <- runifpoint(n, win = owin(c(0, d1), c(0, d2)))
mat <- cbind(coord$x, coord$y)
cutoff <- dnearneigh(mat, 0, 4)
t <- nb2listw(cutoff)
t
# Number of nonzero links: 3119764 
n <- 25357
X <- rnorm(n,5,2)
Y <- rnorm(n,5,3)
d <- as.data.frame(cbind(X,Y))

system.time(r2 <- lagsarlm(Y ~ X, listw = t, type = "lag", method = "Matrix"))
#    user  system elapsed 
# 156.47    8.54  168.25

# v2
d1 <- 200
d2 <- 200
coord <- runifpoint(n, win = owin(c(0, d1), c(0, d2)))
mat <- cbind(coord$x, coord$y)
cutoff <- dnearneigh(mat, 0, 4)
t <- nb2listw(cutoff)
t
# Number of nonzero links: 795054 

system.time(r2 <- lagsarlm(Y ~ X, listw = t, type = "lag", method = "Matrix"))
# user  system elapsed 
# 13.42    2.02   15.61