2

我有一个观察点列表，我希望在这些点上估计地理加权回归。然后，我有一组单独且不同的点，我想使用这个估计的模型在这些点上进行预测。这是使用以下代码模拟的：

library(spgwr)
library(sp)

set.seed(1)
e <- runif(100,0,1)
n <- runif(100,0,1)
y <- runif(100,0,1)
x1 <- runif(100,0,1)
x2 <- runif(100,0,1)
est.df <- data.frame(cbind(y,x1,x2,e,n))
coordinates(est.df) <- ~ e + n

e < -runif(50,0,1)
n <- runif(50,0,1)
x1 <- runif(50,0,1)
x2 <- runif(50,0,1)
pred.df <- data.frame(cbind(x1,x2,e,n))
coordinates(pred.df) <- ~ e + n

remove(list = c("e","n","y","x1","x2"))

S.gwr <- gwr(y~x1+x2,
           data=est.df,
           coords=cbind(est.df$e, est.df$n),
           bandwidth=0.1,
           predictions = TRUE,
           fit.points=pred.df)

dim(S.gwr$SDF@data)
length(S.gwr$lm$fitted.values)

我的问题是 S.gwr 结构中唯一的拟合值是长度 100，并且是用于估计模型的 100 个点，而不是 pred.df 中我想要预测 y_hat 的 50 个点。谁能帮我指出没有为我提供 pred.df 中点的 50 个预测的问题？谢谢。

4

0 回答 0