%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Compute the sample (bias-corrected) skewness. 
% 
% CALL 
% 
% p = Skew(x);
% 
% VERSION:     $Id: Skew.m,v 1.1 2011/02/25 10:04:07 congedo Exp $
%
% HISTORY:     25-02-2011 G. Congedo
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function s = Skew(x)

n = numel(x);
mu3 = sum((x-mean(x)).^3)/n;
sig = sqrt(sum((x-mean(x)).^2)/n);
s = mu3/(sig)^3;

% Bias correction
s = sqrt(n*(n-1))/(n-2)*s;

end