Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: Re: problem with -artsurv-


From   Michael McCulloch <mm@pinest.org>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: Re: problem with -artsurv-
Date   Fri, 24 Oct 2008 14:44:56 -0700

Brilliant! Thank you kindly.

varabbrev on- and the -program- will run :-)

HTH
Martin
_______________________
----- Original Message ----- From: "Michael McCulloch" <mm@pinest.org>
To: <statalist@hsphsun2.harvard.edu>
Sent: Friday, October 24, 2008 10:57 PM
Subject: Re: st: Re: problem with -artsurv-

Thanks Martin. The location of the error is:

----------------------------------------------------------- begin artsurv._condsrv ---
    - args sf csf
    - tempname small
    - scalar `small'=1e-16
    = scalar __000001=1e-16
    - by grp:gen double `csf'=`sf' if _n==1
    = by grp:gen double cesf0=esf0 if _n==1
    - by grp:replace `csf'=1 if _n>1 & `sf'[_n-1]<`small'
    = by grp:replace cesf0=1 if _n>1 & esf0[_n-1]<__000001
    - by grp:replace `csf'=`sf'/`sf'[_n-1] if _n>1 & `sf'[_n-1]>=`small'
    = by grp:replace cesf0=esf0/esf0[_n-1] if _n>1 & esf0[_n-1]>=__000001
    - replace `csf'=`small' if `csf'<`small'
    = replace cesf0=__000001 if cesf0<__000001
------------------------------------------------------------- end artsurv._condsrv ---
  - gen double el0=-log(cesf)
cesf not found





As I said, -set tr on- and show us the location of the error...

HTH
Martin
_______________________
----- Original Message ----- From: "Michael McCulloch" <mm@pinest.org>
To: <statalist@hsphsun2.harvard.edu>
Sent: Friday, October 24, 2008 10:51 PM
Subject: Re: st: Re: problem with -artsurv-

What I'm trying to figure out is the reason for the error message:

cesf not found
r(111);


Works well for me:

**************
ART - ANALYSIS OF RESOURCES FOR TRIALS (version 1.0.4, 13 January 2005)

A sample size program by Abdel Babiker, Patrick Royston & Friederike Barthel,
MRC Clinical Trials Unit, London NW1 2DA, UK.

Type of trial Superiority - time-to-event outcome
Statistical test assumed                 Unweighted logrank test (local)
Number of groups                         2
Allocation ratio                         Equal group sizes

Total number of periods                  3
Length of each period                    One year

Survival probs per period (group 1)      0.794 0.630 0.500
Survival probs per period (group 2)      0.841 0.707 0.595
Number of recruitment periods            2
Number of follow-up periods              1
Method of accrual                        Uniform
Recruitment period-weights               1 1 0

Hazard ratios as entered (groups 1,2)    1, 0.75
Alpha                                    0.050 (two-sided)
Power (designed)                         0.800

Total sample size (calculated)           1624
Expected total number of events          509

*************************

-set tr on- and show us where the error occurs...

HTH
Martin
_______________________
----- Original Message ----- From: "Michael McCulloch" <mm@pinest.org>
To: "Statalist" <statalist@hsphsun2.harvard.edu>
Sent: Friday, October 24, 2008 10:44 PM
Subject: st: problem with -artsurv-

Hello,
I'm using the program -artsurv- which estimates sample size or power for survival data. When I invoke the command, as follows, using syntax confirmed by the program's author:

. artsurv, method(l) edf0(0.5) hr(1,0.75) np(3) lg(1 2) ldf(.1;.1) wg(1 2) wdf(.3;.1) recrt(2 0,1,0) distant(0) detail(0) onesided(0) ni(0) tunit(1)

I get an error message:

cesf not found
r(111);

Could anyone shed some light on this?

--

Best wishes,
Michael McCulloch



Pine Street Foundation
124 Pine St., San Anselmo, CA 94960-2674
Tel: (415) 407-1357
Fax: (415) 485-1065
mcculloch@pinestreetfoundation.org
www.pinestreetfoundation.org
*
*   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/

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

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

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



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index