source: trunk/src/paper01/fig4/significance_test_99.pro @ 70

Last change on this file since 70 was 50, checked in by pinsard, 13 years ago

replace TROPFLUX by PROJECT
File size: 1.2 KB
Line 
1;+
2; .. _significance_test_99.pro:
3;
4; ========================
5; significance_test_99.pro
6; ========================
7;
8; DESCRIPTION
9; ===========
10;
11; SEE ALSO
12; ========
13;
14; :ref:`project_profile.sh`
15;
16; :ref:`fig10_swr_correction.pro`
17; :ref:`fig4_met_var_correction_scatter.pro`
18;
19; EXAMPLES
20; ========
21;
22; ::
23;
24;  IDL> @project_init
25;  IDL> x=++
26;  IDL> y=++
27;  IDL> sig_value=++
28;  IDL> significance_test_99, x, y, sig_value
29;
30; TODO
31; ====
32;
33; coding rules
34;
35; complete description
36;
37; complete example
38;
39; EVOLUTIONS
40; ==========
41;
42; $Id$
43;
44; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
45;
46;   * minimal header
47;
48;-
49pro significance_test_99, x, y, sig_value
50
51   X=x & Y=y
52
53   X=X-total(X)/n_elements(X)
54   Y=Y-total(Y)/n_elements(Y)
55    sx=stdev(X)
56    sy=stdev(Y)
57    reg=correlate(X,Y)*sy/sx
58    sse=total((Y-reg*X)^2)/n_elements(X)
59    t=abs(sqrt(n_elements(X)-2.)*reg/sqrt(sse/sx^2))
60    p=0.01 ;;;(for 99% significant)
61;    tc=t_cvf(p/2.,n_elements(X)-1)
62    tc=t_cvf(p/2.,n_elements(X)/5.-1)
63
64    sig1=1.*(t ge tc)
65if (sig1 eq 1. ) then begin
66   print, 'Statistics are valid at 99% significance level'
67endif else begin
68   print, 'Statistics are not valid at 99% significance level'
69endelse
70
71end
Note: See TracBrowser for help on using the repository browser.