Sorry, I'm withdrawing my question.

I found that correct replication should be:

ROI1 = [1 3 1 23 4 2 0 -2 -20 4 2 9 21]';
ROI2 = [2 3 0 15 2 0 -5 -3 -10 2 0 5 -5]';
W = [1 1 0.5 -0.1 0 1 1 0 -1 1 1 1 0]'; W = abs(W);

b12 = lscov([ROI2 ones(13,1)],ROI1,W);
b21 = lscov([ROI1 ones(13,1)],ROI2,W);
r12 = b12*std(ROI2,W)/std(ROI1,W);
r21 = b21*std(ROI1,W)/std(ROI2,W);
Z12 = atanh(r12);
Z21 = atanh(r21);