

IceS = reshapetimeseries(t, extent_S, 'bin', 52) īy default, the function calculates the mean and upper/lower bounds of the data at each x-location in the data. IceN = reshapetimeseries(t, extent_N, 'bin', 52) The early data uses a two-day interval, while the later is daily we'll bin data weekly to avoid any artifacts from this change: load seaice_extent.mat We'll start by reshaping the two timeseries to year x time-of-year matrices using reshapetimeseries, and concatenating these. Example: Sea ice interannual variabilityįor our example, we'll use the sea ice extent timeseries as our ensemble data. By default, the center lines/bars/points are drawn with the indicated color (saturation = 1), and colors of bounds objects (patch, errorbars, etc) use lighter versions of the color, getting lighter as bounds expand (default limits are ).Įnsemble2bnd(., 'cmap', colors) provides an alternative color order for the datasets plotted.Įnsemble2bnd(., 'axis', hax) adds any plots to the indicated axis rather than the current axis.Įnsemble2bnd(., 'whisker', w) changes te whisker length factor used to calculated outliers for the boxplot (boxplot option only). Įnsemble2bnd(., 'alpha', true) indicates that boundedline patches should be plotted using tranparent patches rather than opaque patches.Įnsemble2bnd(., 'tlim', tlim) changes the color limits associated with the patches or errorbars associated with percentiles. These should come in pairs, from lowest to highest percentile. The default is 'xye'.Įnsemble2bnd(., 'center', center) switches the center-line statistic between the 'mean' (default) and 'median' of each dataset.Įnsemble2bnd(., 'prc', prc) specifies the percentile values corresponding to lower and upper bounds. The default, 'none', indicates no plots will be created.Įnsemble2bnd(., 'dims', dims) allows for permutation of the input matrix y via a 3-letter string consisting of the letters x, y, and e that indicate which dimensions in the input matrix correspond to the independent coordinate, the different datasets, and the ensembles per dataset, respectively. = ensemble2bnd(x,y, 'plot', plottype) plots the datasets as either lines with shaded patches indicating percentiles ('boundedline'), lines with errorbars ('errorbar'), unstacked bar plots with errorbars ('bar'), or boxplots ('boxplot'). (See 'dims' parameter to permute these dimensions) length of x input), ny corresponds to the number of different datasets, and nens is the number of ensemble members in each dataset. The x data should be a vector array, and y is a nx x ny x nens array, where nx corresponds to the number of x-dimension points (i.e. = ensemble2bnd(x,y) calculates the mean and upper/lower bounds of y at each x-coordinate.

= ensemble2bnd(., 'whisker', w) Description Example: Sea ice interannual variability.
