semmap.m

计量工具箱

function [out]=semmap(long,lat,xinp,yinp,W,varargin)
% PURPOSE: adjust a spatial error model on a subregion selected on the map
%--------------------------------------------------------------
% USAGE: out = semmap(long,lat,xinp,yinp,W,carte,label,symbol)
%   where: 
%          long = n x 1 vector of coordinates on the first axis
%          lat = n x 1 vector of coordinates on the second axis
%          xinp = explanatory variable (matrix n x p)  
%          yinp = n x 1 vector of the dependent variable
%          W = contiguity matrix
%          vnames  = an optional vector of variable names,
%                    e.g. vnames = strvcat('y','const','x1','x2');
%          carte = n x 2 optionnal matrix giving the polygons of the edges of the spatial units
%          label = n x 1 optionnal variable used to label selected observations
%           symbol = * symbol=1 : selected spatial units are marked with a different symbol
%                    * symbol=0 : selected spatial units are marked with a different color only (default)
%--------------------------------------------------------------
% OUTPUTS: out = (n x 1) 0-1 variable: selected spatial units are marked with a 1  
%--------------------------------------------------------------
% MANUAL: Select points on the map by clicking with the left mouse button
%         You can select points inside a polygon: - right click to set the first point of the polygon
%                                                 - left click to set the other points
%                                                 - right click to close the polygon
%         Selection is lost when you switch graph
%         To quit, click the middle button or press 'q'
% -----------------------------------------------------------------------
% uses sem.m, selectmap.m
%-----------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit