/usr/share/psychtoolbox-3/Psychometric/ComputeWeibYN.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 37 38 39 40 41 42 43 | function pyes = ComputeWeibYN(x,alpha,beta)
% pyes = ComputeWeibYN(x,alpha,beta)
%
% Compute the value of a yes-no weibull psychometric function.
% If 'x' is a matrix and the parameters are scalars, the same
% psychometric function is applied to the whole matrix.
% If 'x' is a matrix and the parameters are vectors, the
% parameters are taken as applying to separate columns.
%
% pyes = ( 1.0 - exp( - (x./alpha).^beta ) )
%
% 9/15/93 jms Added some special casing to deal with matrix inputs
% with/without vectors of parameters
% Check dimensions of parameters
[nRow,nAlpha] = size(alpha);
if (nRow ~= 1)
error('alpha must be a row vector');
end
[nRow,nBeta] = size(beta);
if (nRow ~= 1)
error('beta must be a row vector');
end
[nInputs, nCol] = size(x);
% Case: nAlpha == nBeta == 1 ~= nCol
if (nAlpha == nBeta && nAlpha == 1)
% disp('Single parameter case');
z = (x./alpha).^beta;
pyes = ( 1.0 - exp( - z ) );
% Case: nAlpha == nBeta == nCol
elseif (nAlpha == nBeta && nAlpha == nCol)
% disp('Vector parameter case');
for (i=1:nCol)
z = (x(:,i)./alpha(i)).^beta(i);
pyes(:,i) = ( 1.0 - exp( - z ) );
end
else
error('Number of input columns not equal to number of parameters');
end
|