Taylor Diagram

Plot a Taylor Diagram

Description

Plot a Taylor diagram from statistics values given by STDs (standard deviations), RMSs (centered root mean square difference) and CORs (correlation)

Ref: K. Taylor
Summarizing multiple aspects of model performance in a single diagram. JGR 2001

• MATLAB release: MATLAB 7.5 (R2007b)

[Pierre Hubertin Andrianirina]

 

[C Schwalm]

I had started to code this up myself, now I don't have to... These plots are wonderful for diagnostics and comparing numerous models. Well done.

[Georgina]

Does anyone have an example of the inputs required for this script. I am having trouble with it.
I get "Something''s wrong with the datas"
I think it is because I don't understand what is supposed to be in the 1st row of of each of the arrays.

[Damian Brady]

I agree with Georgina!

[Guillaume MAZE]

You'll find bellow links to .mat file with datas you can use to test the routine, and with the scripts I used myself to plot these datas (they are daily time series of sea surface net heat flux at a mooring in the North Atlantic, and results from numerical models to compare with)

hope this will be helpful
cheers
guillaume

.mat file is here: http://guillaumemaze.googlecode.com/svn/trunk/matlab/routines/taylordiag_egdata.mat
the example is here: http://guillaumemaze.googlecode.com/svn/trunk/matlab/routines/taylordiag_test.m

ps: I use the subroutine allstats in this peace of code to make things easier, you can get it here:
http://code.google.com/p/guillaumemaze/source/browse/#svn/trunk/matlab/routines

[Ben]

Thank you for the nice script and the examples.
I just have one question. What kind of function is the 'ptable' in the example: ax = ptable([2 3],[2 2;4 6]);
Thanks
Cheers Ben

[Jenny]

Hi Ben!

You can get the ptable script here:
http://code.google.com/p/guillaumemaze/wiki/matlab_routines_ptable

However, I still can't get the script to work. not with my own data and not with the example.

With my own data I get something like this:
??? Error using ==> taylordiag at 85
taylordiag.m : Something's wrong with the datas

With the example the error message is:
??? Index exceeds matrix dimensions.

Error in ==> taylordiag_test at 77
[pp tt axl] =
taylordiag(squeeze(statm(:,2)),squeeze(statm(:,3)),squeeze(statm(:,4)),...

Does anyone know what is wrong?
I think I've misunderstood something about what input to give the script...

Cheers!

[Mariac Valverde Brambila]

I agree with Jenny

[Roberto Venegas]

Hi, Can anybody let me know where I can get the subroutine allstats, I tried to get this from here: http://code.google.com/p/guillaumemaze/source/browse/#svn/trunk/matlab/routines

but doesnt work.
Thanks.
Robert.

[Carlos]

Hi, I noted that the script doesn't use the RMSd value for drawing the lines or placing the dots. What's wrong with it?

[Guillaume MAZE]

The whole idea of a Taylor diagram is the triangulation of a point given statistical properties it represents in a trigonometric context. A point needs only 2 parameters to be placed in the 2D diagram but we have 3 parameters which can be used. RMSd is therefore not used to plot the point but statistics/trigonometry requirements are checked and require RMSd.

[Georgina]

 

[PRASHANT SRIVASTAVA]

i am getting the same problem as jenny????
Any one know how to solve this?????

Prashant

[Guillaume MAZE]

To Jenny and others having the same problem,
If you get the error message about: "Something's wrong with the datas", this is because the statistics are wrong and trigonometric relations are wrong as well, see my comment from Jul.28th.
If you also have an error message running the test script, I can't figure out why it would be. The 'statm' matrix should be 5x4. I tried this test on many different platforms, removing all my custom path, etc. and it always works fine. The only reason why you would have an error and I can think about is that somewhere in your path a function overwrites something done in the script which makes it to behave not as expected.
So, please be sure you have the correct test datas and routines. To run the test you'll also need these two functions:
http://code.google.com/p/guillaumemaze/wiki/matlab_routines_ptable
http://code.google.com/p/guillaumemaze/wiki/matlab_routines_allstats

[Q.]

http://code.google.com/p/guillaumemaze/wiki/matlab_routines_ptable
http://code.google.com/p/guillaumemaze/wiki/matlab_routines_allstats

Page "matlab_routines_allstats" Not Found

So, you move to another places? why not put the test data here?

Thanks.

[Chiyuan]

Sorry, MAZE, I meet the same problem with jenny, and I found the generated 'Statm' matrix is 5x1, not 5x4....The matlab version I used is 7.0

[Guillaume MAZE]

New ptable and allstats location:
http://code.google.com/p/guillaumemaze/wiki/matlab_codes_graphicxPlots_ptable
and
http://code.google.com/p/guillaumemaze/wiki/matlab_codes_statistics_allstats

[Chiyuan]

It succeed now, thanks MAZE!

[Benrique]

aMAZEing. Thanks heaps.

[Shinobue111]

 

[Shinobue111]

I cant get the example to even work.. I have the latest allstats..

>> statm(:,2)
??? Index exceeds matrix dimensions.

[Ming Chen]

This is very handy! Thanks to Guillaume MAZE!

[Lang]

I could plot using example data sets. But there is no information of rmsd on plot other than some green lines. Whether the rmsd intervals label on the plot?

[Ming Chen]

Hi all, I can plot the diagram with only positive correlations. However, when I tried to plot both positive and negative correlations (Npan=2) with the same data, I got this error message:
??? Error using ==> set
Value must be numeric.

Error in ==> xlabel at 43
  set(h, 'FontAngle', get(ax, 'FontAngle'), ...

Error in ==> xlabel at 22
  h = xlabel(gca,varargin{:});

Error in ==> taylordiag at 399
ax(ix).handle =xlabel('Standard deviation','fontweight',axlabweight,...

Error in ==> All_3_years at 601
[hp ht axl] = taylordiag(STD_G21_08,RMS_G21_08,COR_G21_08,'Npan',2,'tickSTD',(0:0.25:1.5));

Is there anybody who has the same problem?

[Aisling]

This is a brilliant function. However, I am struggling to get the correlation lines eminating from the origin to plot when the correlation is negative (the points plot perfectly). Is there a way of doing this? Thanks.