pcamap.m

来自「计量工具箱」· M 代码 · 共 34 行

function [out,inertia,casecoord,varcoord]=pcamap(long,lat,data,varargin)
% PURPOSE: This function links a map and a pca plot
%------------------------------------------------------------------------
% USAGE: [out,inertia,casecoord,varcoord]=pcamap(long,lat,data,direct,weight,metric,center,reduc,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
%           data = n x 1 vector of the variable to study on the first axis
%           direct = 1 x 2 vector of the principal component number to be plotted on the 2 axis
%           weight : optional (n x 1) weight vector. The default is weight=(1/n,...,1/n)'.
%           metric : optional (p x p) matrix. The default is the identity matrix.
%           center : optional boolean: if center=1 : the pca is centered (default).
%                                      if center=0 : the pca is not centered.
%           reduc : optional boolean: if reduc=1 : the pca is standardized (default).
%                                     if reduc=0 : the pca is not standardized.
%           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 pca plot 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 genpca.m
% -----------------------------------------------------------------------
% Christine Thomas-Agnan, Anne Ruiz-Gazen, Julien Moutel
% June 2003
% Universit