|
** 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=OUTF1.
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=OUTF1/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.
|