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 <ovaldia@yahoo.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: RE: Bug in tabodds, cornfield
Date   Tue, 4 Oct 2005 07:53:19 -0700 (PDT)

I am sorry Nick. You are correct I should have been
more explicit. I just wanted someone to verify that
they could not get output with the -cornfield- option.

Yes my version is newer:

. which tabodds
C:\Program Files\Stata9\ado\base\t\tabodds.ado
*! version 6.3.5  17mar2005

Ricardo.

--- n j cox <n.j.cox@durham.ac.uk> wrote:

> I did read your entire post, twice. In
> your original mail you said nothing about your data.
> Giving output is one thing, but the data are
> another.
> 
> There is a long history of people claiming bugs
> on the basis of little or no evidence that can
> be discussed. It is good advice not to claim
> a bug unless you can identify the code segment
> that is problematic. This does not require
> a reporter to be a Stata programmer. Sometimes
> users can -set trace on- and see where code
> stops with an error.
> 
> This reply mail is concrete and lets someone say
> something
> more helpful.
> 
> At this moment I am on a computer away from base
> which
> has just Stata 8.2. I get different output.
> 
> . tabodds   foreign  rep78, or cornfield
> 
>
--------------------------------------------------------------------------
>              |                               odds   
>        Cornfield
>       rep78  |      cases     controls      ratio   
>   [95% Conf. Interval]
>
------------+-------------------------------------------------------------
>            1 |          0            2          .   
>           .         .
>            2 |          0            8          .   
>           .         .
>            3 |          3           27          .   
>     0.04267         .
>            4 |          9            9          .   
>     0.43332         .
>            5 |          9            2          .   
>     1.51393         .
>
--------------------------------------------------------------------------
> 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
> 
> . which tabodds
> *! version 6.2.6  22apr2003
> 
> Does your problem arise with a later or earlier
> version of
> -tabodds-?
> 
> Nick
> n.j.cox@durham.ac.uk
> 
> Ricardo Ovaldia
> 
> --- Nick Cox <n.j.cox@durham.ac.uk> 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
> 
> . 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! Mail - PC Magazine Editors' Choice 2005 
http://mail.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/



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