* Encoding: UTF-8.
*** Análise do poder: Duas médias independentes (T de Student)
*** Valentim Rodrigues Alferes (Universidade de Coimbra, 2003)
*** valferes@fpce.uc.pt

NEW FILE.
INPUT PROGRAM.
LOOP N=2 TO 1000.
END CASE.
END LOOP.
END FILE.
END INPUT PROGRAM.

* Introduza o poder pretendido (por defeito, PODER_P = 0.80).
COMPUTE PODER_P = 0.80.

* Introduza alpha (por defeito, ALFA = 0.05).
COMPUTE ALFA = 0.05.

* Introduza a diferença esperada entre as médias nos dois grupos
* (se não tem uma ideia aproximada, escreva DIF = -999 e preencha
* o campo do tamanho do efeito [D] de acordo com as convenções 
* de Cohen).
COMPUTE DIF = 2.00.

* Introduza o desvio-padrão (ponderado) nos dois grupos.
* (se não tem uma ideia aproximada, escreva DPP = -999 e preencha
* o campo do tamanho do efeito [D] de acordo com as convenções 
* de Cohen).
COMPUTE DPP = 3.00.

* Se escreveu -999 nos dois campos anteriores, introduza o tamanho
* do efeito (D) de acordo com as convenções de Cohen. Caso tenha
* preenchido os dois campos anteriores, não altere esta linha.
COMPUTE D = -999.

DO IF (DIF=-999 OR DPP=-999).
COMPUTE D_COHEN=D.
ELSE.
COMPUTE D_COHEN=DIF/DPP.
END IF.
COMPUTE GL=2*N-2.
COMPUTE NCP=ABS((D_COHEN*SQR(N))/SQR(2)).
COMPUTE T_ALPHA=IDF.T(1-ALFA/2,GL).
COMPUTE PODER1=1-NCDF.T(T_ALPHA,GL,NCP).
COMPUTE PODER2=1-NCDF.T(T_ALPHA,GL,-NCP).
COMPUTE PODER=PODER1+PODER2.
IF (PODER>=PODER_P) K=1.
COMPUTE W=ABS(PODER_P-PODER).
EXECUTE.
SORT CASES BY K(D) W(A).
IF ($CASENUM=1) N_TRAT=N.
IF ($CASENUM=1) N_TOTAL=2*N.
FORMATS DIF(F8.4) DPP(F8.4) D_COHEN(F8.4) N(F8.0) ALFA(F8.4)
PODER(F8.4) PODER_P(F8.4) N_TRAT(F8.0) N_TOTAL(F8.0).
LIST PODER_P ALFA DIF DPP D_COHEN/CASES=FROM 1 TO 1.
LIST N_TRAT N_TOTAL/CASES=FROM 1 TO 1.
SORT CASES BY W(A).
SEL IF($CASENUM<8).
SORT CASES BY N(A).
LIST N PODER.