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

%% Calcolo delle statistiche descrittive per l'intero campione
stat={'mean' 'max' 'min' 'std'};
% il secondo argomento è vuoto poiché voglio le statistiche per tutte le
% unità
grpstats(X(:,{'PREZZO' 'POTENZA'}),[],stat)

%% 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)
varRaggrup={'RDS'};
disp(grpstats(X(:,[1:3 7]),varRaggrup,stat) )

% Osservazione: la variabile di raggruppamento può essere una combinazione
% di variabili. Ad esempio se se volevo avere le statistiche per ogni
% modalità di RDS e ogni modalità di Wireless  la sintassi sarebbe state
disp(grpstats(X(:,[1:3 5 7]),{'RDS' 'WIRELESS'},stat) )

% 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,stat) )


% 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],stat) )

%% Modo alternativo per trovare le statistiche per sottogruppi
% Inserisco come primo input l'intera table contenente tutte le variabili
% Con il name/pair 'DataVars',{'PREZZO' 'POTENZA'} vado a specificare che
% voglio le statistiche contenute dentro stat per le variabili prezzo e
% potenza
grpstats(X,varRaggrup,stat,'DataVars',{'PREZZO' 'POTENZA'})