correltstat.m

来自「toolbox of BVQX, This is the access betw」· M 代码 · 共 47 行

function t = correltstat(r, n)
% correltstat  - convert correlation value into t-statistic
%
% FORMAT:       t = correltstat(r, n)
%
% Input fields:
%
%       r           correlation values (0 ... 1)
%       n           number of points in correlated series
%
% Output fields:
%
%       t           matching t-statistic

% Version:  v0.7a
% Build:    7082110
% Date:     Aug-21 2007, 10:13 AM CEST
% Author:   Andrew Holmes, SPM2
% URL/Info: http://www.fil.ion.ucl.ac.uk/spm/

% argument check
if nargin < 2 || ...
   ~isa(r, 'double') || ...
   ~isa(n, 'double') || ...
    isempty(r) || ...
    isempty(n) || ...
    any(isnan(r(:)) | isinf(r(:)) | abs(r(:)) > 1) || ...
    any(isinf(n(:)) | isnan(n(:)) | n(:) < 2)
    error( ...
        'BVQXtools:BadArgument', ...
        'Missing or invalid argument given.' ...
    );
end
osize = size(r);
if numel(r) == 1 && ...
    numel(n) > 1
    osize = size(n);
end
r = r(:);
n = n(:);

% delete invalid entries
r(r == 1) = 0;

% formula
t = reshape(r .* sqrt((n - 2) ./ (1 - r .^ 2)), osize);