Mailing-List: contact user-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: "Commons Users List" <user@commons.apache.org>
Received-SPF: pass (athena.apache.org: domain of ted.dunning@gmail.com
 designates 209.85.213.175 as permitted sender)
MIME-Version: 1.0
In-Reply-To: 
 <CAKUdAB67f7Ng2pOwu8zuyPDGWLLxhU-pyoTjuPVEyx9hbaSLRg@mail.gmail.com>
References: 
 <CAKUdAB67f7Ng2pOwu8zuyPDGWLLxhU-pyoTjuPVEyx9hbaSLRg@mail.gmail.com>
From: Ted Dunning <ted.dunning@gmail.com>
Date: Wed, 13 Aug 2014 22:22:24 -0700
Message-ID: 
 <CAJwFCa1+p0N42yxaOwc2C_OLfQcg8afaFq+NHFOw6pPmjFNoOw@mail.gmail.com>
Subject: Re: [math] Curve fitting ...
To: Commons Users List <user@commons.apache.org>
Content-Type: multipart/alternative; boundary=001a11c30ad0c6f8050500901842

--001a11c30ad0c6f8050500901842
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: quoted-printable

Have you considered using an interactive system like R, Matlab or Octave?

You might be happier.

Or even have you considered goal search in Excel?


On Wed, Aug 13, 2014 at 6:08 PM, South Light <southlight25@gmail.com> wrote=
:

> Hi,
>
> May be someone can help me with this problem.
>
> Given the follow function: y =3D 10 ^ ((x + 82) / (-10 * A))
>
> I would like to found the A value witch curve fit better for a set of x,y
> values, usually the set is about 20 to 25 x,y values.
>
> I use the CurveFitter class and the ParametricUnivariateFunction
>
>
> ParametricUnivariateFunction function =3D new  ParametricUnivariateFuncti=
on()
> {
>
> =E2=80=8B  =E2=80=8B
> @Override
> =E2=80=8B  =E2=80=8B
> public double[] gradient(double x, double[] params) {
> =E2=80=8B
> (????? comment)
> =E2=80=8B  =E2=80=8B
> }
>
> =E2=80=8B  =E2=80=8B
> @Override
> =E2=80=8B  =E2=80=8B
> public double value(double x, double[] params) {
>
> =E2=80=8B    =E2=80=8B
> double a =3D params[0];
>
> =E2=80=8B    =E2=80=8B
> return Math.pow(10, ((x + 82) /
> =E2=80=8B(=E2=80=8B
> -10 *
> =E2=80=8Ba
> =E2=80=8B)=E2=80=8B
> ));
>
> =E2=80=8B  =E2=80=8B
> }
>
> };
>
> LevenbergMarquardtOptimizer optimizer =3D new LevenbergMarquardtOptimizer=
();
>
> CurveFitter<ParametricUnivariateFunction> fitter =3D new
> CurveFitter<ParametricUnivariateFunction>(optimizer);
>
> double[] x =3D {
> -82
> ,
> =E2=80=8B-85
> ,
> =E2=80=8B-89
>  };
>
> double[] y =3D {
> =E2=80=8B1
> ,
> =E2=80=8B1.4
> ,
> =E2=80=8B2
>  };
>
> for (int i =3D 0; i < x.length; i++)
> =E2=80=8B  =E2=80=8B
> fitter.addObservedPoint(x[i], y[i]);
>
> double[] result =3D fitter.fit(function, new double[] { 1, 10 });
>
>
>
> =E2=80=8BA. =E2=80=8BIs this the best way to solve the problem or there's=
 another better
> way?
>
> B. What do we need to write on the gradient area (????? comment) ?=E2=80=
=8B
>
> Any help will be more then welcome.
>
> Many thanks !!
>

--001a11c30ad0c6f8050500901842--