adjacency_matricies_provincia_lat-long-weight.R

In AnthonyEbert/ItalyCovid19: What the Package Does (One Line, Title Case)

# Create adjacency matrix for countries

library(dplyr)

## Johns Hopkins data -------------

x = read.csv("https://raw.githubusercontent.com/pcm-dpc/COVID-19/master/dati-province/dpc-covid19-ita-province.csv")

## Get average lat and long for all countries ----------

x1 <- x %>%
  group_by(denominazione_provincia) %>%
  summarise(Long = mean(long), Lat = mean(lat))

# Create adjacency matrix

kol = ItalyCovid19::create_adjacency_matrix(x1$denominazione_provincia, x1$Lat, x1$Long)
diag(kol) = max(as.numeric(kol[is.finite(as.numeric(kol))])) + 1

write.csv(kol, "protezione-civile-download/adjacency_matrix_provincia_lat-long-weight.csv", row.names = TRUE)