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

%% Implementazione veloce
eti=({'<=1000]'; '(1000-1200]'; '(1200 1500]'; '(1500-2200]'; '>2200'});

Classi=[-Inf 1000 1200 1500 2200 Inf];
[freq,edges,bin]=histcounts(retr,Classi);
% Il primo ouptut di histcounts (freq) contiene le frequenze in
% corrispondenza di ogni classe
% Il secondo argomento edges contiene gli estremi delle classi (in questo
% caso edges=Classi)
% bin è un vettore che ha la stessa lunghezza di retr che contiene a quale
% classe è stata assegnato ogni elemento di retr. bin assume valori 1, 2,
% ....

% retrclassi è una cell che contiene a quale classe di retribuzione
% appartiene ogni unità statistica
retrclassi=eti(bin);
freqvettorecolonna=freq';
DistFreqtable=array2table(freqvettorecolonna,'RowNames',eti,'VariableNames',{'Distr Frequenza'});
disp(DistFreqtable)


% %% Creazione della distribuzione di frequenze e delle etichette per ogni classe
% % FACOLTATIVO
% % Inizializzazione del cell array che contiene in corrispondeza delle riga
% % i-esima l'etichetta della classe di appartenenza
% retrclassi1=cell(length(retr),1);
% 
% freq=zeros(5,1);
% 
% boo=retr<=1000;
% retrclassi1(boo)=eti(1);
% freq(1)=sum(boo);
%     
% boo=retr>=1000 & retr<=1200;
% retrclassi1(boo)=eti(2);
% freq(2)=sum(boo);
% 
% boo=retr>=1200 & retr<=1500;
% retrclassi1(boo)=eti(3);
% freq(3)=sum(boo);
% 
% boo=retr>1500 & retr<=2200;
% retrclassi1(boo)=eti(4);
% freq(4)=sum(boo);
% 
% boo=retr>2200;
% retrclassi1(boo)=eti(5);
% freq(5)=sum(boo);
% 
% Df=array2table(freq,'RowNames',eti,'VariableNames',{'Distr frequenze'});
% 
% %% Implementazione del vettore retrclassi in maniera avanzata
% % FACOLTATIVO
% % Si cicla su tutte le righe di retr e si trova (tramite find) la prima
% % volta in cui retr(i) è minore di Classi(2:end)
% % Ad es. se retr(i)<Classi(2) allora si ha a che fare con una retribuzione
% % minore di 1000 Euro, in questo caso il risultaton di find è 1
% retrclassi2=retrclassi1;
% for i=1:length(retr)
%     trovaclasse=find(retr(i)<Classi(2:end),1);
%         retrclassi2(i)=eti(trovaclasse);
% end
% 
% disp('Controllo uguaglianza tra retrclassi2 e retrclassi1')
% isequal(retrclassi2,retrclassi)


%% Creazione della tabella a doppia entrata tra retrclassi e sesso

% L'istruzione categorical 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);
disp('Tabella a doppia entrata tra le variabili Retribuzione in classi e sesso')
disp(TB)

%% Analisi della media della retribuzione per sesso e classe di retribuzione
% Xalltable è table che contiene 3 colonne
% Sesso, RetribuzioneEuro e classi di retribuzione
Xalltable=[X(:,[1 3]) table(retrclassic)];

% Osservazione: un modo alternativo per aggiungere una variabile ad una
% table esistente è quello di chiamare la funzione addvars
Xalltablechk=addvars(X(:,[1 3 ]),retrclassic);


% Calcolo della media della retribuzione per ciascuna combinazione di Sesso
% e classe di retribuzione
A=grpstats(Xalltable,{'Sesso' 'retrclassic'},{'mean'});

%% Individuazione delle righe in cui la data di nascita è precedente il 31/12/1967
boo=X{:,2}<'31-Dec-1967';
X1alltable=Xalltable(boo,:);
% la procedura si ripete su X1alltable
% calcolo della media della retribuzione per ciascuna combinazione di Sesso
% e classe di retribuzione
grpstats(X1alltable,{'Sesso' 'retrclassic'},{'mean'})