Dear statalisters,
consider the following output
. xi:poisson count i.A i.B i.C i.A*i.C ,cluster(id)
Poisson regression
Number of obs = 32
Wald
chi2(2) = .
Log pseudo-likelihood = -105.04905 Prob
> chi2 = .
(standard errors
adjusted for clustering on id)
------------------------------------------------------------------------------
| Robust
count | Coef. Std. Err. z P>|z|
[95% Conf. Interval]
-------------+----------------------------------------------------------------
_IA_1 | 1.572489 .0286677 54.85 0.000
1.516302 1.628677
_IB_1 | .0204089 .0449014 0.45 0.649
-.0675963 .108414
_IC_1 | 1.254163 .0339171 36.98 0.000
1.187686 1.320639
_IAXC_1_1 | -1.199721 .014818 -80.96 0.000
-1.228763 -1.170678
_cons | 3.228422 .0355479 90.82 0.000
3.158749 3.298095
------------------------------------------------------------------------------
. lincom _IA_1
( 1) [count]_IA_1 = 0
------------------------------------------------------------------------------
count | Coef. Std. Err. z P>|z|
[95% Conf. Interval]
-------------+----------------------------------------------------------------
(1) | 1.572489 .0286677 54.85 0.000
1.516302 1.628677
------------------------------------------------------------------------------
. lincom _IA_1+_IC_1+ _IAXC_1_1
( 1) [count]_IA_1 + [count]_IC_1 + [count]_IAXC_1_1
= 0
------------------------------------------------------------------------------
count | Coef. Std. Err. z P>|z|
[95% Conf. Interval]
-------------+----------------------------------------------------------------
(1) | 1.626931 .0482076 33.75 0.000
1.532446 1.721416
------------------------------------------------------------------------------
. lincom _IC_1+ _IAXC_1_1
( 1) [count]_IC_1 + [count]_IAXC_1_1 = 0
------------------------------------------------------------------------------
count | Coef. Std. Err. z P>|z|
[95% Conf. Interval]
-------------+----------------------------------------------------------------
(1) | .0544419 .0226218 2.41 0.016
.010104 .0987798
------------------------------------------------------------------------------
A, B, C are binary covariates. As you see the upper
limit of the 95%CI in the first "lincom" is greater
than the beta estimate in the second "lincom" and the
lower limit of the 95% CI in the second "lincom" is
lower than the beta estimate in the first "lincom".
Given this overlap I was expecting a non-significant
(at the 5% level) difference between the first two
estimates. However the third "lincom" gives a p=0.016
for the difference of the first two estimates. Is
there something wrong with this and if not how can one
justify the overlaping in the CIs when the difference
is significant. Maybe I am missing something obvious
but I can't find a good explanation.
Thanks in advance for any comments
Nikos Pantazis
Biostatistician
__________________________________
Do you Yahoo!?
Friends. Fun. Try the all-new Yahoo! Messenger.
http://messenger.yahoo.com/
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
