[R] How do I test against a simple null that two regressions coefficients are equal?

chen jia chen_1002 at fisher.osu.edu
Thu Jul 8 16:58:01 CEST 2010


OK. Thanks again.

I will read the references more.

Best,
Jia

On Thu, Jul 8, 2010 at 10:51 AM,  <markleeds at verizon.net> wrote:
> hi: no. it's not the same. if you read the paper that I referenced last
> night, that explains how to do the following  test :
>
> Ho: R2 = R1
>
> H1: R2 != R1
>
> that's a different test from what you did but i think it's what you want.
>
>
>
>
> On Jul 8, 2010, chen jia <chen_1002 at fisher.osu.edu> wrote:
>
> Thanks, Chuck. I am reading the references, which are helpful.
>
> Just to understand what I have done wrong here,
>
> I proposed an alternative testing strategy:
> I run regressions y = a3 + b1 * x + b2 * z + e3 and test alternative
> hypothesis b1 != b2 against the null hypothesis b1 = b2 in this
> equation.
>
> Is it this the same test as
>
> y = a1 + b1*x + e1
> y = a2 + b2*x + e2
> test alternative hypothesis b1 != b2 against null hypothesis b1 = b2.
>
> Best,
> Jia
>
> On Wed, Jul 7, 2010 at 11:12 PM, Charles C. Berry <cberry at tajo.ucsd.edu>
> wrote:
>> On Wed, 7 Jul 2010, chen jia wrote:
>>
>>> Hi there,
>>>
>>> I run two regressions:
>>>
>>> y = a1 + b1 * x + e1
>>> y = a2 + b2 * z + e2
>>>
>>> I want to test against the null hypothesis: b1 = b2.  How do I design the
>>> test?
>>>
>>
>> You are testing a non-nested hypothesis, which requires special handling.
>>
>> The classical test is due to Hotelling, but see the references (and R code
>> snippets) in this posting:
>>
>>        http://markmail.org/message/egnowmdzpzjtahy7
>>
>> (it is the merest coincidence that the above thread was initiated by Mark
>> Leeds and that the URL is 'markmail' :-) )
>>
>> HTH,
>>
>> Chuck
>>
>>
>>> I think I can add two equations together and divide both sides by 2:
>>> y = 0.5*(a1+a2) + 0.5*b1 * x + 0.5*b2 * z + e3, where e3 = 0.5*(e1 + e2).
>>> or just y = a3 + 0.5*b1 * x + 0.5*b2 * z + e3
>>>
>>> If I run this new regression, I can test against the null b1 = b2 in
>>> this regression.  Is it an equivalent test as the original one? If
>>> yes, how do I do that in R?
>>>
>>> Alternatively, I think I can just test against the null:
>>> correlation(y, x) = correlation(y, z), where correlation(. , .) is the
>>> correlation between two random variables. Is this equivalent too? If
>>> yes, how do I do it in R?
>>>
>>> Thanks.
>>>
>>> Best,
>>> Jia
>>>
>>> --
>>>                        Ohio State University - Finance
>>>                                  248 Fisher Hall
>>>                                   2100 Neil Ave.
>>>                             Columbus, Ohio  43210
>>>                            Telephone: 614-292-2830
>>>                      http://www.fisher.osu.edu/~chen_1002/
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>> Charles C. Berry                            (858) 534-2098
>>                                            Dept of Family/Preventive
>> Medicine
>> E mailto:cberry at tajo.ucsd.edu               UC San Diego
>> http://famprevmed.ucsd.edu/faculty/cberry/  La Jolla, San Diego 92093-0901
>>
>>
>>
>
>
>
> --
> Ohio State University - Finance
> 248 Fisher Hall
> 2100 Neil Ave.
> Columbus, Ohio 43210
> Telephone: 614-292-2830
> http://www.fisher.osu.edu/~chen_1002/
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



-- 
                         Ohio State University - Finance
                                   248 Fisher Hall
                                    2100 Neil Ave.
                              Columbus, Ohio  43210
                             Telephone: 614-292-2830
                       http://www.fisher.osu.edu/~chen_1002/



More information about the R-help mailing list