Re: Automate datafitting to a series of parameterized function
- To: mathgroup at smc.vnet.net
- Subject: [mg70306] Re: [mg70277] Automate datafitting to a series of parameterized function
- From: "Peng Yu" <pengyu.ut at gmail.com>
- Date: Thu, 12 Oct 2006 05:37:37 -0400 (EDT)
- References: <200610110553.BAA19424@smc.vnet.net> <452D0E8E.8000102@wolfram.com>
On 10/11/06, Darren Glosemeyer <darreng at wolfram.com> wrote:
> The model is a linear model in the basis functions 1, x, x^2,..., so you
> can use Regress to easily get this information.
>
> This will give the fitted function and residuals for the model with n terms.
>
> << Statistics`LinearRegression`
>
> fittedAndResiduals[n_] := {BestFit, FitResiduals} /. Regress[data,
> x^Range[0, n], x,
> RegressionReport -> {BestFit, FitResiduals} ];
>
>
> The mentioned results for a constant model can then be computed directly
> as follows, and for higher order models by choosing the appropriate
> value of n.
>
> res0 = fittedAndResiduals[0];
> error0 = res0[[2]]^2;
> avererror0 = Total[error0]/Length[data]
> maxerror = Max[error0]
>
>
> If there is a maximum power for the terms, say 5, that you wish to
> allow, results for all those models could be obtained as follows.
>
> Table[Block[{res = fittedAndResiduals[n], error},
> error = res[[2]]^2;
> {res[[1]], Total[error]/Length[data], Max[error]}], {n, 0, 5}]
>
>
> Other diagnostics could also be obtained from Regress or computed from
> results obtained from Regress if desired.
Thank you for your reply. I gave an linear fitting example because I
just want to describe the problem simply. But that doesn't mean my
real problem can be reduced to a linear fitting problem.
Is there any way to solve the nonlinear fit problem?
Thanks,
Peng
- References:
- Automate datafitting to a series of parameterized function
- From: "Peng Yu" <pengyu.ut@gmail.com>
- Automate datafitting to a series of parameterized function