perl - 关于来自 PDL::GSL::CDF (Perl) 的 gsl_cdf_tdist_P

Question

当 p 值小于 1E-16 时，您能否告诉我如何从 gsl_cdf_tdist_P 函数中获取完整的 p 值？我得到的是 0。

谢谢，伍迪

print "t-test p-value = " . ttest(\@n,\@t) . "\n";

sub ttest{

    my ($n,$t) = @_;
    my @n = @$n;
    my @t = @$t;
    my $nn = pdl(@n);
    my $tt = pdl(@t);
    my ($tstats, $df) = t_test( $nn, $tt );
    use PDL::GSL::CDF; 
    my $p_2tail = 2 * (1 - gsl_cdf_tdist_P( $tstats->abs, $df ));
    return $p_2tail;
}

我的输入值如下：

my @n = qw (1 2 4 2 3 1 2 4 2 1 2 4 2 3 1 2 4 2 1 2 4 2 3 1 2 4 2);
my @t = qw (11 12 13 12 13 11 14 11 12 13 12 13 11 14 11 12 13 12 13 11 14);

Answer 1

我为这个问题找到了一个简单的解决方案。我曾经gsl_cdf_tdist_Q得到p值。p 值不限于小数点后 16 位，因为它们不是重新映射的 p 值，例如(1-gsl_cdf_tdist_P).

伍迪