Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: RE: problem obtaining hazard estimates


From   "Steichen, Thomas J." <SteichT@rjrt.com>
To   <statalist@hsphsun2.harvard.edu>
Subject   st: RE: problem obtaining hazard estimates
Date   Tue, 22 May 2007 10:33:51 -0400

Steiner,

I just tested the example you gave on my Stata SE V9.2.  
No errors occurred and the graph plotted fine.

A trace shows that __ub is used internally in generating the smoothed hazard function. 

Despite your reload, it appears your version still has an error.

Tom

-----------------------------------
Thomas J. Steichen
steicht@rjrt.com
-----------------------------------
 
-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Steinar Fossedal
Sent: Tuesday, May 22, 2007 10:02 AM
To: statalist@hsphsun2.harvard.edu
Subject: st: problem obtaining hazard estimates

Hi, listers

I am experiencing problems obtaining a Kaplan-Meier estimate of the
hazard function. Using -sts graph, haz- returns an error "Variable __ub
not found" r(111). The error does not, however, appear when requesting
the surivival function of the same data. I can reproduce the error using
the leukemia data from stata-press (see printout below), even after a
complete reinstall of Stata 9. How complete it actually was is hard to
tell, my settings were still identical to what they were prior to the
reinstall even though I renamed the old "C:\Program Files\Stata9-folder"
to Stata9_old. Where else does Stata store information?

Until now I have been using -sts gen haz=h- instead, but today, by a one
time fluke, stata suddenly returned different results compared to
earlier. I would normally be willing to write this off as some fault on
my part had it not been for the disturbing fact that the fluke-results
were a lot more in sync with the predictions from my models. The results
are, however, not reproducable in the same session even after doing no
changes to the dataset.

I am now wondering if the two problems are somehow related, and fear
there is something wrong with the hazard-estimates I have been
recieving. Can someone confirm the problem, or suggest a way to handle
it?

Thanks, 

-Steinar




.. use http://www.stata-press.com/data/r8/leukemia, clear
(Leukemia Remission Study)

.. stset weeks, failure(relapse)

     failure event:  relapse != 0 & relapse < .
obs. time interval:  (0, weeks]
 exit on or before:  failure

------------------------------------------------------------------------
------
       42  total obs.
        0  exclusions
------------------------------------------------------------------------
------
       42  obs. remaining, representing
       30  failures in single record/single failure data
      541  total analysis time at risk, at risk from t =         0
                             earliest observed entry t =         0
                                  last observed exit t =        35

.. sts graph

         failure _d:  relapse
   analysis time _t:  weeks

.. sts graph, haz

         failure _d:  relapse
   analysis time _t:  weeks
variable __ub not found
r(111);

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

-----------------------------------------
CONFIDENTIALITY NOTE: This e-mail message, including any
attachment(s), contains information that may be confidential,
protected by the attorney-client or other legal privileges, and/or
proprietary non-public information. If you are not an intended
recipient of this message or an authorized assistant to an intended
recipient, please notify the sender by replying to this message and
then delete it from your system. Use, dissemination, distribution,
or reproduction of this message and/or any of its attachments (if
any) by unintended recipients is not authorized and may be
unlawful.


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