Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: Potential bug in Stata 11
From
Philip Ender <[email protected]>
To
[email protected]
Subject
Re: st: Potential bug in Stata 11
Date
Wed, 19 May 2010 08:53:02 -0700
Babak Oskooei <[email protected]> wrote:
I have been conducting some simulations in Stata and faced a
computational problem in Stata 11. I run the same code in Stata 9 and
it gives me the right answer. I would be grateful for any feedback on
this. Below, I include the code and the output from both version 9 and
version 11. I am simply simulating conditional time-to-event data that
are exponentially distributed with one normally distributed
covariate. ....
--------
Using Stata 11 running on Mac OX 10.6.3, I commented out the line that
generates x like this:
*gen x = uniform()
The code ran fine with these results:
Iteration 0: log likelihood = -82108.928
Iteration 1: log likelihood = -71818.568
Iteration 2: log likelihood = -69286.528
Iteration 3: log likelihood = -69035.66
Iteration 4: log likelihood = -69032.801
Iteration 5: log likelihood = -69032.627
Refining estimates:
Iteration 0: log likelihood = -69032.613
Iteration 1: log likelihood = -69032.613
Cox regression -- no ties
No. of subjects = 10000 Number of obs = 10000
No. of failures = 10000
Time at risk = 1679070425
LR chi2(1) = 26152.63
Log likelihood = -69032.613 Prob > chi2 = 0.0000
------------------------------------------------------------------------------
_t | Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
y | 4.177157 .0343322 121.67 0.000 4.109867 4.244447
------------------------------------------------------------------------------
If the lise is not commented out then I get the same error message as Babak.
--
Phil Ender
UCLA Statistical Consulting Group
*
* 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/