Re: st: punaf error, even with version control

[Roger Newson <r.newson@imperial.ac.uk>]

Thanks to Andrew for alerting me to this catastrophe.

I have reproduced the results that Andrew describes in my 
up-to-the-minute Stata/MP Version 12, with and without first typing

version 11.0

although the examples work perfectly well in Stata/IC Version 11, which 
I also have.

Using Ben Jann's very useful package -tr-, I have traced the problem to 
-nlcom-. When I type, in Stata 12,

use http://www.stata-press.com/data/r11/lbw.dta, clear
describe
logit low i.race i.smoke, or robust
tr: punaf, at(smoke=0) eform

I find that the failure occurs in the program -nlcom-, when it tries to 
execute the line of code

scalar __000007 = log(_b[_cons])

after which Stata complains (in red error letters):

[_cons] not found

Something very strange is happening here inside -nlcom-, especially as 
_b[_cons] appears to be negative (unsurprisingly, as it is a log 
baseline odds from a logit regression). If there is no advice from 
Statalist, and if I cannot work out what is going on myself and do a 
work-around, then I think I will have to report this to Stata Tech Support.

Best wishes

Roger


Roger B Newson BSc MSc DPhil
Lecturer in Medical Statistics
Respiratory Epidemiology and Public Health Group
National Heart and Lung Institute
Imperial College London
Royal Brompton Campus
Room 33, Emmanuel Kaye Building
1B Manresa Road
London SW3 6LR
UNITED KINGDOM
Tel: +44 (0)20 7352 8121 ext 3381
Fax: +44 (0)20 7351 8322
Email: r.newson@imperial.ac.uk
Web page: http://www.imperial.ac.uk/nhli/r.newson/
Departmental Web page:
http://www1.imperial.ac.uk/medicine/about/divisions/nhli/respiration/popgenetics/reph/

Opinions expressed are those of the author, not of the institution.

On 03/10/2011 13:07, graupel75@gmail.com wrote:
> Hello everyone,
>
> I am trying to use _punaf_ in Stata 12; it worked fine in v11, but
> even the examples fail in v12, giving
>
> [_cons] not found
>
> I have tried version control without success- I suspect the logit
> estimates are simply un-findable by _punaf_ ?
> Any ideas would be greatly appreciated.
>
> cheers-
> Andrew
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/