Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: RE: Bug in tabodds, cornfield


From   Ricardo Ovaldia <[email protected]>
To   [email protected]
Subject   Re: st: RE: Bug in tabodds, cornfield
Date   Tue, 4 Oct 2005 07:16:45 -0700 (PDT)

--- Nick Cox <[email protected]> wrote:

> You need to say more about the data. Otherwise 
> "This seems to be a bug" means nothing except 
> "These results surprise me". 

Please read the entire post before replying. I pasted
the output, both with and without the option
specified, to my original email. I Though that was
clear. Here is another example using the auto data and
commands:

. sysuse auto, clear
. tabodds   foreign  rep78
. tabodds   foreign  rep78, or cornfield

Thank you,
Ricardo.


. sysuse auto, clear
(1978 Automobile Data)

. tabodds   foreign  rep78

--------------------------------------------------------------------------
     rep78  |      cases     controls       odds     
[95% Conf. Interval]
------------+-------------------------------------------------------------
          1 |          0            2    0.00000      
       .         .
          2 |          0            8    0.00000      
       .         .
          3 |          3           27    0.11111      
 0.03371   0.36625
          4 |          9            9    1.00000      
 0.39695   2.51919
          5 |          9            2    4.50000      
 0.97229  20.82703
--------------------------------------------------------------------------
Test of homogeneity (equal odds): chi2(4)  =    26.87
                                  Pr>chi2  =   0.0000

Score test for trend of odds:     chi2(1)  =    23.85
                                  Pr>chi2  =   0.0000


. tabodds   foreign  rep78, or cornfield

--------------------------------------------------------------------------
            |                               odds      
    Cornfield
     rep78  |      cases     controls      ratio     
[95% Conf. Interval]
------------+-------------------------------------------------------------
'.' found where integer expected
r(7);



> 
> Ricardo Ovaldia
>  
> > I get an error when I request cornfield CI in
> > -tabodds-. This seems to be a bug. Can anyone
> please
> > verify. Thank you.
> > 
> > Ricardo.
> >  
> > . tabodds vitb1p1 snp if case==1, or
> > 
> >
>
--------------------------------------------------------------
> > -------------
> >          snp |  Odds Ratio       chi2       P>chi2
>    
> > [95% Conf. Interval]
> >
>
-------------+------------------------------------------------
> > -------------
> >           CC |    1.000000          .           . 
>    
> >         .          .
> >           CT |    0.901799       0.26       0.6086
>    
> >  0.607030   1.339707
> >           TT |    0.663223       2.01       0.1564
>    
> >  0.374388   1.174892
> >
>
--------------------------------------------------------------
> > -------------
> > Test of homogeneity (equal odds): chi2(2)  =    
> 2.02
> >                                   Pr>chi2  =  
> 0.3638
> > 
> > Score test for trend of odds:     chi2(1)  =    
> 1.79
> >                                   Pr>chi2  =  
> 0.1811
> > 
> > . tabodds vitb1p1 snp if case==1, or cornf
> > 
> >
>
--------------------------------------------------------------
> > ------------
> >             |                               odds  
>    
> >     Cornfield
> >        snp  |      cases     controls      ratio  
>   
> > [95% Conf. Interval]
> >
>
------------+-------------------------------------------------
> > ------------
> > '.' found where integer expected
> > r(7);
> 
> *
> *   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/
> 


Ricardo Ovaldia, MS
Statistician 
Oklahoma City, OK


	
		
______________________________________________________ 
Yahoo! for Good 
Donate to the Hurricane Katrina relief effort. 
http://store.yahoo.com/redcross-donate3/ 

*
*   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/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index