# Thread: Statistically different tops for 2 parabolas with non-normal distributed data

1. ## Statistically different tops for 2 parabolas with non-normal distributed data

For situation 1 and 2 from my data came two parabola relations between variables X and Y. (the two parabolas have both a maximum)

My question is whether it is possible to test if the two maxima of the parabolas are significantly different from each other, since the tops of the parabolas are unequal to the means of the data. (for both parabolas, the vast majority of the observations is left from the maximum)

If it is possible, how can I measure this? (e.g., which test can do this?)

Thanks a lot

2. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

You could do this a variety of ways. None of which I would describe as 'simple' though. To start with let me as a question - if you fit a quadratic regression (a regression model of the form ) to each set of data do the residuals appear to be approximately normally distributed?

3. ## The Following User Says Thank You to Dason For This Useful Post:

rob1 (08-02-2014)

4. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

Yes, I did the quadratic regressions, and the residuals appear to be normally distributed.

I don't know whether this matters, but I use SPSS for the analysis.

5. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

So I try to continue from Dason's question, with the notation proposed by Dason which is standard.

From the elementary knowledge about the parabola, we know that the maximum is

Note: I am assuming you are not referring to the axis of symmetry , but actually the calculations can be similarly applied

So you will be testing

where representing the parameters of the two situations respectively.

For modelling prospective, you need to specify whether the two situations are independent. If the two variables are the "same" but in two different "situations" then it is likely that the pair in two situations are dependent. Using dummy variable, the model could be written as

where is the dummy variable indicating the two situation. The homogenity of the error can be also addressed as well.

One way to test it is using the bootstrap; another way is using Delta's method. Let's confirm the model set-up first.

6. ## The Following User Says Thank You to BGM For This Useful Post:

rob1 (08-02-2014)

7. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

Thanks so far

I read your part again, and I confirm the model set-up, except that I am comparing the axis of symmetry, instead of the maxima.

8. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

Bootstrapping would probably be the easiest approach but otherwise I would suggest the delta method as bgm mentioned.

9. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

First I must say that I have little experience in implementing the bootstrap procedure in practice, so please correct me if I am wrong.

Now you want to test

which is equivalent to

The basic idea is that as the estimators are consistent estimators, we may use

as the test statistic and reject when it is significantly different from 0.

To determine whether it is significant or not, you need to determine the distribution of under and find out the corresponding quantiles (with the given significance level) and use that to give the rejection/acceptance region.

The steps could be like the following:

1. Suppose you have pairs of data for each situation. Now you re-sample from the original sample with replacement with the same sample size for each situation.

2. Using the generated sample, now estimate all the parameters under the constraint - which you need to jointly estimate for both situation and you may need to use the Lagrange multiplier if you are seeking a closed-form solution.

3. Calculate the test statistic in this sample, and record it.

4. Return to step 1 and repeat for times. Use the sample percentile of these recorded to construct the acceptance region. E.g. if your significance level is then you will use the 2.5 and 97.5 percentile as the pair of end-points for the acceptance region (interval).

Once you obtain the acceptance region, you can calculate for the original sample again and make the decision.

10. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

Actually this would probably be fairly simple to do using non-linear regression... What software are you using?

11. ## Re: Statistically different tops for 2 parabolas with non-normal distributed data

If one set up a model of the form

you could fit that using a non-linear regression routine. If was an indicator that was 0 for group 1 and 1 for group 2 then if you're interested in asking the question "does the max occur at different x-values for these two groups" this translates into testing the null hypothesis .

Really all this is is a way of writing a quadratic function in the form
which is functionally equivalent to the typical way we write quadratics when doing regression but takes a non-linear form in the parameters. We allowed the two groups to have different parameters.

 Tweet

#### Posting Permissions

• You may not post new threads
• You may not post replies
• You may not post attachments
• You may not edit your posts