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

Re: st: matching Agresti's SAS results with Stata


From   Bill Rising <brising@louisville.edu>
To   "Stata Listserve" <statalist@hsphsun2.harvard.edu>
Subject   Re: st: matching Agresti's SAS results with Stata
Date   Tue, 30 Sep 2003 23:44:35 -0400

On 9/30/03 22:59, Bill Rising wrote

>Hello folks,
>
>Has anyone out there tried mimicking the ordinal logistic regression 
>results from the insomnia example given the second edition of Agresti's 
>book on Categorical analysis?
>
>I was trying to learn to use the gllamm package for Stata using Agresti, 
>because his datasets are all on the web
>http://www.stat.ufl.edu/~aa/cda/sas/sas.html
>(and because I have the book).
>
>As one might guess from the URL, SAS is the package he used to get the 
>results in his book.
>
>I'm stuck trying to mimic table 12.7 on page 514, where he compares 
>estimates for a model using ordinal logistic regression (aka cumulative 
>logit) for some insomnia data (found in table 11.4 on p.462).
[snip...]

>Does anyone know whether SAS computes these models 
>differently? 

Well, I answered this question myself by reading the gllamm manual more 
closely (imagine that!). Using the adaptive quadrature option results in 
coefficients which, though not identical, are very close. This makes 
sense, because adaptive quadrature is used by SAS.

Still stumped on the other two questions.

Bill

*
*   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–2022 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index