%% IMPORTAZIONE DATI
X=readtable('ex23_DE.xls','Sheet','Foglio1','Range','D1:F108');
retr=X{:,3};
sesso=X{:,1};

%% Creazione della distribuzione di frequenze e delle etichette per ogni classe
retrclassi=cell(length(retr),1);

eti=({'<=1000]'; '(1000-1200]'; '(1200 1500]'; '(1500-2200]'; '>2200'});
freq=zeros(5,1);

boo=retr<=1000;
retrclassi(boo)=eti(1);
freq(1)=sum(boo);

boo=retr>=1000 & retr<=1200;
retrclassi(boo)=eti(2);
freq(2)=sum(boo);

boo=retr>=1200 & retr<=1500;
retrclassi(boo)=eti(3);
freq(3)=sum(boo);

boo=retr>1500 & retr<=2200;
retrclassi(boo)=eti(4);
freq(4)=sum(boo);

boo=retr>2200;
retrclassi(boo)=eti(5);
freq(5)=sum(boo);

Df=array2table(freq,'RowNames',eti,'VariableNames',{'Distr frequenze'});

%% Creazione della tabella a doppia entrata

% questa istruzione serve per stabilire che l'ordine in cui devono essere
% mostrate le classi nelle successive elaborazioni è quello riportato nella
% cell eti
retrclassic=categorical(retrclassi,eti);
[tbl,chi2,p,labels] =crosstab(retrclassic,sesso);
[r,c]=size(tbl);
etirighe=labels(1:r,1);
eticolonne=labels(1:c,2);
TB=array2table(tbl,'RowNames',etirighe,'VariableNames',eticolonne);

% Xall table che contiene 3 colonne
% Sesso, Retribuzione e classi di retribuzione
Xall=[X(:,[1 3 ]) table(retrclassic)];
% calcolo della media della retribuzione per ciascuna combinazione di Sesso
% e classe di retribuzione
grpstats(Xall,{'Sesso' 'retrclassic'},{'mean'})

%% Individuazione delle righe in cui la data di nascita è precedente il 31/12/1967
boo=X{:,2}<'31-Dec-1967';
X1=Xall(boo,:);
% la procedura si ripete su X1