function [MSE, RMSE, NRMSE, BIAS] = CET_SkillStatistics (xcalc, xmeas)

% GENERAL QUANTITIES
x1m = nanmean(xcalc);
x2m = nanmean(xmeas);

lx1 = length(xcalc);
lx2 = length(xmeas);
ly = lx1;

% MEAN SQUARE ERROR
y = xcalc - xmeas;
y = y.^2;
MSE = sum(y) / ly;

% ROOT MEAN SQUARE ERROR
RMSE = sqrt(MSE);

% NORMALISED ROOT MEAN SQUARE ERRO
NRMSE = RMSE / x2m * 100;

% BIAS
BIAS = x2m - x1m;

return