1. Home
  2. ALIGNCURVE procedure

ALIGNCURVE procedure

Forms an optimal warping to align an observed series of observations with a standard series (D.B. Baird).


PRINT = string tokens What to print (criterion, ss, warps); default * i.e. nothing
PLOT = string tokens What to plot (series, warping); default * i.e. no plots
WARPPENALTY = scalar The relative penalty to add to the criterion when jumping a unit in one series but not the other; default 1
MAXSTEP = scalar The largest jump that can be made between the two series at a single point; default 1
MAXDIFFERENCE = scalar Sets a limit on size of difference between the series to be squared and added to the criterion (differences greater than this are truncated to MAXDIFFERENCE, thus allowing the effects of outliers to be down-weighted); default * i.e. no limit
USEMEANS = string token Whether to use the means of points covered in one step, rather than the final value, when calculating the sums of squares between the two series (yes, no); default no
FORCEENDALIGNMENT = string token Whether to force the ends of the two series to align, so that warping happens only in the middle of the series (yes, no); default no
WINDOW = scalar Window number for the plots; default 1
KEYWINDOW = scalar Window for the key (zero for no key); default 2


Y = variates Series to be aligned with the standard series
STANDARD = variates Standard series for each Y
WEIGHTS = variates Weights for the contribution of each point to the criterion; default * no weighting
UWARP = variates The warped positions of the unit numbers, required to align Y with STANDARD
YWARP = variates The warped series for Y, i.e. the optimally aligned y-values
CRITERIONVALUE = scalars The criterion value (as optimized during the alignment)
TITLE = text Title for the plots


ALIGNCURVE is useful when you have a series of observations of a theoretical curve in which the x-axis may have been distorted by stretching or compression. The aim is to “warp” the units of the curve so that it optimally matches a standard series. (This can be used, for example, to align electrophoresis gels.) The observed and standard series are specified, in variates, by the Y and STANDARD parameters respectively.

The warping takes the form of one of the series jumping values, so it that makes several steps for a single step of the other series. There is no limitation on the number of times that this can happen along the series, but you can set option FORCEENDALIGNMENT=yes, to require the two series to match at their final values.

The dynamic warping algorithm that is used, seeks to optimize a criterion that is based on the sum of squares of residuals between the warped observed series and the standard series at each point, plus a penalty that depends on the amount of warping.

The WEIGHTS parameter can supply a variate of weights to use for the contribution of each point to the criterion. If the standard curve has non-constant variance around it, typically you would set the weights to be the reciprocal of the variance or equivalently, one over the square of the standard deviation. Suppose, for example, if the curve is the mean of a Poisson process. the mean is equal to the variance for a Poisson distribution, so the you could use weights of 1/STANDARD to allow for the non-constant variance. Alternatively, if the points on the curves were means of differing numbers of replicated observations, you could weight by the recriprocals of the replications. By default there is no weighting.

The MAXDIFFERENCE option sets a limit on contribution of each residual: differences greater than this are truncated (to MAXDIFFERENCE), thus allowing the effects of outliers to be down-weighted. By default there is no limit.

The WARPPENALTY option defines the relative penalty to be added to the criterion when the series are warped by one step in either direction; default 1. The penalty is given by WARPPENALTY multiplied by a scale factor equal to the average mean square of the initial (unwarped) differences between Y and STANDARD. (It is thus scale independent, so multiplying both Y and STANDARD by a constant will still give the same solution.) Setting WARPPENALTY to zero finds the optimally aligned series with no constraint on warping. Conversely, setting WARPPENALTY to a large value discourages warps, as a jump will be taken only if the decrease in the sum of squares is greater than the penalty.

The MAXSTEP option sets a limit on the number of units that can be jumped at any one point (default 1).

By default, the y-value after a jump is taken to be the value at the end of the jump, but you can set option USEMEANS=yes to use the mean of the values at the end and intervening (jumped) points.

The UWARP parameter can save a variate containing the warped units, i.e. the unit of Y that corresponds (after warping) to each unit of STANDARD. The YWARP parameter can save the y-values of the optimally aligned series, and the CRITERIONVALUE parameter can save the optimal criterion value.

The PRINT option controls printed output, with settings:

    criterion to print the smoothing weight, maximum step size and optimal criterion value,
    warps to print the location of the warps,
    ss to print sums of squares before and after alignment.

By default nothing is printed.

The PLOT option controls the plots that are displayed, with settings:

    series plots the aligned and standard series,
    warping plots the steps in the warping that have been used to align the two series.

By default nothing is plotted. The WINDOW option specifies the window to use for the plots (default 1), and the KEYWINDOW option specifies the window for their keys (default 2). You can supply a title for the plots using the TITLE parameter.



Action with RESTRICT

Any restrictions are ignored.

See also


Commands for: Calculations and manipulation.


CAPTION    'ALIGNCURVE examples'; STYLE=meta
CALCULATE  X = !(1...1000)/80
&          S = SIN(X)
&          Y = SIN(1.2*X)        " Y has an increased rate from S " 
ALIGNCURVE [PRINT=criterion,ss; PLOT=series,warping; WARPPENALTY=0;\
           UWARP=WarpedunitsY; YWARP=WWarpedY; CRITERIONVALUE=CY
CALCULATE  Z = SIN(X + SIN(X/2)) " Z has an varying rate from S "
ALIGNCURVE [PRINT=criterion,ss; PLOT=series,warping; WARPPENALTY=0]\
" penalize warps "
ALIGNCURVE [PRINT=criterion,ss,warps; PLOT=series,warping; WARPPENALTY=1]\
           Z; STANDARD=S; UWARP=WarpedunitsZp; YWARP=WarpedZp;\
           CRITERIONVALUE=CZp; TITLE='Z aligned to S, penalty=1'
Updated on June 20, 2019

Was this article helpful?