%Reduction in error statistic. For the MBH98 and MM05 studies, y_hat is the
%proxy or simulated obseravtions and y_ref is the instrumental record.

function re=rePH(y_hat,y_ref);

re=1-sum( (y_hat(:)-y_ref(:)).^2 )/sum(y_ref.^2);