fit_sigm - tpcclib 0.8.0 © 2023 by Turku PET Centre
Non-linear fitting of the Hill function (sigmoid function) to a (xi,yi)
curve data. Optionally, a background term B is added to the function:
f(x) = (A * x^n ) / ( x^n + K ) + B , where A>=0, K>0, n>0, B>=0
Function for estimating EC50 from sigmoidal dose-response curve is obtained
by rearranging and marking n=HillSlope, A=Top-Bottom, B=Bottom, and
K=EC50^n :
f(x) = Bottom + (Top-Bottom)/(1 + (EC50/x)^HillSlope)
With HillSlope constrained to 1 (see options below) the function reduces to
unconstrained one site binding hyperbola, f(x)=Top*x/(EC50+x).
Usage: fit_sigm [Options] tacfile [fitfile]
Options:
-EC50
EC50 is estimated by fitting sigmoidal dose-response curve function
when drug concentrations are on a linear scale (not on a log scale);
By default, parameter Top is fitted and Bottom is constrained to 0,
which can be changed with options (see below).
-w1 All weights are set to 1.0 (no weighting); by default, weights in
data file are used, if available.
-wf
Weight by sampling interval.
-n=<value>
Parameter n or HillSlope is constrained to the given value; set -n=1
to fit one site binding hyperbola.
-B Parameter B or Bottom is fitted; by default B=0
-B=<value>
Parameter B or Bottom is constrained to the given value.
-A=<value>
Parameter Top is constrained to the given value; currently this option
does not affect parameter A in the default function.
-MRL
Error is returned if MRL check is not passed.
-res=<Filename>
Fitted parameters are also written in result file format (3).
-svg=<Filename>
Fitted and measured TACs are plotted in specified SVG file.
-fit=<Filename>
Fitted TACs are written in DFT format.
-h, --help
Display usage information on standard output and exit.
-v, --version
Display version and compile information on standard output and exit.
-d[n], --debug[=n], --verbose[=n]
Set the level (n) of debugging messages and listings.
-q, --quiet
Suppress displaying normal results on standard output.
-s, --silent
Suppress displaying anything except errors.
TAC file must contain at least 2 columns, sample times (or concentrations)
and measurement values (or receptor occupancies).
Weights can be specified as usual if data is in DFT format (1).
Program writes the fit start and end times, nr of points, WSS,
and parameters of the fitted function to the fit file (2).
References:
1. http://www.turkupetcentre.net/analysis/doc/format_dft.html
2. http://www.turkupetcentre.net/analysis/doc/format_fit.html
3. http://www.turkupetcentre.net/analysis/doc/format_res.html
See also: fit_hiad, fit_ppf, fit_exp, fit_ratf, fit_bpr, fit2dat, tacinv
Keywords: curve fitting, Hill function, sigmoid, EC50, dose-response curve