sirmap.m

计量工具箱

function [out,Rbar,tstat]=sirmap(long,lat,xinp,yinp,nbcla,varargin)
% PURPOSE: This function links a map and one or two SIR scatterplots
%--------------------------------------------------------------
% USAGE: out = sirmap(long,lat,xinp,yinp,nbcla,direct,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
%          nbcla = number of classes
%          direct = optionnal 2x1 vector of integers or integer corresponding to the order of the edr directions selected for the plot. Default is [1;2]
%          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 either the map or the scatterplots 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 sirf.m, fastbinsmooth.m, setdens4.m, ols.m, selectmap.m, selectstat.m
%------------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit