Hello Statalist!
I am having trouble retrieving the Pr>chi2 for a Test of homogeneity of ORs using mhodds. To clearly explain my situation with completely arbitrary data from online:
. sysuse nlsw88
. mhodds never_married collgrad, by(south)
Maximum likelihood estimate of the odds ratio
Comparing collgrad==1 vs. collgrad==0
by south
-------------------------------------------------------------------------------
----south | Odds Ratio---------chi2(1)--------P>chi2---------[95% Conf. Interval]
----------+--------------------------------------------------------------------
--------0-|---2.076923----------15.94---------0.0001---------1.43906----2.99751
--------1-|---1.483171-----------2.63---------0.1050---------0.91809----2.39607
-------------------------------------------------------------------------------
----Mantel-Haenszel estimate controlling for south
--------------------------------------------------------------------
-----Odds Ratio----chi2(1)-------P>chi2---------[95% Conf.-Interval]
--------------------------------------------------------------------
-----1.830021-------17.14--------0.0000---------1.368687---2.446853
--------------------------------------------------------------------
Test of homogeneity of ORs (approx): chi2(1) = 1.20
-------------------------------------Pr>chi2 = 0.2738
I am interested in printing 0.2738 following Pr>chi2 to excel using putexcel.
However, when I look in the stored results, I cannot for the life of me find 0.2738. I tried the following:
. return list
scalars:
-----------------r(or) = 1.830020550258453
--------------r(lb_or) = 1.368686911814271
--------------r(ub_or) = 2.446852662548659
---------------r(chi2) = 17.13954403711002
------------------r(p) = .0000347316054389
-------------r(df_hom) = 1
-----------r(chi2_hom) = 1.197797587679831
I can't find it. Any help would be greatly appreciated!
Thanks in advance.
—Trenton
Trenton Dailey-Chwalibóg, M.P.H., Ph.D. Fellow
Research Project Manager, Action Against Hunger, France
Research Engineer, AgroParisTech/UGent
I am having trouble retrieving the Pr>chi2 for a Test of homogeneity of ORs using mhodds. To clearly explain my situation with completely arbitrary data from online:
. sysuse nlsw88
. mhodds never_married collgrad, by(south)
Maximum likelihood estimate of the odds ratio
Comparing collgrad==1 vs. collgrad==0
by south
-------------------------------------------------------------------------------
----south | Odds Ratio---------chi2(1)--------P>chi2---------[95% Conf. Interval]
----------+--------------------------------------------------------------------
--------0-|---2.076923----------15.94---------0.0001---------1.43906----2.99751
--------1-|---1.483171-----------2.63---------0.1050---------0.91809----2.39607
-------------------------------------------------------------------------------
----Mantel-Haenszel estimate controlling for south
--------------------------------------------------------------------
-----Odds Ratio----chi2(1)-------P>chi2---------[95% Conf.-Interval]
--------------------------------------------------------------------
-----1.830021-------17.14--------0.0000---------1.368687---2.446853
--------------------------------------------------------------------
Test of homogeneity of ORs (approx): chi2(1) = 1.20
-------------------------------------Pr>chi2 = 0.2738
I am interested in printing 0.2738 following Pr>chi2 to excel using putexcel.
However, when I look in the stored results, I cannot for the life of me find 0.2738. I tried the following:
. return list
scalars:
-----------------r(or) = 1.830020550258453
--------------r(lb_or) = 1.368686911814271
--------------r(ub_or) = 2.446852662548659
---------------r(chi2) = 17.13954403711002
------------------r(p) = .0000347316054389
-------------r(df_hom) = 1
-----------r(chi2_hom) = 1.197797587679831
I can't find it. Any help would be greatly appreciated!
Thanks in advance.
—Trenton
Trenton Dailey-Chwalibóg, M.P.H., Ph.D. Fellow
Research Project Manager, Action Against Hunger, France
Research Engineer, AgroParisTech/UGent
