How to plot Taylor diagram for four time series.

68 visualizaciones (últimos 30 días)
Ojo Olusola
Ojo Olusola el 1 de En. de 2022
Comentada: wentong el 15 de Nov. de 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 comentario
Turlough Hughes
Turlough Hughes el 2 de En. de 2022
There doesn't appear to be any function for making a taylor diagram in MATLAB. I suggest trying on of these FEX submissions.

Iniciar sesión para comentar.

Respuestas (1)

Meg Noah
Meg Noah el 2 de En. de 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 comentarios
Ahmad Bayhaqi
Ahmad Bayhaqi el 28 de Mzo. de 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 el 15 de Nov. de 2023
RMSs(3) - sqrt(STDs(3).^2 + STDs(1)^2 - 2*STDs(3)*STDs(1).*CORs(3)) = 0 Why judge this

Iniciar sesión para comentar.

Categorías

Más información sobre Particle & Nuclear Physics en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by