% FATHOM: A Matlab Toolbox for Ecological & Oceanographic Data Analysis
% by Dave Jones<djones@rsmas.miami.edu>               
% http://www.rsmas.miami.edu/personal/djones
%
%   f_AIC              - compute AIC (or BIC) for least-squares based methods
%   f_ancova           - nonparametric (permutation-based) 1-way ANCOVA
%   f_anosim           - 1-way Analysis of Similarity (ANOSIM)
%   f_anosim2          - 2-way crossed ANOSIM with no replication
%   f_arrow            - Draw an arrow.
%   f_balance          - get indices of elements to create a balanced design
%   f_bartlett         - Bartlett's test for Homogeneity of Variances
%   f_bioenv           - correlation between distance matrix and all possible subsets of x
%   f_biplot           - eigenvector-based 2-d distance biplot
%   f_biplotEnv2       - create environmental vectors for 2-d nMDS ordination distance biplot
%   f_biplotEnv3       - create environmental vectors for 3-d nMDS ordination distance biplot
%   f_biplotPca2       - create 2-d PCA distance biplot
%   f_biplotSpecies    - create species vectors for ordination distance biplot
%   f_boot             - bootstrap resampling with replacement
%   f_braycurtis       - calculate Bray-Curtis symmetric distance matrix
%   f_brokenstick      - determine # of significant ordination dimensions via "Broken-Stick model"
%   f_cap              - Canonical Analysis of Principal Coordinates using ANY distance matrix
%   f_capFindX         - Variable selection for CCorA via backward elimination
%   f_capFindY         - Variable selection for CAP-based CDA via backward elimination
%   f_capScale         - Canonical Analysis of Principal Coordinates using ANY distance matrix
%   f_capValid         - Leave-One-Out Cross-Validation for f_cap
%   f_capVectors       - correlation vector plot for f_cap
%   f_cda              - Canonical Discriminant Analysis
%   f_cda2             - Canonical Discriminant Analysis (legacy version)
%   f_cda632           - estimate CDA error rate using bootstrap .632+ rule
%   f_cdaClass         - classify unknown observations for f_cda
%   f_cdaCV            - leave-one-out cross validation for Canonical Discriminant Analysis
%   f_cdaPlot          - plot results from a Canonical Discriminant Analysis
%   f_cdaValid         - leave-one-out cross-validation for CDA
%   f_cdaVectors       - correlation vector plot for f_cda
%   f_center           - center data on column mean
%   f_centroid         - returns coordinates of the centroid of X, optionally partitioned into groups
%   f_compet           - get index of column containing the maximum value of each row
%   f_confEllip        - bootstrapped confidence ellipse for bivariate data
%   f_confEllipse      - parametric or bootstrapped confidence ellipse for bivariate data
%   f_convHull         - convex hull for bivariate data
%   f_corr             - Pearson's, Spearman's, or Kendall's correlation between 2 vectors
%   f_cov              - returns covariance (= dispersion) matrix from X
%   f_depthCM          - 'Depth of Center of Mass' for plankton survey data
%   f_depthMean        - weighted 'Mean Depth of Occurrence' for plankton survey data
%   f_designMatrix     - create ANOVA design matrix using dummy variables
%   f_distlm           - Nonparametric Multivariate Multiple Regression for ANY distance matrix
%   f_eig              - eigenanalysis of a square matrix
%   f_ekmanDepth       - depth of the Ekman layer
%   f_errRate          - error rate for a classifier
%   f_euclid           - computes Euclidean distance matrix
%   f_export           - write ASCII delimited file.
%   f_exportDods       - export DODS bathymetry for import into Surfer
%   f_exportR          - export data for import into R
%   f_extractFields    - extract structure fields & combine into a single matrix
%   f_figArea          - finds the area of figure objects
%   f_filterMA         - smooth columns of a maxtix using Moving Average
%   f_findCell         - index to rows of Y that match X
%   f_firstOccur       - returns indices of the first occurrence of unique elements of input vector
%   f_getNLOM          - download Navy Layered Ocean Model SSH Nowcast imagery
%   f_gregorian        - convert Julian date to Gregorian
%   f_grpMean          - returns the mean/stdv/se of X (columnwise) separately for groups
%   f_grpRel           - returns the relative proportion of X (column-wise) separately for groups
%   f_gshhs2Shp        - export a shapefile from the GSHHS database for use in Surfer/Arcview
%   f_halfchange       - scale nMDS configuration to half-change
%   f_importSurfer     - import Surfer *.bln blanking file for M_Map
%   f_importSurferBln  - import Surfer *.bln blanking file for M_Map
%   f_importSurferGrd  - import Surfer GRID (*.grd) file
%   f_inv              - matrix inversion via "\" (left division)
%   f_isDST            - determine if dates occur during Daylight Savings Time
%   f_isOdd            - determine if integer is odd or even
%   f_isScalar         - determine if input is a scalar
%   f_issymdis         - determine if input is square symmetric distance matrix
%   f_julian           - convert date vector to Julian date
%   f_julianSub        - get indices to subsample a series of Julian dates
%   f_labelplot        - create 2-d label plot
%   f_latlong          - computes terrestrial distance matrix
%   f_leapYear         - determine if a year is a leap year
%   f_lenfreq          - length-frequency plot (adjusted for sampling effort)
%   f_ll2utm           - convert longitude/latitude coordinates to UTM
%   f_lowpass          - smooth columns of a matrix via a lowpass filter
%   f_mantel           - standardized Mantel statistic for 2 symmetric distance matrices
%   f_modelMatrix      - create model matrix for a Mantel Test
%   f_mregress         - Multiple Linear Regression via Least Squares Estimation
%   f_mst              - Minimum Spanning Tree from a symmetric distance matrix
%   f_multicomb        - generate permutation distribution of grouping labels
%   f_nan2ave          - replace missing values (NaN's) with average value
%   f_nmds             - run Steyvers's nonmetric Multidimensional Scaling program
%   f_nnMLP            - feedforward multilayer perceptron neural network for classification
%   f_nnMLP632         - estimate MLP error rate using bootstrap .632+ rule
%   f_nnMLPcv          - leave-one-out cross validation for f_nnMLP
%   f_normal           - normalize values of X
%   f_npManova         - nonparametric (permutation-based) MANOVA for ANY distance matrix
%   f_npManovaPW       - a posteriori, multiple-comparison tests
%   f_origin           - mark the origin of a graph with horizontal/vertical line(s)
%   f_pca              - Principal Component Analysis of a data matrix
%   f_pcoa             - Principal Coordinates Analysis (PCoA)
%   f_pcoa2            - legacy version of Principal Coordinates Analysis
%   f_pdf              - export current figure to PDF file
%   f_permtest         - two sample permutation test of means
%   f_plotBarsH        - horizontal bar graphs
%   f_plotBarsV        - vertical bar graphs
%   f_plotError        - add error bars to current figure
%   f_plotUSGS         - plot USGS Coastwatch SST image with M_Map Toolbox
%   f_pnn              - Probabilistic Neural Network
%   f_pnn632           - estimate PNN error using bootstrap .632+ rule
%   f_pnnAIC           - stepwise selection of explanatory variables in PNN using AIC (or BIC)
%   f_pnnCV            - leave-one-out cross validation for PNN
%   f_pnnSm            - determine optimal smoothing factor for a PNN
%   f_poly             - create a matrix of polynomials of degree d
%   f_procrustes       - Procrustes rotation of Y to X
%   f_randRange        - returns n random integers ranging from min to max
%   f_randSub          - extract a random subset of N rows from matrix X
%   f_range            - return the min and max values of a vector
%   f_ranging          - scale columns of x to range [0 to 1] or [-1 to 1]
%   f_ranks            - ranks data in x (with averaging of ties)
%   f_rda              - Redundancy Analysis
%   f_rdaAIC           - stepwise selection of explanatory variables in RDA using AIC (or BIC)
%   f_rdaDB            - prepare non-Euclidean distance matrix for db-RDA
%   f_rdaPlot          - ordination distance biplot for a Redundancy Anaysis (RDA)
%   f_rdaStepwise      - stepwise selection of explanatory variables in RDA based on F-stat
%   f_readcwf          - import a Coastwatch Satellite SST file
%   f_recode           - elements of vector as consecutive integers
%   f_rewrap           - wraps vector into symmetric distance matrix (reverses f_unwrap)
%   f_rgb              - utility program for selecting color of plot symbols
%   f_round            - round data to specified number of decimal places
%   f_shadeBox         - shade subsets of a time series plot
%   f_shuffle          - randomly sorts vector, matrix, or symmetric distance matrix
%   f_smooth           - smooth columns of a maxtix using a Moving Average filter
%   f_stdErr           - returns the standard error
%   f_stnd             - standardize values of a matrix, column-wise (= z-scores)
%   f_struct2flat      - turn a structure into a flat table & optionally export
%   f_sub              - subsample data every 6 hrs
%   f_subsetDisPW      - extract subsets of distance matrix based on all pairs of grouping factor
%   f_svd              - singular value decomposition of a matrix
%   f_symb             - utility program for selecting plot symbols
%   f_trajectory       - get distance along a trajectory
%   f_transform        - transform data by 7 methods
%   f_unwrap           - unwrap lower tri-diagonal (w/o diag) of symmetric distance matrix into a column vector
%   f_vecAngle         - counter-clockwise angle between 2 points
%   f_vecDiagram       - plot Progressive Vector Diagram
%   f_vecMagDir        - get magnitude & direction from u,v vector components
%   f_vecPlot          - plot time series of velocity vectors
%   f_vecRot           - rotate vectors (U,V) by angle THETA
%   f_vectorfit        - plot environmental ordination vectors via multiple linear regression
%   f_vecTrans         - transform 2d vector coordinates
%   f_vecTrans3d       - transform 3d vector coordinates
%   f_vecUV            - returns U,V components of a vector given its magnitude & direction
%   f_vonBert          - von Bertalanffy growth function
%   f_vonBertAge       - use VBGF to estimate age of a fish of known length
%   f_vonBertModel     - nonlinear regression estimate of von Bertalanffy growth model
%   f_vonBertModelInit - initial parameter estimate of a von Bertalanffy growth model
%   f_wascores         - weighted-averages scores of species for a site ordination
%   f_windCman         - process CMAN (or NDBC) historical wind data
%   f_windstress       - wind stress in dynes/cm^2
%   f_writeTxt         - write character array to a text file
%   f_wtMean           - weighted mean
%   f_xdiss            - calculate Extended Dissimilarities from a symmetric distance matrix
%   m_2D_surf          - Draws a 2D surface on a map
%   m_extContours      - extract contours from an M_Map map
%   m_ginput           - get lon/lat coordinates from an M_Map figure