Correction: normrnd in script with no correlation

script_workflow_PcTK_NoCorrelation.m has been updated to version 1.1. The new file can be downloaded from here and will be included in the package version 3.22.

There was a mistake on how to use a Matlab function normrnd(mu, sigma). The ‘sigma’ must be the standard deviation of data, sqrt(m2_yt), not the variance of data, m2_yt. The corrected code in lines 260–268 in version 1.1 appears:

•             %   Step 5. Generate a noisy 4-D count sinogram, nRm(Nl,ch,row,view)
•             switch iflag_noise
•                 case flag_noise_normal  % gen Gaussian data
•                     m2_yn = normrnd(m2_yt, sqrt(m2_yt));        % ver 1.1
•                 case flag_noise_Poisson % gen Poisson data
•                     m2_yn = random(‘Poisson’, m2_yt);
•             end
•             m3_yt(:,:,irow) = m2_yt;
•             m3_yn(:,:,irow) = m2_yn;

We thank Danyang Li of The Southern Medical University for drawing our attention to it.