У меня есть большой набор растровых данных, который я преобразовал в 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)