There is a script by Guillaume Maze available http://code.google.com/p/guillaumemaze/

In chase you can't open it either here is the copy of it:

% TAYLORDIAG Plot a Taylor Diagram

%

% [hp ht axl] = taylordiag(STDs,RMSs,CORs,['option',value])

%

% Plot a Taylor diagram from statistics of different series.

%

% INPUTS:

% STDs: Standard deviations

% RMSs: Centered Root Mean Square Difference

% CORs: Correlation

%

% Each of these inputs are one dimensional with same length. First

% indice corresponds to the reference serie for the diagram. For exemple

% STDs(1) is the standard deviation of the reference serie and STDs(2:N)

% are the standard deviations of the other series.

%

% Note that by definition the following relation must be true for all series i:

% RMSs(i) - sqrt(STDs(i).^2 + STDs(1)^2 - 2*STDs(i)*STDs(1).*CORs(i)) = 0

% This relation is checked and if not verified an error message is sent. Please see

% Taylor's JGR article for more informations about this.

% You can use the ALLSTATS function to avoid this to happen, I guess ;-). You can get

% it somewhere from: http://codes.guillaumemaze.org/matlab

%

% OUTPUTS:

% hp: returns handles of plotted points

% ht: returns handles of the text legend of points

% axl: returns a structure of handles of axis labels

%

% LIST OF OPTIONS:

% For an exhaustive list of options to customize your diagram, please call the function

% without arguments:

% >> taylordiag

%

% SHORT TUTORIAL (see taylordiag_test.m for more informations):

% An easy way to get compute inputs is to use the ALLSTATS function you can get from:

% http://codes.guillaumemaze.org/matlab

% Let's say you gathered all the series you want to put in the Taylor diagram in a

% single matrix BUOY(N,nt) with N the number of series and nt their (similar) length.

% If BUOY(1,:) is the serie of reference for the diagram:

% for iserie = 2 : size(BUOY,1)

% S = allstats(BUOY(1,:),BUOY(iserie,:));

% MYSTATS(iserie,:) = S(:,2); % We get stats versus reference

% end%for iserie

% MYSTATS(1,:) = S(:,1); % We assign reference stats to the first row

% Note that the ALLSTATS function can handle NaNs, so be careful to compute statistics

% with enough points !

% Then you're ready to simply run:

% taylordiag(MYSTATS(:,2),MYSTATS(:,3),MYSTATS(:,4));

%

% REF: K. Taylor

% Summarizing multiple aspects of model performance in a single diagram

% Journal of Geophysical Research-Atmospheres, 2001, V106, D7.

%

% Rev. by Guillaume Maze on 2010-02-10: Help more helpful ! Options now displayed by call.

% Copyright (c) 2008 Guillaume Maze.

% http://codes.guillaumemaze.org

% All rights reserved.

%

% Redistribution and use in source and binary forms, with or without

% modification, are permitted provided that the following conditions are met:

% * Redistributions of source code must retain the above copyright notice, this list of

% conditions and the following disclaimer.

% * Redistributions in binary form must reproduce the above copyright notice, this list

% of conditions and the following disclaimer in the documentation and/or other materials

% provided with the distribution.

% * Neither the name of the Laboratoire de Physique des Oceans nor the names of its contributors may be used

% to endorse or promote products derived from this software without specific prior

% written permission.

%

% THIS SOFTWARE IS PROVIDED BY Guillaume Maze ''AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES,

% INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A

% PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL Guillaume Maze BE LIABLE FOR ANY

% DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT

% LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR

% BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,

% STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE

% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

%

function varargout = taylordiag(varargin)

%%

if nargin == 0

disp_optionslist;

return

else

narg = nargin - 3;

if mod(narg,2) ~=0

error('taylordiag.m : Wrong number of arguments')

end

end

STDs = varargin{1};

RMSs = varargin{2};

CORs = varargin{3};

%% CHECK THE INPUT FIELDS:

apro = 100;

di = fix(RMSs*apro)/apro - fix(sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs)*apro)/apro;

if find(di~=0)

% help taylordiag.m

ii = find(di~=0);

if length(ii) == length(di)

error(sprintf('taylordiag.m : Something''s wrong with ALL the datas\nYou must have:\nRMSs - sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs) = 0 !'))

else

error(sprintf('taylordiag.m : Something''s wrong with data indice(s): [%i]\nYou must have:\nRMSs - sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs) = 0 !',ii))

end

end

%% IN POLAR COORDINATES:

rho = STDs;

theta = real(acos(CORs));

dx = rho(1); % Observed STD

%%

%% BEGIN THE PLOT HERE TO GET AXIS VALUES:

hold off

cax = gca;

tc = get(cax,'xcolor');

%ls = get(cax,'gridlinestyle');

ls = '-'; % DEFINE HERE THE GRID STYLE

next = lower(get(cax,'NextPlot'));

%% LOAD CUSTOM OPTION OF AXE LIMIT:

nopt = narg/2; foundrmax = 0;

for iopt = 4 : 2 : narg+3

optvalue = varargin{iopt+1};

switch lower(varargin{iopt}), case 'limstd', rmax = optvalue; foundrmax=1; end

end

% make a radial grid

hold(cax,'on');

if foundrmax==0

maxrho = max(abs(rho(:)));

else

maxrho = rmax;

end

hhh = line([-maxrho -maxrho maxrho maxrho],[-maxrho maxrho maxrho -maxrho],'parent',cax);

set(cax,'dataaspectratio',[1 1 1],'plotboxaspectratiomode','auto')

v = [get(cax,'xlim') get(cax,'ylim')];

ticks = sum(get(cax,'ytick')>=0);

delete(hhh);

% check radial limits and ticks

rmin = 0;

if foundrmax == 0;

rmax = v(4);

end

rticks = max(ticks-1,2);

if rticks > 5 % see if we can reduce the number

if rem(rticks,2) == 0

rticks = rticks/2;

elseif rem(rticks,3) == 0

rticks = rticks/3;

end

end

rinc = (rmax-rmin)/rticks;

tick = (rmin+rinc):rinc:rmax;

%% LOAD DEFAULT PARAMETERS:

if find(CORs=0 & th= i ,1);

ig = find(i*cos(th)+dx

I'm also facing similar problem

Hi Asid,

You can download from here!!

http://guillaumemaze.googlecode.com/svn/trunk/matlab/codes/statistics/taylordiag.m

Hi Jafar

Thank you for your reply. However, I am unable to reach mentioned web page.

There is a script by Guillaume Maze available http://code.google.com/p/guillaumemaze/

In chase you can't open it either here is the copy of it:

% TAYLORDIAG Plot a Taylor Diagram

%

% [hp ht axl] = taylordiag(STDs,RMSs,CORs,['option',value])

%

% Plot a Taylor diagram from statistics of different series.

%

% INPUTS:

% STDs: Standard deviations

% RMSs: Centered Root Mean Square Difference

% CORs: Correlation

%

% Each of these inputs are one dimensional with same length. First

% indice corresponds to the reference serie for the diagram. For exemple

% STDs(1) is the standard deviation of the reference serie and STDs(2:N)

% are the standard deviations of the other series.

%

% Note that by definition the following relation must be true for all series i:

% RMSs(i) - sqrt(STDs(i).^2 + STDs(1)^2 - 2*STDs(i)*STDs(1).*CORs(i)) = 0

% This relation is checked and if not verified an error message is sent. Please see

% Taylor's JGR article for more informations about this.

% You can use the ALLSTATS function to avoid this to happen, I guess ;-). You can get

% it somewhere from: http://codes.guillaumemaze.org/matlab

%

% OUTPUTS:

% hp: returns handles of plotted points

% ht: returns handles of the text legend of points

% axl: returns a structure of handles of axis labels

%

% LIST OF OPTIONS:

% For an exhaustive list of options to customize your diagram, please call the function

% without arguments:

% >> taylordiag

%

% SHORT TUTORIAL (see taylordiag_test.m for more informations):

% An easy way to get compute inputs is to use the ALLSTATS function you can get from:

% http://codes.guillaumemaze.org/matlab

% Let's say you gathered all the series you want to put in the Taylor diagram in a

% single matrix BUOY(N,nt) with N the number of series and nt their (similar) length.

% If BUOY(1,:) is the serie of reference for the diagram:

% for iserie = 2 : size(BUOY,1)

% S = allstats(BUOY(1,:),BUOY(iserie,:));

% MYSTATS(iserie,:) = S(:,2); % We get stats versus reference

% end%for iserie

% MYSTATS(1,:) = S(:,1); % We assign reference stats to the first row

% Note that the ALLSTATS function can handle NaNs, so be careful to compute statistics

% with enough points !

% Then you're ready to simply run:

% taylordiag(MYSTATS(:,2),MYSTATS(:,3),MYSTATS(:,4));

%

% REF: K. Taylor

% Summarizing multiple aspects of model performance in a single diagram

% Journal of Geophysical Research-Atmospheres, 2001, V106, D7.

%

% Rev. by Guillaume Maze on 2010-02-10: Help more helpful ! Options now displayed by call.

% Copyright (c) 2008 Guillaume Maze.

% http://codes.guillaumemaze.org

% All rights reserved.

%

% Redistribution and use in source and binary forms, with or without

% modification, are permitted provided that the following conditions are met:

% * Redistributions of source code must retain the above copyright notice, this list of

% conditions and the following disclaimer.

% * Redistributions in binary form must reproduce the above copyright notice, this list

% of conditions and the following disclaimer in the documentation and/or other materials

% provided with the distribution.

% * Neither the name of the Laboratoire de Physique des Oceans nor the names of its contributors may be used

% to endorse or promote products derived from this software without specific prior

% written permission.

%

% THIS SOFTWARE IS PROVIDED BY Guillaume Maze ''AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES,

% INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A

% PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL Guillaume Maze BE LIABLE FOR ANY

% DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT

% LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR

% BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,

% STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE

% OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

%

function varargout = taylordiag(varargin)

%%

if nargin == 0

disp_optionslist;

return

else

narg = nargin - 3;

if mod(narg,2) ~=0

error('taylordiag.m : Wrong number of arguments')

end

end

STDs = varargin{1};

RMSs = varargin{2};

CORs = varargin{3};

%% CHECK THE INPUT FIELDS:

apro = 100;

di = fix(RMSs*apro)/apro - fix(sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs)*apro)/apro;

if find(di~=0)

% help taylordiag.m

ii = find(di~=0);

if length(ii) == length(di)

error(sprintf('taylordiag.m : Something''s wrong with ALL the datas\nYou must have:\nRMSs - sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs) = 0 !'))

else

error(sprintf('taylordiag.m : Something''s wrong with data indice(s): [%i]\nYou must have:\nRMSs - sqrt(STDs.^2 + STDs(1)^2 - 2*STDs*STDs(1).*CORs) = 0 !',ii))

end

end

%% IN POLAR COORDINATES:

rho = STDs;

theta = real(acos(CORs));

dx = rho(1); % Observed STD

%%

%% BEGIN THE PLOT HERE TO GET AXIS VALUES:

hold off

cax = gca;

tc = get(cax,'xcolor');

%ls = get(cax,'gridlinestyle');

ls = '-'; % DEFINE HERE THE GRID STYLE

next = lower(get(cax,'NextPlot'));

%% LOAD CUSTOM OPTION OF AXE LIMIT:

nopt = narg/2; foundrmax = 0;

for iopt = 4 : 2 : narg+3

optvalue = varargin{iopt+1};

switch lower(varargin{iopt}), case 'limstd', rmax = optvalue; foundrmax=1; end

end

% make a radial grid

hold(cax,'on');

if foundrmax==0

maxrho = max(abs(rho(:)));

else

maxrho = rmax;

end

hhh = line([-maxrho -maxrho maxrho maxrho],[-maxrho maxrho maxrho -maxrho],'parent',cax);

set(cax,'dataaspectratio',[1 1 1],'plotboxaspectratiomode','auto')

v = [get(cax,'xlim') get(cax,'ylim')];

ticks = sum(get(cax,'ytick')>=0);

delete(hhh);

% check radial limits and ticks

rmin = 0;

if foundrmax == 0;

rmax = v(4);

end

rticks = max(ticks-1,2);

if rticks > 5 % see if we can reduce the number

if rem(rticks,2) == 0

rticks = rticks/2;

elseif rem(rticks,3) == 0

rticks = rticks/3;

end

end

rinc = (rmax-rmin)/rticks;

tick = (rmin+rinc):rinc:rmax;

%% LOAD DEFAULT PARAMETERS:

if find(CORs=0 & th= i ,1);

ig = find(i*cos(th)+dx