A randomised trail studied the effect of a drug on mortality in low-birth-weight babies. The drug had no effect but, if it *had* been effective, how would I use Stata 15 (Windows) to test whether the effect of the drug in babies weighing <1000 g (baseline) differed from (1) the effect in babies weighing 1000-1499 g and (2) from the effect in babies weighing 1500-1999 g?
The Stata file has
died: 0 (survived), 1 (died)
drug: 0 (control), 1 (treatment)
wtgp: 0 (0-999g), 1 (1000-1499g), 2 (1500-1999g).
Using the Excel spreadsheet for ratios from http://www.statpages.info/Confidence...0P-values.xlsx
wtgp 0 vs 1: 0.86 (0.72-1.04) vs 1.03 (0.83-1.29) gives p = 0.22.
wtgp 0 vs 2: 0.86 (0.72-1.04) vs 0.95 (0.69-1.33) gives p = 0.60.
P = 0.22 (from Excel) differs from p = 0.120 (margins)
P = 0.60 (from Excel) differs from p = 0.133 (margins).
What is the correct way to compare the effect on drug of wtgp 0 vs 1, and wtgp 0 vs 2?
I apologise if this problem is very basic, but I’ve spent many hours trying to solve it.
The Stata file has
died: 0 (survived), 1 (died)
drug: 0 (control), 1 (treatment)
wtgp: 0 (0-999g), 1 (1000-1499g), 2 (1500-1999g).
Code:
. cs died drug, by(wtgp) wtgp | RR [95% Conf. Interval] M-H Weight -----------------+------------------------------------------------- 0-999g | .8643293 .7181096 1.040322 32.19632 1000-1499g | 1.034783 .831774 1.287339 57.5 1500-1999g | .9549525 .6873575 1.326725 34.01655 ------------------------------------------------------------------------------
wtgp 0 vs 1: 0.86 (0.72-1.04) vs 1.03 (0.83-1.29) gives p = 0.22.
wtgp 0 vs 2: 0.86 (0.72-1.04) vs 0.95 (0.69-1.33) gives p = 0.60.
Code:
. logistic died i.drug##i.wtgp ------------------------------------------------------------------------------ died | Odds Ratio Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- 1.drug | .5721154 .2068869 -1.54 0.123 .281627 1.162232 wtgp | 1 | .0979067 .0287701 -7.91 0.000 .0550408 .1741566 2 | .0188347 .0056566 -13.23 0.000 .0104549 .0339312 drug#wtgp | 1 1 | 1.832185 .7197436 1.54 0.123 .8483771 3.956852 1 2 | 1.663846 .6716229 1.26 0.207 .7542577 3.670342 _cons | 3.764706 1.02721 4.86 0.000 2.205355 6.426634 ------------------------------------------------------------------------------ . margins wtgp, dydx(drug) pwcompare(effects) Pairwise comparisons of conditional marginal effects Model VCE : OIM Expression : Pr(died), predict() dy/dx w.r.t. : 1.drug ------------------------------------------------------------------------------ | Contrast Delta-method Unadjusted Unadjusted | dy/dx Std. Err. z P>|z| [95% Conf. Interval] -------------+---------------------------------------------------------------- 0.drug | (base outcome) -------------+---------------------------------------------------------------- 1.drug | wtgp | 1 vs 0 | .1165643 .074964 1.55 0.120 -.0303625 .2634911 2 vs 0 | .1042139 .0693238 1.50 0.133 -.0316582 .2400861 2 vs 1 | -.0123504 .0323953 -0.38 0.703 -.075844 .0511433 ------------------------------------------------------------------------------
P = 0.60 (from Excel) differs from p = 0.133 (margins).
What is the correct way to compare the effect on drug of wtgp 0 vs 1, and wtgp 0 vs 2?
I apologise if this problem is very basic, but I’ve spent many hours trying to solve it.