How to plot Taylor diagram for four time series.

69 views (last 30 days)
A =readtable
Warning: Column headers from the file were modified to make them valid MATLAB identifiers before creating variable names for the table. The original column headers are saved in the VariableDescriptions property.
Set 'VariableNamingRule' to 'preserve' to use the original column headers as table variable names.
A = 166×6 table
Time insitu era5 ncar_ncep merra2 ceres ___________ ______ ______ _________ ______ _____ 01-Mar-2000 101 242.59 145.1 123.7 135.5 01-Apr-2000 117.44 211.69 162.1 113.5 130.5 01-May-2000 126.48 187.73 155.6 121.5 133.9 01-Jun-2000 120.7 133.91 143.1 131.7 123.5 01-Jul-2000 106.31 123.15 132.8 134 103.3 01-Aug-2000 108.62 133.91 135.5 135.7 104.8 01-Sep-2000 114.76 147.34 136.8 139.1 121.9 01-Oct-2000 125.73 182.99 139.9 131.2 140.1 01-Nov-2000 125.48 223.84 141.3 97.5 123.7 01-Dec-2000 111.72 210.65 131.9 89.5 108.5 01-Jan-2001 101.74 222.8 132 92.9 104.2 01-Feb-2001 100.56 233.1 137.1 104 122.3 01-Mar-2001 109 235.42 159.7 107.4 144.6 01-Apr-2001 125.2 196.64 152.1 113.2 130.8 01-May-2001 130.99 168.17 157.5 114.9 137.6 01-Jun-2001 121.22 142.71 143 111.1 114.9
('Sample data.xlsx')
I have the above five different datasets, i want matlab code to execute Taylor diagram for the dataset taking correlation between insitu and each of the other four model data in the sample data. Kindly help. Thanks.
  1 Comment
Turlough Hughes
Turlough Hughes on 2 Jan 2022
There doesn't appear to be any function for making a taylor diagram in MATLAB. I suggest trying on of these FEX submissions.

Sign in to comment.

Answers (1)

Meg Noah
Meg Noah on 2 Jan 2022
Here's how:
Download this:
Then use this code:
myData = readtable('Sample data.xlsx');
% STATM(1,:) => Mean
% STATM(2,:) => Standard Deviation (scaled by N)
% STATM(3,:) => Centered Root Mean Square Difference (scaled by N)
% STATM(4,:) => Correlation
stat1 = allstats(myData.insitu,myData.era5);
stat2 = allstats(myData.insitu,myData.ncar_ncep);
stat3 = allstats(myData.insitu,myData.merra2);
stat4 = allstats(myData.insitu,myData.ceres);
STDs = horzcat(stat1(2,:),stat2(2,2),stat3(2,2),stat4(2,2));
RMSs = horzcat(stat1(3,:),stat2(3,2),stat3(3,2),stat4(3,2));
CORs = horzcat(stat1(4,:),stat2(4,2),stat3(4,2),stat4(4,2));
[pp tt axl] = taylordiag(STDs,RMSs,CORs);
for ii = 1 : length(tt)
set(tt(ii),'fontsize',9,'fontweight','bold')
set(pp(ii),'markersize',12)
set(tt(ii),'String',strrep(myData.Properties.VariableNames{ii+1},'_',' '));
end
  5 Comments
Ahmad Bayhaqi
Ahmad Bayhaqi on 28 Mar 2023
Hi @Meg Noah , how about the ticklabels of STD, Cor Coef and RMSD? do you know how to change their fontsize and fontweight?
Thank you
wentong
wentong on 15 Nov 2023
RMSs(3) - sqrt(STDs(3).^2 + STDs(1)^2 - 2*STDs(3)*STDs(1).*CORs(3)) = 0 Why judge this

Sign in to comment.

Categories

Find more on Particle & Nuclear Physics in Help Center and File Exchange

Products


Release

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!