Plotting Morans I values over time: the hard way

In this notebook the graph of the global Moran’s I and the differential Moran’s I over time is created.

1 Opening the libraries

knitr::opts_chunk$set(
  echo = TRUE,
  message = FALSE,
  warning = FALSE
)
library(ExPanDaR)
library(tidyverse)
library(readr)
library(reshape2)

# Change the presentation of decimal numbers to 4 and avoid scientific notation
options(prompt="R> ", digits=7, scipen=999)

2 Importing the data generated in Geoda

I created the moran.csv file manually by typing the output from Geoda, this was not too difficult given that there were just 17 different years. (that means I had to click less than 100 times!)

moran <- read_csv("moran.csv")

# Define function for  interactive table
create_dt <- function(x){
  DT::datatable(x,
                extensions = 'Buttons',
                options = list(dom = 'Blfrtip',
                               buttons = c('copy', 'csv', 'excel', 'pdf', 'print'),
                               lengthMenu = list(c(10,25,50,-1),
                                                 c(10,25,50,"All"))))
}

create_dt(moran)

3 Plotting the graph

plot1<-moran %>% 
  ggplot() +
   theme_minimal()+
    theme(text=element_text(family="Palatino"))+
     geom_line(aes(x = year, y = Morans_RGDPpc), alpha = 1) +
  geom_line(aes( x = year, y = differential_Morans), linetype= "dashed") +
    theme(text=element_text(family="Palatino"))+
  labs(subtitle = "",
     x="year", y="Moran's I")+
    geom_label(
    label="original", 
    x=2004,
    y=0.22,
    label.size = 0.35,
    color = "black",
  )+
   geom_label(
    label="differential", 
    x=2012,
    y=-0.03,
    label.size = 0.35,
    color = "black",
  )

plot(plot1)

4 appendix

other chunks of code to create similar graphs

##Subset the necessary columns

df<- moran
dd_sub = df[,c(1, 2,4)]

names(dd_sub)[2] <- "standard"
names(dd_sub)[3] <- "differential"

dd_sub

##Then rearrange your data frame


dd = melt(dd_sub, id=c("year"))

plot2<- ggplot(dd) + 
    theme_minimal()+
    theme(text=element_text(family="Palatino"))+
  geom_line(aes(x=year, y=value, colour=variable)) +
  scale_colour_manual(values=c("black","blue"))+
  labs(x="year", y="Moran's I")

plot(plot2)

##Subset the necessary columns
df<- df %>% 
  filter(year>1999)
dd_sub2 = df[,c(1, 2,4  )]
##Then rearrange your data frame
dd2 = melt(dd_sub2, id=c("year"))

plot3<- ggplot(dd2) +
  geom_line(aes(x=year, y=value, colour=variable)) +
  scale_colour_manual(values=c("black","purple"))+
  labs(x=" final year (t)", y="  Moran's I")

plot(plot3)