How to plot Taylor diagram for four time series.

39 Ansichten (letzte 30 Tage)
Ojo Olusola
Ojo Olusola am 1 Jan. 2022
Kommentiert: wentong am 15 Nov. 2023
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 Kommentar
Turlough Hughes
Turlough Hughes am 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.

Melden Sie sich an, um zu kommentieren.

Antworten (1)

Meg Noah
Meg Noah am 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 Kommentare
Ahmad Bayhaqi
Ahmad Bayhaqi am 28 Mär. 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 am 15 Nov. 2023
RMSs(3) - sqrt(STDs(3).^2 + STDs(1)^2 - 2*STDs(3)*STDs(1).*CORs(3)) = 0 Why judge this

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu Detection, Range and Doppler Estimation finden Sie in Help Center und File Exchange

Produkte


Version

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by