Fixes about computing pair-wise and multi-correlations for the Hessian case.
parent
574f94eafa
commit
1e5e214503
|
@ -65,14 +65,15 @@ if hess_flag==0,
|
||||||
McoJ(ind1(ii),:) = cosn([JJ1(:,ii),JJ1(:,find([1:1:size(JJ1,2)]~=ii))]);
|
McoJ(ind1(ii),:) = cosn([JJ1(:,ii),JJ1(:,find([1:1:size(JJ1,2)]~=ii))]);
|
||||||
end
|
end
|
||||||
else
|
else
|
||||||
deltaJ = sqrt(diag(JJ));
|
deltaJ = sqrt(diag(JJ(ind1,ind1)));
|
||||||
tildaJ = JJ./((deltaJ)*(deltaJ'));
|
tildaJ = JJ(ind1,ind1)./((deltaJ)*(deltaJ'));
|
||||||
McoJ(:,1)=(1-1./diag(inv(tildaJ)));
|
McoJ(ind1,1)=(1-1./diag(inv(tildaJ)));
|
||||||
rhoM=sqrt(1-McoJ);
|
rhoM=sqrt(1-McoJ);
|
||||||
% PcoJ=inv(tildaJ);
|
% PcoJ=inv(tildaJ);
|
||||||
PcoJ=inv(JJ);
|
PcoJ=NaN(npar,npar);
|
||||||
|
PcoJ(ind1,ind1)=inv(JJ(ind1,ind1));
|
||||||
sd=sqrt(diag(PcoJ));
|
sd=sqrt(diag(PcoJ));
|
||||||
PcoJ = PcoJ./((sd)*(sd'));
|
PcoJ = abs(PcoJ./((sd)*(sd')));
|
||||||
end
|
end
|
||||||
|
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue