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
I don't have emotions and sometimes that makes me very sad.
rob1 (08-02-2014)
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.
Last edited by rob1; 07-30-2014 at 06:14 PM.
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.
rob1 (08-02-2014)
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.
Last edited by rob1; 08-01-2014 at 03:21 PM.
Bootstrapping would probably be the easiest approach but otherwise I would suggest the delta method as bgm mentioned.
I don't have emotions and sometimes that makes me very sad.
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.
Actually this would probably be fairly simple to do using non-linear regression... What software are you using?
I don't have emotions and sometimes that makes me very sad.
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.
I don't have emotions and sometimes that makes me very sad.
Tweet |