EOF_STAT
EOF_STAT(A,FRAC_TIMESER) Used with EOF_SPACE and/or EOF_TFUNC. Return statistics related to an EOF solution for a given set of parameters. Results are on the x-axis j = 1: # EOFscomputed and scaled, j = 2: % percentage of total variance accounted for by each eigenvector, j = 3: the eigenvalues.
|
Arguments: |
A |
Variable in any spatial dimensions, and time |
|
FRAC_TIMESER |
Use only those time series with this fraction valid data, e.g. 0.8 to require that 80% of the data be present to use the data at a location. |
|
|
Result Axes: |
X |
ABSTRACT: 1 to NEOF |
|
Y |
ABSTRACT: 1 through 3 as outlined in the description. |
|
|
Z |
NORMAL (no axis) |
|
|
T |
NORMAL (no axis) |
Please see the discussion under EOF_SPACE, and see the demonstration ef_eof_demo.jnl for examples of this function.
Example results:
For a simple sample function, eof_stat called to decompose it into eigenfunctions. We allow data to be used if the time series at the point has at least 80% valid data.
Request the number of eigenvalues computed for this choice of parameters.
yes? list/i=1/j=1 eofstat
VARIABLE : EOF_STAT(SST[X=67W:1W,Y=11S:11N], 0.8)0
DATA SET : COADS Monthly Climatology (1946-1989)
FILENAME : coads_climatology.des
FILEPATH : /home/ja9/tmap/fer_dsets/descr/
X : 1
Y : 1
284.0
Now get the percent variance explained by the eigenfunctions which were computed.
yes? list/i=1:10/j=2 eofstat
VARIABLE : EOF_STAT(SST[X=67W:1W,Y=11S:11N], 0.8)
DATA SET : COADS Monthly Climatology (1946-1989)
FILENAME : coads_climatology.des
FILEPATH : /home/ja9/tmap/fer_dsets/descr/
SUBSET : 10 points (X)
Y : 2
2
2
1 / 1: 86.95
2 / 2: 5.82
3 / 3: 3.87
4 / 4: 1.51
5 / 5: 0.56
6 / 6: 0.38
7 / 7: 0.31
8 / 8: 0.23
9 / 9: 0.15
10 / 10: 0.11
And finally the eigenvalues associated with these eigenfunctions.
yes? list/i=1:10/j=3 eofstat
VARIABLE : EOF_STAT(SST[X=67W:1W,Y=11S:11N], 0.8)
DATA SET : COADS Monthly Climatology (1946-1989)
FILENAME : coads_climatology.des
FILEPATH : /home/ja9/tmap/fer_dsets/descr/
SUBSET : 10 points (X)
Y : 3
3
3
1 / 1: 249.4
2 / 2: 16.7
3 / 3: 11.1
4 / 4: 4.3
5 / 5: 1.6
6 / 6: 1.1
7 / 7: 0.9
8 / 8: 0.7
9 / 9: 0.4
10 / 10: 0.3