I'm wondering why SAS v9.4 proc seqdesign is giving me a substantively different sample size than EAST6.  The clinical trial design for both is as follows:   Type I error:  0.05 Type II error:  0.10 one sample test of proportions:  null = 0.10, alternative = 0.265 upper test of the alternative (proportion > 0.265 reference proportion is the null (proportion = 0.10) stop for both futility or efficacy; binding futility boundary two stages using an error spending group sequential design with an O'Brien-Fleming type spending function interim analysis at 50% of information   I get n=29 with SAS v9.4 proc seqdesign with these design parameters; EAST6 gives n=42.  I've pasted in the SAS code I used but I do not have the EAST6 equivalent.  I don't have access to EAST6 to investigate this, either.   Any help vastly appreciated! swannie 　 proc seqdesign plots=boundary (hscale=samplesize) boundaryscale=mle errspend; OneSidedOBrienFleming: design nstages=2 method=errspendobf alt=upper stop=both (betaboundary=binding) alpha=0.05 beta=0.10 info=equal; samplesize model(ceiladjdesign=include) =onesamplefreq(nullprop=0.10 prop=0.265 ref=nullprop); ods output AdjustedBoundary=Bnd_Prop4; run;   ... View more

Hi,   I'm using SAS 9.4 with the attached code and data.  I would like to be able to set the order of the legend items to be something other than alphabetic (which is not meaningful in this context).  It should be Complete Response, Partial Response, Unconfirmed Partial Response, Stable Disease, and Progressive Disease.  My last attempt was to give the responses numeric values with formats, but the legend uses the formats and not the numeric values to determine the order.  There are some workarounds for other types of graphs, but I didn't find anything for a waterfall plot.  Any help vastly appreciated!   The site keeps telling me that the contents of my .sas and .sas7bdat files don't match their file types (????) so I'm pasting in the code and data below. Regards, Suzanne   　   * Set today's date*; data _null_; tday=put(today(),date9.); call symput('today',tday); run;   * Set libname; libname ads204 'path'; libname bor204 'path';   * Assign a macro variable to the data cut-off date; %let cutoff=09APR2018;   　 *******************************************************************************************************************************************; * Assign format for borcd ; *******************************************************************************************************************************************; proc format;     value dos 1='IMO2125 04 mg' 2='IMO2125 08 mg' 3='IMO2125 16 mg' 4='IMO2125 32 mg'   ; value rr 1='Complete Response' 2='Partial Response' 2.5='Unconfirmed Partial Response' 3='Stable Disease' 4='Progressive Disease' 5='Not Evaluable'   ; run;   　 proc sort data=ads204.waterfall out=waterfall; by subject; run; proc sort data=bor204.bor out=bor ; by usubjid; run;   　 ods graphics / reset=all border=off WIDTH=10.0in; ods powerpoint file = "path\Waterfall_Plot_&today..PPTX" ; title ' '; title1 "Maximum Percentage Decrease in Target Lesion Diameters"; title2 "Subjects Treated with IMO-2125 8mg + Ipilimumab " ; footnote1 j=left height=1 "Data cut-off date: &cutoff Produced on &SYSDATE9." ; proc sgplot data=waterfall; format borcd rr.; styleattrs datacolors=(CXC0C0C0 CX4682B4 cx7edb8b CXA3E5AC CX19BF30 ) datacontrastcolors=(black); vbar n / response=rschg group=borcd dataskin=sheen; refline 20 -30 / lineattrs=(pattern=shortdash); xaxis display=none; yaxis values=(60 to -100 by -20) grid label="Change from Basline (%)"; keylegend / title='' sortorder=ascending location=inside position=bottomleft across=1 border; run;   　 ods powerpoint close; ods graphics off; proc printto; run;       n borcd RSCHG   1 4 40.0   2 4 21.2   3 3 16.2   4 3 12.6   5 4 12.0   6 4 11.7   7 3 10.3   8 4 9.0   9 3 0.0   10 3 -1.1   11 3 -12.5   12 2 -39.5   13 2 -51.5   14 3 -57.1   15 2 -58.2   16 2 -73.3   17 2 -81.8   18 2.5 -93.2   19 1 -100 ... View more

      Hi, I need some help figuring out how to use the loss=exp(hazard=hazard rate) option in proc seqdesign correctly.  If the trial is expected to have 10% drop out total, but most is early (within 4 months of randomization), what's the proper (or best) way to calculate the hazard rate?  Any help vastly appreciated.   swannie   %let drop_out=0.0027; ods graphics on; proc seqdesign plots=boundary(hscale=samplesize) boundaryscale=stdz errspend ; OneSidedOBrienFleming: design nstages=2 method=errfuncobf alpha=0.02 alt=upper beta=0.10 info=cum(0.50 1) stop=both(betaboundary=nonbinding); samplesize model(ceiladjdesign=include)=twosamplesurvival (nullhazard=0.0608 hazard=0.0385 weight= 1 1 accrual=uniform acctime=&rectime. loss=exp(hazard=&droprate.) totaltime=&total_time. ceiling=n); run; ods graphics off;   ... View more

I'm using proc sgpanel to create scatterplots with regression lines for 2 variables in an overlay.  I've jittered the data to see it all but the resulting plot has, in addition to the filled red squares and black circles that I've specified mystery data (small, unfilled circles).  What are these, why are they there, and most importantly, how do I make them go away? Code included;  *Assign libraries and formats; libname in 'U:\'; proc format ; value tmm 1='predose' 2='24hr' 8='C3W8' ; =run; options orientation=landscape nodate nonumber; ods pdf file="U:\test_plot.pdf" style=styles.arial; ods graphics / reset=all border=off HEIGHT = 193.8867mm WIDTH = 266.7mm border=off; title ''; proc sgpanel data=in.check noautolegend; panelby respond /rows=1 columns=1; scatter x=time2 y=ki67t_pct / jitter jitterwidth=1 markerattrs=(symbol=circlefilled size=12 color=black ); reg x=time2 y=ki67t_pct; scatter x=time2 y=ki67b_pct / jitter markerattrs=(symbol=squarefilled size=12 color=red ); reg x=time2 y=ki67b_pct; where respond='N'; rowaxis label ='% KI67+ on CD8+ T Cells' labelattrs=(size=14 weight=bold) VALUEATTRS= (Weight=Bold) values=(0 to 100 by 20) offsetmin=0.10; colaxis label=' ' labelattrs=(size=14 weight=bold) VALUEATTRS= (size=10 Weight=Bold) label=" " values=(1, 2, 8) offsetmin=0.15; format time2 tmm.; run; ods pdf close; ods graphics off; data, and pdf of the plot attached.     Any help vastly appreciated. ... View more