*** Planos de Blocos Aleatorizados
*** Valentim Rodrigues Alferes (Universidade de Coimbra, 2002)
*** valferes@fpce.uc.pt
**
** Esta sintaxe distribui aleatoriamente as unidades (sujeitos)
** pelas condições ou tratamentos experimentais e pode ser utilizada
** com Planos Aleatorizados de Blocos (Simples ou Generalizados) e com
** Planos Completamente Aleatórios com n igual por tratmento [para
** Planos Completamente Aleatorizados com n igual ou diferente por
** tratamento, pode utilizar outra rotina neste mesmo site].
**
** Nos Planos de Blocos Aleatorizados (Simples), o número de unidades
** por bloco é igual ao número de tratamentos (p). O plano necessita 
** de w blocos de p unidades.
**
** Nos Planos Generalizados de Blocos Aleatorizados (ver exemplo), o 
** número de unidades por bloco deve ser múltiplo do número de trata-
** mentos (p). O plano requer w blocos de kp unidades (k?2).
**
** Nos Planos Completamente Aleatorizados com n igual por tratamento,
** deve introduzir número de blocos = 1.
**
** No exemplo (Planos Generalizados de Blocos Aleatorizados), temos:
** 18 sujeitos (unidades);
** 3 tratamentos (p);
** 2 blocos (w);
** 9 sujeitos por bloco (9=3p).

DATA LIST FREE /TRAT(A20).
* Introduza as designações dos tratamentos (A, B e C no exemplo).
BEGIN DATA
A B C 
END DATA.
COMPUTE TRATM=$CASENUM.
SAVE OUTFILE='%temp%\fout1'.

INPUT PROGRAM.
* Introduza o número de sujeitos (18 no exemplo).
LOOP N_CASOS=1 TO 18.

* Introduza o número de tratamentos (3 no exemplo).
COMPUTE N_TRATM = 3.

* Introduza o número de blocos (2 no exemplo).
COMPUTE N_BLOCOS = 2.
END CASE.
END LOOP.
END FILE.
END INPUT PROGRAM.
COMPUTE TRATM=MOD($CASENUM,N_TRATM)+1.
SORT CASES BY TRATM(A).
MATCH FILES/FILE=*/TABLE='%temp%\fout1'/BY TRATM.
COMPUTE BLOCO=MOD($CASENUM,N_BLOCOS)+1.
SORT CASES BY BLOCO(A) TRAT(A).
SET SEED RANDOM.
COMPUTE RANDOM=RV.UNIFORM(0,1).
SORT CASES BY BLOCO(A) RANDOM(A).
RANK VARIABLES=RANDOM(A) BY BLOCO.
COMPUTE ID=RRANDOM.
FORMATS ID(F8.0) BLOCO(F8.0).
LIST ID BLOCO TRAT.
SORT CASES BY TRAT(A) BLOCO(A) ID(A).
SPLIT FILE SEPARATE BY TRAT.
LIST BLOCO ID.