# Moran scatterplot of the regional unemployment rate (u5EX)
u5EX.lm = lm(Lu5EX~u5EX)
u5EX.lm
Lu5EX.fit = u5EX.lm$fitted.values
mu5EX = mean(u5EX)
su5EX = sqrt((n-1)/n)*sd(u5EX)
zu5EX = (u5EX-mu5EX)/su5EX
Lzu5EX = lag.listw(W5EX.lw, zu5EX)
Lzu5EX.fit = (Lu5EX.fit-mu5EX)/su5EX
plot(zu5EX, Lzu5EX, main="Moran scatterplot (5 regions)",xlim=c(-1.5, 1.5), ylim=c(-1, 1), cex.main=1.15, font.main=1)
abline(h=0, v=0)
lines(zu5EX, Lzu5EX.fit)
box(which="figure")