Nonlinear regression
Part of a series on 
Regression analysis 

Models 
Estimation 
Background 

In statistics, nonlinear regression is a form of regression analysis in which observational data are modeled by a function which is a nonlinear combination of the model parameters and depends on one or more independent variables. The data are fitted by a method of successive approximations.
General
In nonlinear regression, a statistical model of the form,
 [math]\displaystyle{ \mathbf{y} \sim f(\mathbf{x}, \boldsymbol\beta) }[/math]
relates a vector of independent variables, [math]\displaystyle{ \mathbf{x} }[/math], and its associated observed dependent variables, [math]\displaystyle{ \mathbf{y} }[/math]. The function [math]\displaystyle{ f }[/math] is nonlinear in the components of the vector of parameters [math]\displaystyle{ \beta }[/math], but otherwise arbitrary. For example, the Michaelis–Menten model for enzyme kinetics has two parameters and one independent variable, related by [math]\displaystyle{ f }[/math] by:^{[loweralpha 1]}
 [math]\displaystyle{ f(x,\boldsymbol\beta)= \frac{\beta_1 x}{\beta_2 + x} }[/math]
This function is nonlinear because it cannot be expressed as a linear combination of the two [math]\displaystyle{ \beta }[/math]s.
Systematic error may be present in the independent variables but its treatment is outside the scope of regression analysis. If the independent variables are not errorfree, this is an errorsinvariables model, also outside this scope.
Other examples of nonlinear functions include exponential functions, logarithmic functions, trigonometric functions, power functions, Gaussian function, and Lorentz distributions. Some functions, such as the exponential or logarithmic functions, can be transformed so that they are linear. When so transformed, standard linear regression can be performed but must be applied with caution. See Linearization§Transformation, below, for more details.
In general, there is no closedform expression for the bestfitting parameters, as there is in linear regression. Usually numerical optimization algorithms are applied to determine the bestfitting parameters. Again in contrast to linear regression, there may be many local minima of the function to be optimized and even the global minimum may produce a biased estimate. In practice, estimated values of the parameters are used, in conjunction with the optimization algorithm, to attempt to find the global minimum of a sum of squares.
For details concerning nonlinear data modeling see least squares and nonlinear least squares.
Regression statistics
The assumption underlying this procedure is that the model can be approximated by a linear function, namely a firstorder Taylor series:
 [math]\displaystyle{ f(x_i,\boldsymbol\beta) \approx f(x_i,0) + \sum_j J_{ij} \beta_j }[/math]
where [math]\displaystyle{ J_{ij} = \frac{\partial f(x_i,\boldsymbol\beta)}{\partial \beta_j} }[/math]. It follows from this that the least squares estimators are given by
 [math]\displaystyle{ \hat{\boldsymbol{\beta}} \approx \mathbf { (J^TJ)^{1}J^Ty}, }[/math]
compare generalized least squares with covariance matrix proportional to the unit matrix. The nonlinear regression statistics are computed and used as in linear regression statistics, but using J in place of X in the formulas.
When the function [math]\displaystyle{ f(x_i,\boldsymbol\beta) }[/math] itself is not known analytically, but needs to be linearly approximated from [math]\displaystyle{ n+1 }[/math], or more, known values (where [math]\displaystyle{ n }[/math] is the number of estimators), the best estimator is obtained directly from the Linear Template Fit as ^{[1]}[math]\displaystyle{ \hat{\boldsymbol\beta} = ((\mathbf{Y\tilde{M}})^\mathsf{T} \boldsymbol\Omega^{1} \mathbf{Y\tilde{M}})^{1}(\mathbf{Y\tilde{M}})^\mathsf{T}\boldsymbol\Omega^{1}(\mathbf{d}\mathbf{Y\bar{m})} }[/math] (see also linear least squares).
The linear approximation introduces bias into the statistics. Therefore, more caution than usual is required in interpreting statistics derived from a nonlinear model.
Ordinary and weighted least squares
The bestfit curve is often assumed to be that which minimizes the sum of squared residuals. This is the ordinary least squares (OLS) approach. However, in cases where the dependent variable does not have constant variance, a sum of weighted squared residuals may be minimized; see weighted least squares. Each weight should ideally be equal to the reciprocal of the variance of the observation, but weights may be recomputed on each iteration, in an iteratively weighted least squares algorithm.
Linearization
Transformation
Some nonlinear regression problems can be moved to a linear domain by a suitable transformation of the model formulation.
For example, consider the nonlinear regression problem
 [math]\displaystyle{ y = a e^{b x}U \,\! }[/math]
with parameters a and b and with multiplicative error term U. If we take the logarithm of both sides, this becomes
 [math]\displaystyle{ \ln{(y)} = \ln{(a)} + b x + u, \,\! }[/math]
where u = ln(U), suggesting estimation of the unknown parameters by a linear regression of ln(y) on x, a computation that does not require iterative optimization. However, use of a nonlinear transformation requires caution. The influences of the data values will change, as will the error structure of the model and the interpretation of any inferential results. These may not be desired effects. On the other hand, depending on what the largest source of error is, a nonlinear transformation may distribute the errors in a Gaussian fashion, so the choice to perform a nonlinear transformation must be informed by modeling considerations.
For Michaelis–Menten kinetics, the linear Lineweaver–Burk plot
 [math]\displaystyle{ \frac{1}{v} = \frac{1}{V_\max} + \frac{K_m}{V_{\max}[S]} }[/math]
of 1/v against 1/[S] has been much used. However, since it is very sensitive to data error and is strongly biased toward fitting the data in a particular range of the independent variable, [S], its use is strongly discouraged.
For error distributions that belong to the exponential family, a link function may be used to transform the parameters under the Generalized linear model framework.
Segmentation
The independent or explanatory variable (say X) can be split up into classes or segments and linear regression can be performed per segment. Segmented regression with confidence analysis may yield the result that the dependent or response variable (say Y) behaves differently in the various segments.^{[2]}
The figure shows that the soil salinity (X) initially exerts no influence on the crop yield (Y) of mustard, until a critical or threshold value (breakpoint), after which the yield is affected negatively.^{[3]}
See also
 Nonlinear least squares
 Curve fitting
 Generalized linear model
 Local regression
 Response modeling methodology
 Genetic Programming
 Multi expression programming
 Linear or quadratic template fit
References
 ↑ Britzger, Daniel (2022). "The Linear Template Fit". Eur. Phys. J. C 82: 731. doi:10.1140/epjc/s1005202210581w.
 ↑ R.J.Oosterbaan, 1994, Frequency and Regression Analysis. In: H.P.Ritzema (ed.), Drainage Principles and Applications, Publ. 16, pp. 175224, International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands. ISBN:9070754339 . Download as PDF : [1]
 ↑ R.J.Oosterbaan, 2002. Drainage research in farmers' fields: analysis of data. Part of project “Liquid Gold” of the International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands. Download as PDF : [2]. The figure was made with the SegReg program, which can be downloaded freely from [3]
Notes
 ↑ This model can also be expressed in the conventional biological notation:
 [math]\displaystyle{ v = \frac{V_\max\ [\mbox{S}]}{K_m + [\mbox{S}]} }[/math]
Further reading
 Bethea, R. M.; Duran, B. S.; Boullion, T. L. (1985). Statistical Methods for Engineers and Scientists. New York: Marcel Dekker. ISBN 082477227X.
 Meade, N.; Islam, T. (1995). "Prediction Intervals for Growth Curve Forecasts". Journal of Forecasting 14 (5): 413–430. doi:10.1002/for.3980140502.
 Schittkowski, K. (2002). Data Fitting in Dynamical Systems. Boston: Kluwer. ISBN 1402010796.
 Seber, G. A. F.; Wild, C. J. (1989). Nonlinear Regression. New York: John Wiley and Sons. ISBN 0471617601.
Original source: https://en.wikipedia.org/wiki/Nonlinear regression.
Read more 