PeakModel¶

PeakModel()¶

PeakModel Parametric model of a single spectral peak

obj = PeakModel(modelFunClass, paramLinkFuns, paramLinkFunParams, modelOptions, peakName)

Input:

modelFunClass: string, model function class: ‘gaussian’, ‘box’,: ‘harmonic’, ‘skewgamma’, ‘expdecay’, ‘expdoubledecay’, ‘smoothdoubledecay’, ‘linear’
paramLinkFuns: 1 x numPeakParams cell array of strings of link: function classes: ‘sigmoid’ (default), ‘exp’, ‘identity’, ‘poly’
paramLinkFunParams: 1 x numPeakParams cell array of: strings of link function parameters

modelOptions: model specific options peakName: string, descriptive name of peak

PEAK FUNCTIONS:

Gaussian (‘gaussian’)

Function:
y(x) = A * exp(-(w-F)^2 / (2*B))

modelParams:
F - central frequency A - maximum amplitude B - bandwidth (variance)

Box Function (‘box’);

Function:
y(x) = A * exp(-(w-F)^p / (2*B))

modelParams:
F - central frequency A - maximum amplitude B - bandwidth (variance)

modelOptions:
p - polynomial order

Gaussian with Harmonics (‘harmonic’);

Function:
y(x) = SUM h = 0 to H_peak: (Beta^h) * A * exp(-(w - (h+1)*F)^2 / (2*B));

modelParams:
F - central frequency A - maximum amplitude B - bandwidth (variance) Beta - harmonic ampiltude multiplier

modelOptions:
H_peak - number of harmonics

Gamma with Skewness (‘skewgamma’);

Function:
A gamma function with fixed mode, variance, and skew

modelParams:
F - central frequency (mode) A - maximum amplitude B - bandwidth (variance) S - skewness

Exponential Decay (‘expdecay’);

Function:
y(x) = BA.*(1-BR).^x +BO;

modelParams:
BA - maximum amplitude BR - decay rate BO - offset/baseline

Double Exponential Decay (‘expdecay’);

Function:
y(x) = BA1.*(1-BR1).^x + BA2.*(1-BR2).^x +BO;

modelParams:
BA1 - maximum amplitude of first exponential BR1 - decay rate of first exponential BA2 - maximum amplitude of second exponential BR2 - decay rate of second exponential BO - offset/baseline

Piecewise Exponential with Smooth Transition (‘smoothdoubledecay’);

Function:
see: Zelterman et. al, “Piecewise exponential survival curves with smooth transitions”, Mathematical Biosciences, 1994

modelParams:
L1mL2 - difference between the rates of the two exponentials L2 - The rate of the second exponential K - maximum amplitude O - offset/baseline

Last modified 5/24/2019 by Michael Prerau