Built-in Mathematica Symbol

FindFit

FindFit[data, expr, pars, vars]
finds numerical values of the parameters pars that make expr give a best fit to data as a function of vars. The data can have the form {{x₁, y₁, ..., f₁}, {x₂, y₂, ..., f₂}, ...}, where the number of coordinates x, y, ... is equal to the number of variables in the list vars. The data can also be of the form {f₁, f₂, ...}, with a single coordinate assumed to take values 1, 2, ....

FindFit[data, {expr, cons}, pars, vars]
finds a best fit subject to the parameter constraints cons.

MORE INFORMATION

FindFit returns a list of replacements for par₁, par₂, ....

The expression expr must yield a numerical value when pars and vars are all numerical.

The expression expr can depend either linearly or nonlinearly on the par_i.

In the linear case, FindFit finds a globally optimal fit.

In the nonlinear case, it finds in general only a locally optimal fit.

FindFit[data, expr, {{par₁, p₁}, {par₂, p₂}, ...}, vars] starts the search for a fit with {par₁->p₁, par₂->p₂, ...}.

FindFit by default finds a least-squares fit.

The option NormFunction->f specifies that the norm f[residual] should be minimized.

The constraints cons can contain equations, inequalities or logical combinations of these.

The following options can be given:

AccuracyGoal	Automatic	the accuracy sought
EvaluationMonitor	None	expression to evaluate whenever expr is evaluated
Gradient	Automatic	the list of gradient components for expr
MaxIterations	Automatic	maximum number of iterations to use
Method	Automatic	method to use
NormFunction	Norm	the norm to minimize
PrecisionGoal	Automatic	the precision sought
StepMonitor	None	expression to evaluate whenever a step is taken
WorkingPrecision	Automatic	the precision used in internal computations

The default settings for AccuracyGoal and PrecisionGoal are WorkingPrecision/2.

The settings for AccuracyGoal and PrecisionGoal specify the number of digits to seek in both the values of the parameters returned, and the value of the NormFunction.

FindFit continues until either of the goals specified by AccuracyGoal or PrecisionGoal is achieved.

Possible settings for Method include "ConjugateGradient", "Gradient", "LevenbergMarquardt", "Newton", "NMinimize" and "QuasiNewton", with the default being Automatic.

EXAMPLES

CLOSE ALL

Basic Examples (1)

Find a nonlinear fit to a list of primes:

Fit the best-fit parameters

,

, and

:

Evaluate the fitted function:

Find a nonlinear fit to a list of primes:

In[1]:=

Out[1]=

Fit the best-fit parameters

,

, and

:

In[2]:=

Out[2]=

Evaluate the fitted function:

In[3]:=

Out[3]=

Scope (7)

Fit data to a model of exponential decay:

The resulting model function:

Show the data with the model:

Separate the time from the response:

Looking at residuals gives a good indication whether the model is a good fit:

Use a linear fit on the logarithm of the data for a model of exponential decay:

The logarithm of the exponential model is a linear model for the logarithm of the data:

The exponential model formed by exponentiating the model for the logarithm:

Exponential fit as parameters in a differential equation:

The model function is the solution of a differential equation:

Using caching is faster:

Find a constrained fit:

Fit to the model with positive amplitude and frequency between 1 and 2:

Compare the fitted model to the data:

The residuals show a pattern, indicating that the frequency constraint is too strict:

Give starting values for some parameters to get a better fit:

Search with all parameters starting at 1:

Search with a better starting value for the parameter c:

Compare the quality of the two fits:

Fit a surface to data in two-dimensions:

Find the fit starting from an approximate position for the peak:

Show the fitted surface with the data:

Fit a model to data in four dimensions:

The residual is comparable in size to the randomness added to the data:

Options (7)

Count evaluations of the model with numerical values of the parameters:

Specify the model gradient to avoid problems with a removable singularity:

With symbolic derivatives FindFit fails since the derivative for Sinc is given as a generic formula:

The gradient has a singularity when

since there is a point in the data with

:

Specify finite differences to avoid the removable singularity:

Show the data with the fit:

When convergence is slow, increasing MaxIterations may allow convergence to the best fit:

With the default, convergence is not reached:

Convergence is reached before 1000 iterations are used:

Use Method->NMinimize to search more exhaustively for a global minimum:

Find a local minimum:

Search for a global minimum:

Compare the quality of the two fits (the global best fit is shown in red):

The default is to find the best least-squares fit:

Use the

-norm instead:

Show steps taken in parameter space to fit the Michaelis-Menten model to experimental data:

Data from an experiment on the rate of an enzymatic reaction:

Show how the parameters evolve during the search:

Use a higher WorkingPrecision to get a more precise parameter estimate:

Applications (7)

Fit a list of Fibonacci numbers:

The exponent is approximately the golden ratio:

Generate data with some noise present:

Fit the data:

Determine the approximate complexity of a computation:

Fit to the theoretical model for the complexity of an FFT:

Show the timing data and the fit in a log-log plot:

Fit a model for lubricant viscosity to experimental data:

Model for lubricant viscosity as a function of temperature

and pressure/1000

:

Scale the pressure to improve the conditioning of the model:

Find the parameters

that best fit the scaled data:

Show a plot of the fit with the scaled experimental data:

Fit a model for the rate of isomerization of n-pentane to isopentane with experimental data:

Nonlinear model for rate as a function of hydrogen

, n-pentane

, and isopentane

:

Find the parameters

that best fit the data starting with values having the right size:

Find parameters in a differential equation model that gives the best fit:

Define a model so the solution for a set of parameters

is reused for all

values:

Find the fit with tolerance appropriate to the precision of the solutions of the ODE:

Show the fitted model with the data and the residuals:

Find how a nonsoliton initial condition for the KdV equation breaks down into two solitons:

The initial wave breaks into two separate waves:

Get the spatial grid used by NDSolve to solve the PDE:

The discretized solution at the final time:

Model of two solitons with different velocities:

Find fit with starting conditions estimated from the graph of the discretized solution:

Plots of the two-soliton model and residuals at intermediate times:

Properties & Relations (4)

FindFit gives the parameter estimates:

NonlinearModelFit allows for extraction of additional information about the fitting:

Extract the parameter estimates:

Extract additional results and diagnostics:

For linear models, Fit and FindFit are equivalent and use LeastSquares:

Use Fit with the basis {1, x, x²}:

This is equivalent to FindFit with the model function a linear combination of the basis:

Construct a design matrix and use LeastSquares to find the coefficients:

LinearModelFit fits linear models and gives additional results about the fitting:

FindFit effectively uses FindMinimum to minimize a norm of the residual function:

The residual is the difference between response and model prediction:

FindFit gives a result consistent with normally distributed responses by default:

LogitModelFit assumes binomially distributed responses:

The fits are not identical:

The same is true for ProbitModelFit:

GeneralizedLinearModelFit allows for binomial and other distributional structures:

Model gamma-distributed responses:

Compare with the result for normally distributed responses:

Possible Issues (1)

Good starting values may be necessary to get a reasonable fit:

Find the fit with the default starting values of 1. The local minimum is far from correct:

Typically, better starting values are required to get a good fit: