/******************************************************************************
 * Programme : Puissance Statistique : Ne laissez plus le hasard décider de la taille de vos échantillons (T-test & GLM)
 * Reference : ANALYS8217
 * Source    : https://www.wearecas.eu/en/sampleCode/ANALYS8217
 ******************************************************************************/

/* --- BLOC 1 --- */
data sleep;
  do hours = 1 to 2;
    input recall;
    output;
  end;
  datalines;
30
33
;
run;

/* --- BLOC 2 --- */
proc glmpower data=sleep;
  class hours;
  model recall = hours;
  power
    stddev = 30
    ntotal = .
    power = 0.7 0.8 0.9;
  plot x=power min=0.7 max=0.95;
run;

/* --- BLOC 3 --- */
proc power;
  twosamplemeans
    meandiff=3
    alpha=0.05
    ntotal = .
    power = 0.7 0.8 0.9
    stddev = 30
  ;
run;