Проблема с проекцией координат для определения пространственной автокорреляции в R

Question

У нас есть набор данных, который содержит координаты широты и долготы, а также информацию об атрибутах, каждый в отдельном столбце, хранящийся как числовой. Эти координаты были геокодированы на основе географической системы координат WGS 1984. Мы знаем, что в наших данных имеется значительная пространственная автокорреляция, которую мы надеемся визуализировать на графике пузырьков с использованием пакета «sp». Мы моделируем наш пример другими пользователями, например, здесь: https://beckmw.wordpress.com/2013/01/07/breaking-the-rules-with-spatial-correlation/. Однако, когда мы пытаемся использовать координаты команды в "зр", мы продолжаем получать сообщение об ошибке:

Пример кода: Coords < - data.frame (свинец $ X, свинец $ Y) координаты (Coords) < - с ("lead6.X", "lead6.Y")

Ошибка в случае (NCHAR (projargs) == 0) projargs < - as.character (NA), где отсутствующее значение TRUE/FALSE, необходимый

Мы не можем загрузить наш прямой код, потому что он чувствителен и размещен в виртуальной среде без доступа в Интернет. У кого-нибудь есть идеи, почему это может произойти? Мы рассмотрели пакет proj4, но не можем понять, как указать систему проецирования (или это даже ошибка, которую мы получаем?). Если кто-либо знает о каких-либо других пакетах в R или способах визуализации пространственной автокорреляции, это тоже будет очень полезно.

Answer 1

Ваш код немного «странный»: кажется, вы пытаетесь построить набор данных, содержащий только координаты. AFAIU, вам может понадобиться что-то в этой строке:

data <- data.frame(lead$X, lead$Y, lead$Z) 

с свинцом $ Z, соответствующим общим «переменный» Вы хотите, чтобы проверить, то

coordinates(data) <- c('X','Y')` 
proj4string(data) <- "+init=epsg:4326" 

, который должен дать вам правильный " SpatialPointsDataframe "с полными географическими координатами WGS84 (первая строка также может быть отброшена, и вы сохраните все переменные в данных пространственной точки данных).

НТН