Тест - прогноз разделения поездов на пространственную модель в R - PullRequest
Новая
       42

Тест - прогноз разделения поездов на пространственную модель в R

0 голосов
/ 31 марта 2020

Итак, мне нужно создать прогноз расщепления тест-поезда для модели пространственной ошибки.

Мы оцениваем модель в 80% наблюдений, а затем позволяем модели прогнозировать 20%, однако я не могу получить прогноз для остальные 20%

, поэтому я бы не стал просить совета

Пример: 



library(spdep)

data(oldcol)

index_train <- sample(1:nrow(COL.OLD), round(0.80 * nrow(COL.OLD)))
df_train <- COL.OLD[index_train, ]

df_test <- COL.OLD[-index_train, ]


# Creating W matrix for train subset
CORD = cbind(COL.OLD$X, COL.OLD$Y)
cns <- knearneigh(CORD[index_train, ], k=4, longlat=T) 
scnsn <- knn2nb(cns, row.names = NULL, sym = T) 
W <- nb2listw(scnsn)


model <- lagsarlm(CRIME ~ INC + HOVAL, data=df_train, W)
summary(model)


# Creating network matrix for test subset
cns <- knearneigh(CORD[-index_train, ], k=4, longlat=T) 
scnsn <- knn2nb(cns, row.names = NULL, sym = T) 
W <- nb2listw(scnsn)

predict(model, newdata = df_test, listw = W)

возвращает ошибку, однако прогноз по линейной регрессии lm довольно прямой вперед

...