Выбор полосы пропускания для GWR
Я использую пакет spgwr в R для запуска модели географически взвешенной регрессии. Для модели я хочу автоматически вычислить пропускную способность, для чего я использую метод gwr.sel().
bandwidth <- gwr.sel(y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8,
data= data_sf,
coords=cbind(data_sf$x, data_sf$y),
longlat=T,
adapt=F)
При запуске кода всегда выдает предупреждение
...
na/inf replaced by maximum positive value. Bandwidth: 7859.3 CV score: NA
na/inf replaced by maximum positive value. Bandwidth: 7860.6 CV score: NA
na/inf replaced by maximum positive value. Bandwidth: 7861.1 CV score: NA
Bandwidth converged to upper bound: 7861.09
В описании пакета (http://rspatial.r-forge.r-project.org/spgwr/reference/gwr.cv.html) говорится, что: "Если регрессия содержит небольшой шаблон, пропускная способность будет сходиться к верхняя граница линейного поиска, которая является диагональю ограничивающего прямоугольника координат точки данных для "adap =FALSE" и 1 для "adap =TRUE";
Одна из переменных x имеет только 0 значений, что означает, что ни в одной из моих точек наблюдения нет. Когда я запускаю выбор полосы пропускания без этой переменной, я не получаю предупреждения.
Вопрос 1: должна ли функция gwr.sel() включать только те переменные, которые не содержат нулевых значений?
Проблема в том, что я хочу запустить регрессию автоматически для разных отметок времени, и я не могу быть уверен, что одна из переменных может иметь все нулевые значения. Если я исключу конкретную переменную, которая теперь имеет только нули, возможно, для следующей отметки времени это будет другая переменная только с нулевыми значениями. Но не могли бы вы сначала проверить все переменные, действительно ли они имеют значения?
Вопрос 2: И еще один вопрос, кстати, если мы говорим о выборе полосы пропускания. В каких единицах указан результат выбора полосы пропускания? Если у меня в качестве входных данных указаны координаты широты и долготы (например, -100, 40 для точки в США), будет ли ширина полосы в метрах, километрах, градусах?