How to plot Taylor diagram for four time series.

52 次查看(过去 30 天)
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 个评论
Turlough Hughes
Turlough Hughes 2022-1-2
There doesn't appear to be any function for making a taylor diagram in MATLAB. I suggest trying on of these FEX submissions.

请先登录,再进行评论。

回答(1 个)

Meg Noah
Meg Noah 2022-1-2
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 个评论
Ahmad Bayhaqi
Ahmad Bayhaqi 2023-3-28
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 2023-11-15
RMSs(3) - sqrt(STDs(3).^2 + STDs(1)^2 - 2*STDs(3)*STDs(1).*CORs(3)) = 0 Why judge this

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by