[43] | 1 | ;+ |
---|
[46] | 2 | ; .. _statistics.pro: |
---|
[43] | 3 | ; |
---|
| 4 | ; ============== |
---|
| 5 | ; statistics.pro |
---|
| 6 | ; ============== |
---|
| 7 | ; |
---|
| 8 | ; DESCRIPTION |
---|
| 9 | ; =========== |
---|
| 10 | ; |
---|
| 11 | ; SEE ALSO |
---|
| 12 | ; ======== |
---|
| 13 | ; |
---|
[50] | 14 | ; :ref:`project_profile.sh` |
---|
[43] | 15 | ; |
---|
[46] | 16 | ; :ref:`as_validation_net_1994_95_v10.pro` |
---|
| 17 | ; :ref:`coare_validation_net_1992_93_all_v10.pro` |
---|
| 18 | ; :ref:`stratus_validation_net_2001_2002_all_v10.pro` |
---|
| 19 | ; :ref:`subdctn_cs_validation_net_1991_93_all_v10.pro` |
---|
| 20 | ; :ref:`subdctn_ses_validation_net_1991_93_all_v10.pro` |
---|
| 21 | ; :ref:`subdctn_sws_validation_net_1991_93_all_v10.pro` |
---|
| 22 | ; :ref:`whots_validation_net_2004_2007_all_v10.pro` |
---|
| 23 | ; :ref:`icoads_q2m_stats_paper.pro` |
---|
| 24 | ; :ref:`icoads_sst_stats_paper.pro` |
---|
| 25 | ; :ref:`icoads_t2m_stats_paper.pro` |
---|
| 26 | ; :ref:`icoads_ws_stats_paper.pro` |
---|
| 27 | ; |
---|
[43] | 28 | ; EXAMPLES |
---|
| 29 | ; ======== |
---|
| 30 | ; |
---|
| 31 | ; :: |
---|
| 32 | ; |
---|
[180] | 33 | ; moor = ++ |
---|
| 34 | ; ++ |
---|
| 35 | ; statistics, moor, erai, cor, bias, std, rmsd |
---|
[43] | 36 | ; |
---|
[46] | 37 | ; TODO |
---|
| 38 | ; ==== |
---|
| 39 | ; |
---|
| 40 | ; coding rules |
---|
| 41 | ; |
---|
| 42 | ; complete description |
---|
| 43 | ; |
---|
| 44 | ; complete example |
---|
| 45 | ; |
---|
[43] | 46 | ; EVOLUTIONS |
---|
| 47 | ; ========== |
---|
| 48 | ; |
---|
| 49 | ; $Id$ |
---|
| 50 | ; |
---|
[85] | 51 | ; $URL$ |
---|
| 52 | ; |
---|
[43] | 53 | ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) |
---|
| 54 | ; |
---|
| 55 | ; * minimal header |
---|
| 56 | ; |
---|
| 57 | ;- |
---|
| 58 | ; |
---|
[41] | 59 | pro statistics, moor, erai, $ |
---|
| 60 | cor, bias, std, rmsd |
---|
[50] | 61 | @cm_project |
---|
[97] | 62 | ; mean1= mean of mooring variable std1=std of mooring variable |
---|
| 63 | ; mean2= mean of erai variable std2=std of erai variable |
---|
| 64 | x=moor |
---|
[94] | 65 | y=erai |
---|
[97] | 66 | ind1=where(finite(x,/nan)) |
---|
[94] | 67 | ind2=where(finite(y,/nan)) |
---|
[41] | 68 | ind=inter(ind1,ind2) |
---|
| 69 | if (ind ne -1 ) then begin |
---|
[97] | 70 | x(ind)=!Values.f_nan |
---|
[94] | 71 | y(ind)=!Values.f_nan |
---|
[41] | 72 | endif |
---|
[97] | 73 | ; |
---|
[41] | 74 | ind=where(finite(x)) |
---|
[97] | 75 | x=x(ind) |
---|
[94] | 76 | y=y(ind) |
---|
[97] | 77 | ; |
---|
[41] | 78 | nn=n_elements(ind) |
---|
[97] | 79 | ; |
---|
[41] | 80 | cor=correlate(x,y) |
---|
| 81 | xxx=moment(x) |
---|
[97] | 82 | mean1=xxx(0) |
---|
[94] | 83 | std1=sqrt(xxx(1)) |
---|
[41] | 84 | xxx=moment(y) |
---|
[97] | 85 | mean2=xxx(0) |
---|
[94] | 86 | std2=sqrt(xxx(1)) |
---|
[97] | 87 | bias=mean2-mean1 |
---|
[94] | 88 | std=std2/std1 |
---|
[41] | 89 | rmsd=sqrt(total((x-y)*(x-y))/n_elements(x)) |
---|
[97] | 90 | ; |
---|
[41] | 91 | end |
---|