SSFA residuals Moran plot

Share:

Description

This function allows to plot the residuals of the object against their spatially lagged values, augmented by reporting the summary of influence measures for the linear relationship between the data and the lag.

Usage

1
plot_moran(x, main, xlab, ylab, labels, listw, ...)

Arguments

x

an object of class ssfa.

main

an overall title for the plot.

xlab

a label for the x axis.

ylab

a label for the y axis.

labels

character labels for points with high influence measures, if set to FALSE, no labels are plotted for points with large influence.

listw

a listw object from nb2listw (see nb2listw).

...

arguments to be passed to methods, such as graphical parameters (see par).

References

Anselin, L. (1995). Local indicators of spatial association, Geographical Analysis, 27, 93-115.

Anselin, L. (1996). The Moran scatterplot as an ESDA tool to assess local instability in spatial association. pp. 111-125 in M. M. Fischer, H. J. Scholten and D. Unwin (eds) Spatial analytical perspectives on GIS, London, Taylor and Francis.

See Also

moran.plot

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
library(ssfa) 
data(SSFA_example_data)
data(Italian_W)

#### SFA and SSFA comparison ###
sfa <- ssfa(log_y ~ log_x, data = SSFA_example_data, data_w=Italian_W, 
             form = "production", par_rho=FALSE)
ssfa <- ssfa(log_y ~ log_x, data = SSFA_example_data, data_w=Italian_W, 
             form = "production", par_rho=TRUE)

moran.test(residuals.ssfa(sfa), sfa$list_w)
moran.test(residuals.ssfa(ssfa), ssfa$list_w)

plot_moran(sfa, listw=sfa$list_w)
plot_moran(ssfa, listw=ssfa$list_w)