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
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
// =============================================================================
// Tests for cdfchn function
// =============================================================================
// Run with test_run('statistics','cdfchn',['no_check_error_output']);
prec = 1.e-5;
x = [0.01,0.025,0.05,0.1,0.9,0.95];
y1 = cdfchi("PQ",x,ones(x));
y2 = cdfchn("PQ",x,ones(x),0*ones(x));
if norm(y1-y2) > prec then bugmes();quit;end
df = [1,2,3,4,5,6];
pno = df/10;
[P,Q] = cdfchn("PQ",x,df,pno);
[x1] = cdfchn("X",df,pno,P,Q);
[df1] = cdfchn("Df",pno,P,Q,x);
if norm(x1-x) > prec then bugmes();quit;end
if norm(df1-df) > prec then bugmes();quit;end
// IEEE support
// See http://bugzilla.scilab.org/show_bug.cgi?id=7296
Df = 1;
Pnonc = 0.1;
X = %inf; // Inf
[P,Q] = cdfchn("PQ", X, Df, Pnonc);
assert_checkequal(P, 1);
assert_checkequal(Q, 0);
X = %nan; // NaN
[P,Q] = cdfchn("PQ", X, Df, Pnonc);
assert_checkequal(P, %nan);
assert_checkequal(Q, %nan);
|
