As noted earlier, the data in
Figure 5indicate that thinning rates of both the PCTF and the PLTF are often in a rather narrow range of approximately 0 to 2 μm/min, with a considerably broader distribution of more rapid thinning rates.
Figures 9A and 9Bshow histograms of thinning rates for the PCTF and PLTF, respectively. These histograms are obviously asymmetrical, with narrow peaks corresponding to slow thinning of approximately 1 μm/min but with many instances of much more rapid thinning. PDFs were fitted to these histograms by maximizing likelihood (
L) defined by
\[L{=}{{\prod}_{i{=}1}^{N}}\ p(r_{i}),\]
where
r i is the rate of thickness change for the
ith of
N recordings (
N = 80 for the PCTF, 76 of the PLTF), and
p(r i) is the PDF equation to be fitted to the histogram. The following two PDFs were fitted: first, a bimodal PDF, being the sum of two Gaussian functions
\[p(r){=}(f_{1}/\sqrt{2{\pi}}{\sigma}_{1})e^{{-}(r{-}r_{1})^{2}/2{\sigma}_{1}^{2}}{+}{[}(1{-}f_{1})/\sqrt{2{\pi}}{\sigma}_{2}{]}e^{{-}(r{-}r_{2})^{2}/2{\sigma}_{2}^{2}},\]
where
r is rate of thickness change, and the parameters
f 1, σ
1, r 1, σ
2, and
r 2 were adjusted to maximize the likelihood. The second type of fit was a skewed PDF with a single peak (at
r =
r 0) of the form
\[p(r){=}(f/\sqrt{2{\pi}}{\sigma})e^{{-}(r{-}r_{0})^{2}/2{\sigma}^{2}}{+}{[}(1{-}f)k/(ae^{b(r{-}r_{0})}{+}be^{{-}a(r{-}r_{0})}){]},\]
where the parameters
f, σ
, r 0 , a, and
b were adjusted to maximize likelihood. The parameter
k was adjusted so that, for any values of
a and
b \[{{\int}_{{-}{\infty}}^{{\infty}}}k/(ae^{b(r{-}r_{0})}{+}be^{{-}a(r{-}r_{0})})dr{=}1.\]
Both fits involve the same number (
n = 5) of adjustable parameters, all of which have considerable impact on the final fit. For the PCTF, the skewed PDF of equation gave a maximum likelihood approximately eight times greater than that for the bimodal PDF of
equation 3 , and the former PDF is given by the curve in
Figure 9A . The peak of this PDF,
r 0, corresponding to slow thinning, occurred at a thinning rate of 0.79 μm/min; rapid thinning corresponds to the long tail on the left of the PDF. For the PLTF, the bimodal PDF of
equation 3gave a maximum likelihood approximately 4000 times greater than that for the skewed PDF of
equation 4 , and the former PDF is given by the curve in
Figure 9B . The peaks of the two Gaussian functions correspond to a slow thinning rate (
r 1) of 1.25 μm/min and a rapid thinning rate (
r 2) of 9.19 μm/min. The probability of slow thinning (
f 1 in
equation 3 ) was 0.30. (In considering the relative likelihood of bimodal and skewed PDFs, it should be noted that the rate of thickness change data,
r i in
equation 2 , are not all independent data but are obtained from four recordings from each of 20 subjects. This tends to overestimate the likelihood ratios. Also, the likelihood-ratio estimates, of course, depend on the exact form of the assumed PDF equations.)