Re: NonlinearModelFit and "ANOVATable" and "ParameterConfidenceIntervalTable"
- To: mathgroup at smc.vnet.net
- Subject: [mg121311] Re: NonlinearModelFit and "ANOVATable" and "ParameterConfidenceIntervalTable"
- From: Darren Glosemeyer <darreng at wolfram.com>
- Date: Fri, 9 Sep 2011 05:50:19 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <201109070939.FAA21479@smc.vnet.net>
On 9/7/2011 4:39 AM, HansChristopf wrote:
> Dear friends,
>
> I'm a Mathematica beginner and use in the course of a small project at my college in Germany. I'm using NonlinearModelFit in order to do fit a model to experimental data. I post the code in the following:
>
> ----------------------------------------------------------------------------------
>
> In[1]:= data = {{1., 0.75}, {2., 0.89}, {3., 0.42}, {4., 0.99}, {5.,0.84}, {6., 0.34}, {7., 0.83}, {8., 0.93}, {9., 0.76}, {10.,0.11}};
> nlm = NonlinearModelFit[data, {Exp[a + b x^2], b< -1/2}, {{a, 0}, {b, -1}}, x]
>
> In[2]:= nlm[{"ANOVATable", "ParameterConfidenceIntervalTable"}]
>
> During evaluation of In[9]:= FittedModel::constr: The property values {ANOVATable,ParameterConfidenceIntervalTable} assume an unconstrained model. The results for these properties may not be valid, particularly if the fitted parameters are near a constraint boundary.>>
>
> ----------------------------------------------------------------------------------
>
> Now my questions to my problem:
>
> I cannot interpret these results ?
>
> 1. Are the results from the ANOVATable the same as the results from the ANOVA (Analysis of Variance) ?
> 2. Which formulas uses Mathematica to obtain these results ?
> 3. How Mathematica calculates the ANOVATable EXACTLY ???????
>
>
> Greatings, Hans
>
For question 1, if you're asking about the Analysis of Variance model
(linear model with categorical predictors), the answer is that they are
conceptually similar but the models are different. The ANOVATable of
nonlinear models decomposes the sources of variation in the data as does
the ANOVA table for an Analysis of Variance model.
Here are the formulas for the terms in the SS (sum of squares) column:
In[10]:= error = Total[nlm["FitResiduals"]^2]
Out[10]= 4.60397
In[11]:= uncorrected = Total[nlm["Response"]^2]
Out[11]= 5.4758
In[12]:= corrected = Total[(nlm["Response"] - Mean[nlm["Response"]])^2]
Out[12]= 0.76984
In[13]:= model = uncorrected - error
Out[13]= 0.871827
For the DF (degrees of freedom) terms, the model DF is the number of
parameters in the model, uncorrected DF is the number of data points,
and the error DF is their difference. The corrected DF is reduced by 1
because of the correction by the mean estimate. The MS (mean square)
values are obtained by dividing the SS terms by the associated DF values.
Darren Glosemeyer
Wolfram Research
- References:
- NonlinearModelFit and "ANOVATable" and "ParameterConfidenceIntervalTable"
- From: HansChristopf <hanschristopf@yahoo.de>
- NonlinearModelFit and "ANOVATable" and "ParameterConfidenceIntervalTable"