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
|
function [f,g] = WeibAlphTAFCFitFun(x,beta,level,nCorrect,nError)
% [f,g] = WeibAlphTAFCFitFun(x,,beta,level,nCorrect,nError)
%
% 8/26/94 dhb Wrote it.
% Unpack vector x
alpha = x(1);
% Force really big error if alpha or beta <= 0
if (alpha <= 0)
f = 1e25;
g = -1;
return;
end
% Compute Weibull function on input levels
pCorrect = ComputeWeibTAFC(level,alpha,beta);
% Handle range problem, can't take log(0);
tol = 1e-4;
o_index = find(pCorrect == 1);
if (~isempty(o_index))
pCorrect(o_index) = (1-tol)*ones(length(o_index),1);
end
% Compute error function, which is -log(likelihood).
tmp = nCorrect.*log(pCorrect) + nError.*log(1 - pCorrect);
f = -sum(tmp);
g = -1;
|
