Блог

Главная - Вопросы по программированию - Создание матрицы смежности на уровне округа для разных состояний

Создание матрицы смежности на уровне округа для разных состояний

#r #matrix

#r #матрица

Вопрос:

Я пытаюсь создать матрицу смежности для округов в AL, SC и GA, и я попробовал один подход, который, казалось, хорошо работал. Я использовал следующие библиотеки и код для создания матрицы смежности (Adj_mat).

 library(tidyverse)
library(spdep)
library(urbnmapr)

counties_sf <- get_urbn_map("counties", sf = TRUE)

counties_south <- filter(counties_sf, state_abbv %in% c("AL", "GA", "SC" ))

south_counties_polylist <- poly2nb(counties_south)

Adj_mat <- nb2mat(south_counties_polylist, style = "B", zero.policy = T)           # Adjacency matrix

Общие смежности для этой матрицы 1526 получены следующим образом
m1 <- apply(Adj_mat, 2, sum) .

Это было немного тревожно, поскольку у меня есть другая матрица смежности для AL, GA и SC, которую я использую уже некоторое время и имеет 1528 полные смежности. У меня нет кода, который использовался для создания этой матрицы, и я немного не уверен, был ли мой подход неправильным или существующая матрица смежности, которую я использовал, неверна.

На основе описания пакета urbnmapr библиотека использует шейп-файлы из Бюро переписи населения США. Я хотел бы иметь возможность самостоятельно создавать матрицы для разных состояний и был бы признателен за любые указания, чтобы убедиться, что я делаю это правильно. Спасибо!

Комментарии:

1. Матрицы смежности очень чувствительны к условиям, например, когда округа разделяют только одну точку, или когда четыре округа встречаются в одной точке, или когда оцифрованные границы отличаются на крошечную величину (например, сантиметры на земле). Можете ли вы найти 2 дополнительных смежности в наборе 1528 и отобразить округа? (Это может быть только одна смежность и ее обратная). Это может объяснить это.

2. @Spacedman Спасибо за отзывы и предложения! Я считаю, что это одна смежность и ее обратная из-за симметрии. Я посмотрю, смогу ли я найти, какая это смежность.

3.Я попробовал полигоны округа из counties() функции в tigris r пакете, и это сработало хорошо. Мои две матрицы смежности для AL, GA и SC теперь согласуются. Tigris использует шейп-файлы TIGER / Lines из переписи населения США.