Re: st: Re: problem with -artsurv-

[Michael McCulloch <mm@pinest.org>]

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/