% PURPOSE: An example of using arc_moranplot() 
%          with an ArcView shape file containing 
%          30 polygons for China provinces 
%---------------------------------------------------
% USAGE: arc_moranplotd
%---------------------------------------------------

clear all;

filename = '..\shape_files\eu138';

results = shape_read(filename);

nobs = results.nobs;
% pull out some variables to map
ecu95 = results.data(:,20);
ecu80 = results.data(:,5);
igrwth = ecu95./ecu80;
igrwth = igrwth - ones(nobs,1);        
igrwth = igrwth/16; % annual percentage growth in per capita income

vname = 'per capita income growth';

variable = igrwth - mean(igrwth);

latt = results.xc;
long = results.yc;

[j,W,j] = xy2cont(long,latt);

options.vname = vname;
%options.labels = 1;
options.mapmenu = 1;
options.legendmenu = 1;

arc_moranplot(variable,W,results,options);