function [dret] = CohensD( distroA, distroB ) n1 = numel( distroA ); n2 = numel( distroB ); s = sqrt( ( (n1 - 1)*(std(distroA)^2) + (n2 - 1)*(std(distroB)^2) ) / (n1 + n2 - 2) ); dret = ( mean(distroA) - mean(distroB) ) / s; end