correlpvalue.m

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

function p = correlpvalue(r, n)
% correlpvalue  - pvalue of correlation coefficient
%
% FORMAT:       p = correlpvalue(r, n)
%
% Input fields:
%
%       r           correlation values (0 ... 1)
%       n           number of points in correlated series
%
% Output fields:
%
%       p           alpha probability of false alarm for significance
%
% See also correltstat.

% Version:  v0.7a
% Build:    7082110
% Date:     Aug-21 2007, 10:19 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

% compute p value from 2-tailed t-statistic
p = 2 - 2 .* custom_tcdf(correltstat(abs(r), n), n - 2);