%% Importazioni dati in formato table 
% (la prima colonna contiene i nomi delle
% righe) e la prima riga i nomi delle variabili.
% Osservazione: non è necessario specificare 'ReadVariableNames',true
% poiché questo è già il default. Al contrario il default di ReadRowNames è
% false di conseguenza è necessario settarlo a true
Xtable=readtable('Ind_ec.xlsx','Sheet','tutte_le_unità','Range','A1:I15','ReadRowNames',true,...
    'ReadVariableNames',true);

Xdouble=table2array(Xtable(:,4:end));

% Definisco la variabile peso (gli addetti)
% w=table2array(Xtable(:,1));
w=Xtable.ADDETTI;

% Calcolo delle medie ponderate
meanweighted=sum(Xdouble.*w)/sum(w);

% Calcolare gli scostamenti quadratici medi ponderati 
stdweighted=sqrt(sum(((Xdouble-meanweighted).^2).*w)/sum(w));

%% Display results
% Mostro i risulati utilizzando le table
disp('Medie ponderate')
meanT=array2table(meanweighted,'VariableNames',Xtable.Properties.VariableNames(4:end));
disp(meanT)

disp('Standard deviation ponderate')
stdT=array2table(stdweighted,'VariableNames',Xtable.Properties.VariableNames(4:end));
disp(stdT)

%% Esclusione del settore delle raffinerie
X1table=Xtable;
% Cancello la riga riferita alle raffinerie
X1table('raffinerie',:)=[];
X1double=table2array(X1table(:,4:end));
w=table2array(X1table(:,1));

% Calcolo delle medie ponderate
meanweighted=sum(X1double.*w)/sum(w);

% Calcolare gli scostamenti quadratici medi ponderati 
stdweighted=sqrt(sum(((X1double-meanweighted).^2).*w)/sum(w));

% Display results
% Mostro i risulati utilizzando le table
disp('Medie ponderate')
meanT=array2table(meanweighted,'VariableNames',Xtable.Properties.VariableNames(4:end));
disp(meanT)

disp('Standard deviation ponderate')
stdT=array2table(stdweighted,'VariableNames',Xtable.Properties.VariableNames(4:end));
disp(stdT)