MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: NonlinearModelFit on correlated data

  • To: mathgroup at
  • Subject: [mg104911] Re: [mg104861] NonlinearModelFit on correlated data
  • From: Darren Glosemeyer <darreng at>
  • Date: Fri, 13 Nov 2009 05:53:31 -0500 (EST)
  • References: <>

L.J.A.Vasquez at wrote:
> Hello.
> I want to perform a nonlinear model fitting on data points with some
> significant correlations in them.
> For these correlated data points, the definition of chi^2 is
> chi^2=Sum_i  Sum_j  [ y_i - f(y_i) ] * [ y_j - f(y_j) ] / Covariance(y_i,
> y_j) .
> Currently, I am using the Mathematica function "NonlinearModelFit" which
> unfortunately only  assumes statistically independent data points and does
> not allow to incorporate correlation in the data points when estimating
> the best fit parameters.
> I have looked at other data fitting functions such as FindFit and even
> GeneralizedLinearFit but i haven't found one yet that takes into account
> correlated data.
> Maybe i am missing something.  Is there a way to perform a nonlinear fit
> for a correlated data points? Is there a way to change how the chi^2 is
> being defined in any of the Mathematica fitting function? Perhaps, there
> is an additional Mathematica package that is available and deals with this
> problem.  Any help is most welcome.
> With my best regards,

You could try approaching this as a direct optimization problem. One way 
to do that is to construct the chi^2 expression and use FindMinimum or 
NMinimize to minimize the chi^2 to obtain parameter estimates.

Darren Glosemeyer
Wolfram Research

  • Prev by Date: Re: Solving ODE for rotational, irrotational vector field
  • Next by Date: Text on a curve, revisited and modified
  • Previous by thread: NonlinearModelFit on correlated data
  • Next by thread: Re: NonlinearModelFit on correlated data