Step 2: Assuming nothing is known about \theta, a uniform prior in the interval [0,1] is employed.
Step 4: The distance between the observed and simulated transition frequencies \rho(\omega_{S,i}, \omega_E) = |\omega_{S,i} - \omega_{E}| is computed for all parameter points (Table 1, column 5).