Med2top

From Eigenvector Research Documentation Wiki
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

Purpose

Fits a constant to top/(bottom) of data.

Synopsis

[yf,residual,options] = med2top(y,options)

Description

MED2TOP is similar to LSQ2TOP with a 0 order polynomial, it can be considered an asymmetric estimate of the mean.

For fitting to the bottom:

>> tsq = residual/res; % (res) is an input

>> tsqst = ttestp(1-options.tsqlim,5000,2); % T-test limit from table

>> ii = find(tsq>-tsqst); % finds samples below the line

The ii samples are kept for the next estimate of (yf):

>> yf = median(y(ii));

Inputs

  • y = trace to be filtered, Mx1 vector.

Outputs

  • yf = scalar, estimate of filtered data.
  • residual = y - yf.
  • options = input options echoed back, the field initwt may have been modified.

Options

  • ''' options = a structure array with the following fields.
  • display: [ {'off'} | 'on'] Governs screen display to command line.
  • trbflag: [ {'top'} | 'bottom' | 'middle'] flag that tells algorithm to fit to the top, bottom, or middle of the data cloud.
  • tsqlim: [ 0.99 ] limit that govers whether a data point is outside the fit residual defined by input (res).
  • initwt: [ ] empty or Mx1 vector of initial weights (0<=w<=1).

See Also

baseline, baslinew, fastnnls, lsq2top