%% Caricamento dati
X=readtable('htgruppi.xlsx','Sheet','dati','Range','A1:H26','ReadRowNames',true,'Format','auto');


%% Calcolo delle statistiche richieste in base alla presenza/assenza di RDS
% Calcolo delle statistiche descrittive delle variabili [1:3] utilizzando
% come variabile di raggrupamento la variabile 7 ('RDS')
% Osservazione: la funzione grpstats deve avere come input una table che
% contiene le variabili su cui calcolare le statistiche descrittive e la
% variabile di raggruppamento (in questo caso RDS)
disp(grpstats(X(:,[1:3 7]),{'RDS'},{'mean' 'max' 'min' 'std'}) )

% Modo alternativo per identificare la variabile di raggrupamento.
% Il numero 4 poiché nella table X(:,[1:3 7]) la variabile RDS è in
% posizione 4
disp(grpstats(X(:,[1:3 7]),4,{'mean' 'max' 'min' 'std'}) )


% Modo alternativo per identificare la variabile di raggrupamento.
% [false false false true]
% true in posizione 4 poiché nella table X(:,[1:3 7]) la variabile RDS è in
% posizione 4
disp(grpstats(X(:,[1:3 7]),[false false false true],{'mean' 'max' 'min' 'std'}) )