Выполнение географической взвешенной регрессии на подмножестве данных из большого растра с использованием пакета spgwr
У меня есть большой набор растровых данных, который я преобразовал в SpatialPointsDataFrame и оттуда в обычный фрейм данных. Мне нужно запустить географическую взвешенную регрессию (GWR) с помощью пакета spgwr. Я успешно вычислил пропускную способность GWR с помощью функции gwr.sel(), но когда я пытаюсь запустить GWR с помощью gwr(), я получаю сообщение об ошибке, в котором говорится, что R не может выделить вектор размером 25 Гб. Я думал, что, возможно, можно запустить локальную регрессию GWR только на случайной подвыборке всего фрейма данных, сохраняя все остальные точки в качестве значений окрестности для регрессии.
Вот что у меня есть, но я не уверен, как запустить GWR на подвыборке, не избавляясь от всех остальных вопросов.
##Calculate kernel bandwidth
gwr.bw <- gwr.sel(merge.df$nsEVI_YAproj ~ merge.df$aDEM_YAproj, data = merge.df, coords = cbind(merge.df$x, merge.df$y), adapt = TRUE)
##Run GWR model
gwr <- gwr(merge$nsEVI_YAproj ~ merge$aDEM_YAproj, data = merge.df, coords = cbind(merge.df$x, merge.df$y), adapt = gwr.bw, hatmatrix = TRUE, se.fit = TRUE)