Hi all,
I am trying to run a seemingly unrelated regression based on two ordered probit models with interaction terms.
Do I need to include margins after running suest?
Here is my current code:
svy:oprobit rs1p1_w6 wife_fininf wife_cr i.wife_interact wife_mp wife_cmt prb8p1_w6 racep1_w2 religp1_w3 edu1p1_w6 wrk3p1_w6 inc1p1_w6
estimates store WIFE
margins, dydx(*)
svy:oprobit rs1p2_w6 husband_fininf husband_cr i.husband_interact husband_mp husband_cmt prb8p2_w6 racep2_w2 religp2_w3 edu1p2_w6 wrk3p2_w6 inc1p2_w6
estimates store HUSBAND
margins, dydx(*)
suest WIFE HUSBAND
The code works fine if that is all I do, but I think I might need to run "margins, dydx(*)" after the suest command and when I do that here is the output:
Help?
----------------------------------------------------------------------------------
| Delta-method
| dy/dx std. err. t P>|t| [95% conf. interval]
-----------------+----------------------------------------------------------------
wife_fininf | -4.179195 .382902 -10.91 0.000 -4.930424 -3.427966
wife_cr | .0666802 .0173444 3.84 0.000 .0326516 .1007088
|
wife_interact |
1 | -1.214107 .290087 -4.19 0.000 -1.783238 -.6449749
2 | 4.692742 .6423256 7.31 0.000 3.432541 5.952942
3 | 3.2315 .6572719 4.92 0.000 1.941975 4.521024
4 | -.7777807 .27968 -2.78 0.006 -1.326495 -.2290668
5 | 3.994617 .7977942 5.01 0.000 2.429396 5.559837
6 | 4.909395 .4040208 12.15 0.000 4.116733 5.702058
7 | 3.449989 .3839863 8.98 0.000 2.696633 4.203345
8 | 3.99036 .4934313 8.09 0.000 3.02228 4.95844
9 | 4.305126 .3618916 11.90 0.000 3.595118 5.015133
10 | 3.855901 .4935441 7.81 0.000 2.8876 4.824203
11 | 3.772383 .3795927 9.94 0.000 3.027647 4.517119
12 | 4.256579 .4556861 9.34 0.000 3.362553 5.150605
13 | 3.858597 .3658492 10.55 0.000 3.140824 4.576369
14 | 4.495018 .3720911 12.08 0.000 3.765 5.225037
15 | 4.270707 .4037713 10.58 0.000 3.478534 5.06288
16 | 4.021922 .4131977 9.73 0.000 3.211255 4.832588
17 | 4.257633 .397539 10.71 0.000 3.477688 5.037579
18 | 4.071482 .423954 9.60 0.000 3.239712 4.903252
19 | 4.095519 .426118 9.61 0.000 3.259504 4.931535
20 | 3.754885 .4237966 8.86 0.000 2.923424 4.586346
21 | 4.097495 .4400121 9.31 0.000 3.23422 4.96077
22 | 4.160238 .4560379 9.12 0.000 3.265522 5.054955
23 | 3.905599 .4986291 7.83 0.000 2.927321 4.883876
24 | 3.919222 .5023273 7.80 0.000 2.933689 4.904755
|
wife_mp | .0843045 .0127529 6.61 0.000 .0592841 .1093248
wife_cmt | .0695423 .0084608 8.22 0.000 .0529428 .0861419
prb8p1_w6 | -.1695322 .0650855 -2.60 0.009 -.2972257 -.0418387
racep1_w2 | .0334887 .0348548 0.96 0.337 -.0348942 .1018715
religp1_w3 | -.0086093 .0115435 -0.75 0.456 -.0312568 .0140382
edu1p1_w6 | -.0140671 .0364919 -0.39 0.700 -.0856618 .0575276
wrk3p1_w6 | .0105946 .0188262 0.56 0.574 -.0263412 .0475304
inc1p1_w6 | -.0109769 .0106545 -1.03 0.303 -.0318802 .0099265
husband_fininf | 0 (omitted)
husband_cr | 0 (omitted)
|
husband_interact |
1 | 0 (omitted)
2 | 0 (omitted)
3 | 0 (omitted)
4 | 0 (omitted)
5 | 0 (omitted)
6 | 0 (omitted)
7 | 0 (omitted)
8 | 0 (omitted)
9 | 0 (omitted)
10 | 0 (omitted)
11 | 0 (omitted)
12 | 0 (omitted)
13 | 0 (omitted)
14 | 0 (omitted)
15 | 0 (omitted)
16 | 0 (omitted)
17 | 0 (omitted)
18 | 0 (omitted)
19 | 0 (omitted)
20 | 0 (omitted)
21 | 0 (omitted)
22 | 0 (omitted)
23 | 0 (omitted)
24 | 0 (omitted)
|
husband_mp | 0 (omitted)
husband_cmt | 0 (omitted)
prb8p2_w6 | 0 (omitted)
racep2_w2 | 0 (omitted)
religp2_w3 | 0 (omitted)
edu1p2_w6 | 0 (omitted)
wrk3p2_w6 | 0 (omitted)
inc1p2_w6 | 0 (omitted)
----------------------------------------------------------------------------------
Note: dy/dx for factor levels is the discrete change from the base level.
I am trying to run a seemingly unrelated regression based on two ordered probit models with interaction terms.
Do I need to include margins after running suest?
Here is my current code:
svy:oprobit rs1p1_w6 wife_fininf wife_cr i.wife_interact wife_mp wife_cmt prb8p1_w6 racep1_w2 religp1_w3 edu1p1_w6 wrk3p1_w6 inc1p1_w6
estimates store WIFE
margins, dydx(*)
svy:oprobit rs1p2_w6 husband_fininf husband_cr i.husband_interact husband_mp husband_cmt prb8p2_w6 racep2_w2 religp2_w3 edu1p2_w6 wrk3p2_w6 inc1p2_w6
estimates store HUSBAND
margins, dydx(*)
suest WIFE HUSBAND
The code works fine if that is all I do, but I think I might need to run "margins, dydx(*)" after the suest command and when I do that here is the output:
Help?
----------------------------------------------------------------------------------
| Delta-method
| dy/dx std. err. t P>|t| [95% conf. interval]
-----------------+----------------------------------------------------------------
wife_fininf | -4.179195 .382902 -10.91 0.000 -4.930424 -3.427966
wife_cr | .0666802 .0173444 3.84 0.000 .0326516 .1007088
|
wife_interact |
1 | -1.214107 .290087 -4.19 0.000 -1.783238 -.6449749
2 | 4.692742 .6423256 7.31 0.000 3.432541 5.952942
3 | 3.2315 .6572719 4.92 0.000 1.941975 4.521024
4 | -.7777807 .27968 -2.78 0.006 -1.326495 -.2290668
5 | 3.994617 .7977942 5.01 0.000 2.429396 5.559837
6 | 4.909395 .4040208 12.15 0.000 4.116733 5.702058
7 | 3.449989 .3839863 8.98 0.000 2.696633 4.203345
8 | 3.99036 .4934313 8.09 0.000 3.02228 4.95844
9 | 4.305126 .3618916 11.90 0.000 3.595118 5.015133
10 | 3.855901 .4935441 7.81 0.000 2.8876 4.824203
11 | 3.772383 .3795927 9.94 0.000 3.027647 4.517119
12 | 4.256579 .4556861 9.34 0.000 3.362553 5.150605
13 | 3.858597 .3658492 10.55 0.000 3.140824 4.576369
14 | 4.495018 .3720911 12.08 0.000 3.765 5.225037
15 | 4.270707 .4037713 10.58 0.000 3.478534 5.06288
16 | 4.021922 .4131977 9.73 0.000 3.211255 4.832588
17 | 4.257633 .397539 10.71 0.000 3.477688 5.037579
18 | 4.071482 .423954 9.60 0.000 3.239712 4.903252
19 | 4.095519 .426118 9.61 0.000 3.259504 4.931535
20 | 3.754885 .4237966 8.86 0.000 2.923424 4.586346
21 | 4.097495 .4400121 9.31 0.000 3.23422 4.96077
22 | 4.160238 .4560379 9.12 0.000 3.265522 5.054955
23 | 3.905599 .4986291 7.83 0.000 2.927321 4.883876
24 | 3.919222 .5023273 7.80 0.000 2.933689 4.904755
|
wife_mp | .0843045 .0127529 6.61 0.000 .0592841 .1093248
wife_cmt | .0695423 .0084608 8.22 0.000 .0529428 .0861419
prb8p1_w6 | -.1695322 .0650855 -2.60 0.009 -.2972257 -.0418387
racep1_w2 | .0334887 .0348548 0.96 0.337 -.0348942 .1018715
religp1_w3 | -.0086093 .0115435 -0.75 0.456 -.0312568 .0140382
edu1p1_w6 | -.0140671 .0364919 -0.39 0.700 -.0856618 .0575276
wrk3p1_w6 | .0105946 .0188262 0.56 0.574 -.0263412 .0475304
inc1p1_w6 | -.0109769 .0106545 -1.03 0.303 -.0318802 .0099265
husband_fininf | 0 (omitted)
husband_cr | 0 (omitted)
|
husband_interact |
1 | 0 (omitted)
2 | 0 (omitted)
3 | 0 (omitted)
4 | 0 (omitted)
5 | 0 (omitted)
6 | 0 (omitted)
7 | 0 (omitted)
8 | 0 (omitted)
9 | 0 (omitted)
10 | 0 (omitted)
11 | 0 (omitted)
12 | 0 (omitted)
13 | 0 (omitted)
14 | 0 (omitted)
15 | 0 (omitted)
16 | 0 (omitted)
17 | 0 (omitted)
18 | 0 (omitted)
19 | 0 (omitted)
20 | 0 (omitted)
21 | 0 (omitted)
22 | 0 (omitted)
23 | 0 (omitted)
24 | 0 (omitted)
|
husband_mp | 0 (omitted)
husband_cmt | 0 (omitted)
prb8p2_w6 | 0 (omitted)
racep2_w2 | 0 (omitted)
religp2_w3 | 0 (omitted)
edu1p2_w6 | 0 (omitted)
wrk3p2_w6 | 0 (omitted)
inc1p2_w6 | 0 (omitted)
----------------------------------------------------------------------------------
Note: dy/dx for factor levels is the discrete change from the base level.
