data a;   input z p q e;   cards; 1.96 0.95 0.05 0.05 ; data b;   set a;   n=(z*z*p*q)/(e*e);   run;   proc print; run; proc power;    onesamplefreq test=z method=normal       sides          = 2       alpha          = 0.05       nullproportion = 0.90       proportion     = 0.95       ntotal         = .       power          = 0.9; run; ... View more

I'm trying to program sample size for a bayesian adaptive design, with a binary performance goal of say, AE rate less than 20%.  Data will be analyzed after X subjects, and if needed after each subsequent Y subjects.  There are three possible outcomes after each stage until success, futility or max sample size of N is reached.  So I need to program simulations and probability required for determining study success at each analysis stage. Note that success decision can be something like: "Termination of enrollment for success if the posterior probability using observed data was sufficient to meet the performance goal (99% probability or greater for 50 and 75 subject interim checks, 97.5% thereafter)" Continued accrual due to reasonable predictive probability (10% or greater) that performance goal can be met with more subjects. Futility if predictive probability the performance goal could be met with enrollment of more subjects was less than 10% Thank you for any insight.... ... View more

want to left justify the column headers for the variables 'characteristic' and 'stat_cat_'.  As you can see from my output below, this is NOT happening - see my code after the output.  Please help! Here's my output:              Disposition                              Value Group A Group B Total (n=7) (n=6) (n=13) Randomized Yes 7 (100.0 %) 6 (100.0 %) 13 (100.0 %) Prematurely Discontinued No 4 (57.1 %) 3 (50.0 %) 7 (53.8 %) Yes 3 (42.9 %) 3 (50.0 %) 6 (46.2 %) Enrolled in Next One No 1 (14.3 %) 1 (16.7 %) 2 (15.4 %) No 2 (28.6 %) 1 (16.7 %) 3 (23.1 %) Missing Values 4 (57.1 %) 4 (66.7 %) 8 (61.5 %) Reason Discontinued: Adverse Event 0 (0.0%) 1 (33.3 %) 1 (16.7 %) Ulcer 2 (66.7 %) 0 (0.0%) 2 (33.3 %) Lost to follow-up 1 (33.3 %) 0 (0.0%) 1 (16.7 %) Other 0 (0.0%) 1 (33.3 %) 1 (16.7 %) 0 (0.0%) 1 (33.3 %) 1 (16.7 %) Here's my code: /** RTF **/ %RTFOPTNS; %OUTFILE(NAME=T1_1_SubjectDisposition); ** Name the Output file ** ; OPTION MISSING=' '; ODS LISTING CLOSE; ODS RTF FILE=RTFFILE STYLE=STYLES.RTFSTYLE BODYTITLE; ODS NORESULTS; PROC REPORT DATA=final NOWD HEADLINE HEADSKIP SPLIT='|' SPACING=2 MISSING ; COLUMN characteristic stat_cat_  var1 var2  var99;     DEFINE characteristic/DISPLAY WIDTH=15 left "Disposition"             STYLE(COLUMN)=[CELLWIDTH=1.25IN ];     DEFINE stat_cat_/DISPLAY WIDTH=15 left "Value"             STYLE(COLUMN)=[CELLWIDTH=3.75IN ];     DEFINE var1/DISPLAY WIDTH=15 CENTER "PluroGel N|(n=&grpPRand.)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     DEFINE var2/DISPLAY WIDTH=15 CENTER "Vehicle|(n=&grpPRand.)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     DEFINE var99/DISPLAY WIDTH=15 CENTER "Total|(n=&grpTotRand)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     %RTFTITLES(t3=Table 14.1.1 Subject Disposition); RUN; OPTIONS MISSING='.'; FOOTNOTE; TITLE; ODS RTF CLOSE; ODS pdf CLOSE; ODS LISTING ; PROC PRINTTO; RUN; ... View more

I want to left justify the column headers for the variables 'characteristic' and 'stat_cat_'.  As you can see from my output below, this is NOT happening - see my code after the output.  Please help! Here's my output:              Disposition                              Value Group A Group B Total (n=7) (n=6) (n=13) Randomized Yes 7 (100.0 %) 6 (100.0 %) 13 (100.0 %) Prematurely Discontinued No 4 (57.1 %) 3 (50.0 %) 7 (53.8 %) Yes 3 (42.9 %) 3 (50.0 %) 6 (46.2 %) Enrolled in Next One No 1 (14.3 %) 1 (16.7 %) 2 (15.4 %) No 2 (28.6 %) 1 (16.7 %) 3 (23.1 %) Missing Values 4 (57.1 %) 4 (66.7 %) 8 (61.5 %) Reason Discontinued: Adverse Event 0 (0.0%) 1 (33.3 %) 1 (16.7 %) Ulcer 2 (66.7 %) 0 (0.0%) 2 (33.3 %) Lost to follow-up 1 (33.3 %) 0 (0.0%) 1 (16.7 %) Other 0 (0.0%) 1 (33.3 %) 1 (16.7 %) 0 (0.0%) 1 (33.3 %) 1 (16.7 %) Here's my code: /** RTF **/ %RTFOPTNS; %OUTFILE(NAME=T1_1_SubjectDisposition); ** Name the Output file ** ; OPTION MISSING=' '; ODS LISTING CLOSE; ODS RTF FILE=RTFFILE STYLE=STYLES.RTFSTYLE BODYTITLE; ODS NORESULTS; PROC REPORT DATA=final NOWD HEADLINE HEADSKIP SPLIT='|' SPACING=2 MISSING ; COLUMN characteristic stat_cat_  var1 var2  var99;     DEFINE characteristic/DISPLAY WIDTH=15 left "Disposition"             STYLE(COLUMN)=[CELLWIDTH=1.25IN ];     DEFINE stat_cat_/DISPLAY WIDTH=15 left "Value"             STYLE(COLUMN)=[CELLWIDTH=3.75IN ];     DEFINE var1/DISPLAY WIDTH=15 CENTER "PluroGel N|(n=&grpPRand.)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     DEFINE var2/DISPLAY WIDTH=15 CENTER "Vehicle|(n=&grpPRand.)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     DEFINE var99/DISPLAY WIDTH=15 CENTER "Total|(n=&grpTotRand)"             STYLE(COLUMN)=[CELLWIDTH=1.0IN ];     %RTFTITLES(t3=Table 14.1.1 Subject Disposition); RUN; OPTIONS MISSING='.'; FOOTNOTE; TITLE; ODS RTF CLOSE; ODS pdf CLOSE; ODS LISTING ; PROC PRINTTO; RUN; ... View more

How can I use proc power to calculate the sample size needed for the following scenario: 2 independent groups 1:2 randomization 2-sided 95% confidence interval lower bound of the mean difference must exceed +2,SD=3.1   hel How can I specify the lower bound in the SAS code?   Thanks! ... View more

In phreg, using the output statement with survival= option, what is the formula to recreate the survival= calculation?  I cannot recreate it. Thank you very much! ... View more