/usr/share/psychtoolbox-3/PsychGamma/FitGammaPoly.m is in psychtoolbox-3-common 3.0.14.20170103+git6-g605ff5c.dfsg1-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 [fit_out,x,err] = FitGammaPoly(values_in,measurements,values_out)
% [fit_out,x,err] = FitGammaPoly(values_in,measurements,values_out)
%
% Fit homogeneous polynomial to gamma data.
%
% 3/15/94 dhb, jms Ignore low values on fit.
% 3/4/05 dhb Removed old commented out code.
% Allocate return space
[mOut,nOut] = size(values_out);
fit_out = zeros(mOut,nOut);
% Find first non-zero entry
[mIn,nIn] = size(values_in);
index = find(measurements ~= 0.0);
if (~isempty(index))
useIndex = index(1);
useIn = values_in(useIndex:mIn,1);
useMeas = measurements(useIndex:mIn,1);
% Get initial values
x0 = InitialXPoly(useIn,useMeas);
x = x0;
% Compute fit values and error to data for return
zeroThresh = values_in(useIndex);
fitOutIndex = find(values_out >= zeroThresh);
fitInIndex = find(values_in >= zeroThresh);
fit_out(fitOutIndex) = ComputeGammaPoly(x,values_out(fitOutIndex));
fit_in = ComputeGammaPoly(x,values_in(fitInIndex));
err = ComputeFSSE(fit_in,measurements(fitInIndex));
else
x = [];
err = 0.0;
end
|