% PURPOSE: An example of using semo()
%          spatial error model
%          on a small data set                   
%---------------------------------------------------
% USAGE: semo_d, (see semo_d2 for a large data set)
%---------------------------------------------------

% load Anselin (1988) Columbus neighborhood crime data
load anselin.dat;

y = anselin(:,1);
n = length(y);
x = [ones(n,1) anselin(:,2:3)];

% load Anselin (1988) 1st order contiguity matrix
load wmat.dat; % standardized
W = wmat;

vnames = strvcat('crime','constant','income','hvalue');
          
% do sem regression
res = semo(y,x,W);
prt(res,vnames);