%% 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
X=readtable('Ind_ec.xlsx','Sheet','tutte_le_unità','Range','A1:I15','ReadRowNames',true,...
    'ReadVariableNames',true,'Format','auto');

Xdouble=table2array(X(:,4:end));
w=table2array(X(:,1));

% 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',X.Properties.VariableNames(4:end));
disp(meanT)

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

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

% 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',X.Properties.VariableNames(4:end));
disp(meanT)

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