function [X,moya,ecaa] = cenred(Z,moya,eca)
% Fonction qui centree par moya et eca 
%avec la formule X=(2/3)*(Z-moya)/eca
%les invivudus etant en ligne
UN = ones(size(Z));
Moy =moya ;
Std =eca;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Moy = diag(Moy);
Moy = UN * Moy;
taille=size(Z,1);
Ecart_type = diag(sqrt((taille-1)/taille)*Std);
Ecart_type = UN * Ecart_type;
 
 
% Calcul de la matrice X des données centrées réduites
X = (2/3).*(Z - Moy)./Ecart_type ;
