Incorect PGD Index equation
Opened this issue · 1 comments
dervinism commented
Are you using the correct PGD Index equation in the code? The one that is in the paper is different. It seems that the following lines in the GradientMethod function
Mag_of_MeanGrad = sqrt(sum(nansum(Gx(:,:,TrailTime)))^2+sum(nansum(Gy(:,:,TrailTime)))^2);
Mean_of_MagGrad = sum(nansum(sqrt(Gx(:,:,TrailTime).^2+Gy(:,:,TrailTime).^2)));
PGD(TrailTime) = Mag_of_MeanGrad/Mean_of_MagGrad;
should be replaced with this line
PGD(TrailTime) = (sum(nansum(Gx(:,:,TrailTime).*Gy(:,:,TrailTime))))/sqrt(sum(nansum(Gx(:,:,TrailTime).^2))*sum(nansum(Gy(:,:,TrailTime).^2)));
OccassionalBioinformatician commented
@dervinism Sorry? The equation of PGD is given in Rubino (2006; p. 1555) and the original code is accurate.
Zabeh, et al. (2022; p. 10) incorrectly labels the Pearson Correlation Coefficient as PGD. For what it's worth, I do not see any code in this package that allows us to evaluate the goodness of fit by such a calculation; we need to do it by ourselves.