 |
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: st: Re: problem with -artsurv-
Careful, though, if it has been -off- for a long time you may get surprising
(for you) behaviour in moments when it really matters. Make sure you are on
high alert or -set varabbrev off- again after you are done with -artsurv-
You might also like http://www.stata-journal.com/article.html?article=dm0016
HTH
Martin
_______________________
----- Original Message -----
From: "Michael McCulloch" <[email protected]>
To: <[email protected]>
Sent: Friday, October 24, 2008 11:44 PM
Subject: Re: st: Re: problem with -artsurv-
Brilliant! Thank you kindly.
varabbrev on- and the -program- will run :-)
HTH
Martin
_______________________
----- Original Message ----- From: "Michael McCulloch" <[email protected]>
To: <[email protected]>
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" <[email protected]>
To: <[email protected]>
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" <[email protected]>
To: "Statalist" <[email protected]>
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
[email protected]
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/
*
* 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/