Statalist The Stata Listserver


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

st: RE: RE: problem obtaining hazard estimates


From   "Steinar Fossedal" <Steinar.Fossedal@nhh.no>
To   <statalist@hsphsun2.harvard.edu>
Subject   st: RE: RE: problem obtaining hazard estimates
Date   Tue, 22 May 2007 17:28:21 +0200

Thank you for checking, Tom. 

My trace shows that the error when renaming the __ub variable at the
very end of the do-file, but apparently (and not surprisingly) this is
not a program error. Does anyone know how to perform a complete,
remnant-free, uninstallation of Stata so I can start with a clean slate?

Thanks again for helping

-Steinar


(output omitted)

    - if "`by'"!="" {
    = if ""!="" {
      qui replace `id' = _n
      reshape long __y __yy __ub __lb __by, i(`id')
      }
    - drop `time' `deltah' `deltaVh'
    = drop _t __000007 __000008
    - if "`by'"!="" {
    = if ""!="" {
      drop `by'
      rename __y `deltah'
      rename __yy `deltaVh'
      rename __by `by'
      rename __ub `ub'
      rename __lb `lb'
      }
    - else {
    - rename __y1 `deltah'
    = rename __y1 __000007
    - rename __yy1 `deltaVh'
    = rename __yy1 __000008
    - rename __ub `ub'
    = rename __ub __00000F
variable __ub not found
      rename __lb `lb'
      }

(output omitted)


-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Steichen,
Thomas J.
Sent: 22. mai 2007 16:34
To: statalist@hsphsun2.harvard.edu
Subject: st: RE: problem obtaining hazard estimates

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/

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