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

Last change on this file since 182 was 182, checked in by pinsard, 12 years ago

fix some svn propset
  • Property svn:keywords set to Id URL
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;  @project_init
25;  x = ++
26;  y = ++
27;  sig_value = ++
28;  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; $URL$
45;
46; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
47;
48;   * minimal header
49;
50;-
51pro significance_test_99, x, y, sig_value
52;
53   X=x
54   Y=y
55;
56   X=X-total(X)/n_elements(X)
57   Y=Y-total(Y)/n_elements(Y)
58    sx=stdev(X)
59    sy=stdev(Y)
60    reg=correlate(X,Y)*sy/sx
61    sse=total((Y-reg*X)^2)/n_elements(X)
62    t=abs(sqrt(n_elements(X)-2.)*reg/sqrt(sse/sx^2))
63    p=0.01 ;(for 99% significant)
64;    tc=t_cvf(p/2.,n_elements(X)-1)
65    tc=t_cvf(p/2.,n_elements(X)/5.-1)
66;
67    sig1=1.*(t ge tc)
68if (sig1 eq 1. ) then begin
69   print, 'Statistics are valid at 99% significance level'
70endif else begin
71   print, 'Statistics are not valid at 99% significance level'
72endelse
73;
74end
Note: See TracBrowser for help on using the repository browser.