1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
function output = MakeGammaMonotonic(input)
% output = MakeGammaMonotonic(input)
%
% Make input monotonically increasing from lowest value.
%
% This version also forces the last value to be 1, and then
% enforces decreasing downwards from there, after the
% enforcement of increasing from 0.
%
% 3/1/99 dhb Handle multiple columns.
% 8/03/07 dhb Old routine just enforced non-decreasing. Fixed to make strictly increasing.
% 3/07/10 dhb Wrote from MakeMonotonic.
% Did not want to change behavior of MakeMonotonic in case it is called from
% programs unrelated to gamma fitting.
% 3/08/10 dhb Actually do the enforcement of 1.
% 5/27/10 dhb Use a larger bump, 100*eps
[m,n] = size(input);
output = input;
for j = 1:n
for i = 1:m-1
if (output(i,j) >= output(i+1,j))
output(i+1,j) = output(i,j)+100*eps;
end
end
output(m,j) = 1;
for i = m:-1:2
if (output(i,j) <= output(i-1,j))
output(i-1,j) = output(i,j)-100*eps;
end
end
end
|
