The SEQDESIGN Procedure

Example 103.5 Creating Designs Using Haybittle-Peto Methods

(View the complete code for this example.)

This example requests two 3-stage group sequential designs for normally distributed statistics. Each design uses a Haybittle-Peto method with a two-sided alternative and early stopping to reject the hypothesis. One design uses the specified interim boundary Z values and derives the final-stage boundary value for the specified and errors. The other design uses the specified boundary Z values and derives the overall and errors.

The following statements specify the interim boundary Z values and derive the final-stage boundary value for the specified and :

ods graphics on;
proc seqdesign altref=0.25
               errspend
               stopprob
               plots=errspend
               ;
   OneSidedPeto: design nstages=3
                 method=peto( z=3)
                 alt=upper   stop=reject
                 alpha=0.05  beta=0.10;
run;

The "Design Information" table in Output 103.5.1 displays design specifications and maximum information in percentage of its corresponding fixed-sample design.

Output 103.5.1: Haybittle-Peto Design Information

The SEQDESIGN Procedure
Design: OneSidedPeto

Design Information
Statistic DistributionNormal
Boundary ScaleStandardized Z
Alternative HypothesisUpper
Early StopReject Null
MethodHaybittle-Peto
Boundary KeyBoth
Alternative Reference0.25
Number of Stages3
Alpha0.05
Beta0.1
Power0.9
Max Information (Percent of Fixed Sample)100.2466
Max Information137.3592
Null Ref ASN (Percent of Fixed Sample)100.1192
Alt Ref ASN (Percent of Fixed Sample)87.35


The "Method Information" table in Output 103.5.2 displays the and errors and the derived drift parameter, which is the standardized alternative reference at the final stage.

Output 103.5.2: Method Information

Method Information
BoundaryMethodAlphaBetaAlternative
Reference
Drift
Upper AlphaHaybittle-Peto0.050000.100000.252.930009


With the STOPPROB option, the "Expected Cumulative Stopping Probabilities" table in Output 103.5.3 displays the expected stopping stage and cumulative stopping probability to reject the null hypothesis at each stage under various hypothetical references , where is the alternative reference and are the default values in the CREF= option.

Output 103.5.3: Stopping Probabilities

Expected Cumulative Stopping Probabilities
Reference = CRef * (Alt Reference)
CRefExpected
Stopping Stage
SourceStopping Probabilities
Stage_1Stage_2Stage_3
0.00002.996Reject Null0.001350.002460.05000
0.50002.941Reject Null0.015610.043720.42762
1.00002.614Reject Null0.095380.290570.90000
1.50001.944Reject Null0.321850.734420.99698


The "Boundary Information" table in Output 103.5.4 displays information level, alternative references, and boundary values. By default (or equivalently if you specify BOUNDARYSCALE=STDZ), the standardized Z scale is used to display the alternative references and boundary values. The resulting standardized alternative reference at stage k is given by , where is the alternative reference and is the information level at stage k, .

Output 103.5.4: Boundary Information

Boundary Information (Standardized Z Scale)
Null Reference = 0
_Stage_ AlternativeBoundary Values
Information LevelReferenceUpper
ProportionActualUpperAlpha
10.333345.78641.691643.00000
20.666791.572812.392343.00000
31.0000137.35922.930011.65042


At each interim stage, if the standardized statistic , the trial is stopped and the null hypothesis is rejected. If the statistic , the trial continues to the next stage. At the final stage, the null hypothesis is rejected if the statistic . Otherwise, the hypothesis is accepted. Note that the boundary values at the final stage, 1.65, are close to the critical values 1.645 in the corresponding fixed-sample design.

The "Error Spending Information" in Output 103.5.5 displays cumulative error spending at each stage for each boundary. The stage 1 spending 0.00135 corresponds to the one-sided p-value for a standardized Z statistic, .

Output 103.5.5: Error Spending Information

Error Spending Information
_Stage_Information
Level
Cumulative Error Spending
Upper
ProportionBetaAlpha
10.33330.000000.00135
20.66670.000000.00246
31.00000.100000.05000


With ODS Graphics enabled, a detailed boundary plot with the rejection and acceptance regions is displayed, as shown in Output 103.5.6. With the STOP=REJECT option, the interim rejection boundaries are displayed.

Output 103.5.6: Boundary Plot

Boundary Plot


With the PLOTS=ERRSPEND option, the procedure displays a plot of error spending for each boundary, as shown in Output 103.5.7. The error spending values in the "Error Spending Information" in Output 103.5.4 are displayed in the plot. As expected, the error spending at each of the first two stages is small, with the standardized Z boundary value 3.

Output 103.5.7: Error Spending Plot

Error Spending Plot


The following statements specify the boundary Z values and derive the and errors from these completely specified boundary values:

ods graphics on;
proc seqdesign altref=0.25
               maxinfo=200
               errspend
               stopprob
               plots=errspend
               ;
   OneSidedPeto: design nstages=3
                 method=peto(z=3 2.5 2)
                 alt=upper  stop=reject
                 boundarykey=none
                 ;
run;

The "Design Information" table in Output 103.5.8 displays design specifications and derived and error levels.

Output 103.5.8: Design Information

The SEQDESIGN Procedure
Design: OneSidedPeto

Design Information
Statistic DistributionNormal
Boundary ScaleStandardized Z
Alternative HypothesisUpper
Early StopReject Null
MethodHaybittle-Peto
Boundary KeyNone
Alternative Reference0.25
Number of Stages3
Alpha0.02532
Beta0.06035
Power0.93965
Max Information (Percent of Fixed Sample)101.6769
Max Information200
Null Ref ASN (Percent of Fixed Sample)101.3933
Alt Ref ASN (Percent of Fixed Sample)73.74031


The "Method Information" table in Output 103.5.9 displays the and errors and the derived drift parameter for each boundary.

Output 103.5.9: Method Information

Method Information
BoundaryMethodAlphaBetaAlternative
Reference
Drift
Upper AlphaHaybittle-Peto0.025320.060350.253.535534


With the STOPPROB option, the "Expected Cumulative Stopping Probabilities" table in Output 103.5.10 displays the expected stopping stage and cumulative stopping probability to reject the null hypothesis at each stage under various hypothetical references , where is the alternative reference and are the default values in the CREF= option.

Output 103.5.10: Stopping Probabilities

Expected Cumulative Stopping Probabilities
Reference = CRef * (Alt Reference)
CRefExpected
Stopping Stage
SourceStopping Probabilities
Stage_1Stage_2Stage_3
0.00002.992Reject Null0.001350.007020.02532
0.50002.826Reject Null0.023890.150300.41775
1.00002.176Reject Null0.168840.655440.93965
1.50001.508Reject Null0.524660.967080.99954


The "Boundary Information" table in Output 103.5.11 displays information level, alternative references, and boundary values.

Output 103.5.11: Boundary Information

Boundary Information (Standardized Z Scale)
Null Reference = 0
_Stage_ AlternativeBoundary Values
Information LevelReferenceUpper
ProportionActualUpperAlpha
10.333366.666672.041243.00000
20.6667133.33332.886752.50000
31.00002003.535532.00000


The "Error Spending Information" in Output 103.5.12 displays cumulative error spending at each stage for each boundary. The first-stage spending 0.00135 corresponds to the one-sided p-value for a standardized Z statistic, .

Output 103.5.12: Error Spending Information

Error Spending Information
_Stage_Information
Level
Cumulative Error Spending
Upper
ProportionBetaAlpha
10.33330.000000.00135
20.66670.000000.00702
31.00000.060350.02532


With ODS Graphics enabled, a detailed boundary plot with the rejection and acceptance regions is displayed, as shown in Output 103.5.13. With the STOP=REJECT option, the interim rejection boundaries are displayed.

Output 103.5.13: Boundary Plot

Boundary Plot


With the PLOTS=ERRSPEND option, the procedure displays a plot of error spending for each boundary, as shown in Output 103.5.14. The error spending values in the "Error Spending Information" table in Output 103.5.10 are displayed in the plot.

Output 103.5.14: Error Spending Plot

Error Spending Plot