new file.
input Program.
compute #reps = 10.
compute #na = 5.
compute #nb = 16.




loop j = 1 to #nb.
compute b = rv.normal(0,1).
loop i = 1 to #na.
compute y = 6*i.
compute d = rv.normal(0,2).
loop k = 1 to #reps.
compute e = rv.normal(0,.5).
leave i, j, k, b, d , e, y.
end case.
end loop.

end loop.
end loop.
end file.
end input program.
execute.

compute y = y + b + d + e.
execute.


UNIANOVA
  y  BY i j
  /RANDOM = j
  /METHOD = SSTYPE(3)
  /INTERCEPT = INCLUDE
  /CRITERIA = ALPHA(.05)
  /EMMEANS = TABLES(i)  
  /DESIGN = i j i*j .


VARCOMP
  y  BY j   i
  /RANDOM = j
  /METHOD = MINQUE (1)
  /DESIGN
  /INTERCEPT = INCLUDE .

MIXED
  y  BY i j
  /CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1)
  SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
  PCONVERGE(0.000001, ABSOLUTE)
  /FIXED = i  | SSTYPE(3)
  /METHOD = REML
  /EMMEANS = TABLES(i)  
/test = "I" i 1 -1 0 0 0;
               i 1 0 -1 0 0;
               i 1 0 0 -1 0;
               i 1 0 0 0 -1
/test "L1" intercept 1 i 1 0 0 0 0
 /RANDOM j i*j  | COVTYPE(VC) .