IF (region = 1) y = 5+ 2*gr83 +x+ rv.normal(0,2) .
IF (region = 2) y = rv.normal(110,2) .
IF (region = 3) y = 150  - 2*gr83 + rv.normal(0,2) .

EXECUTE .



UNIANOVA
  y  BY region  WITH gr83 x
  /METHOD = SSTYPE(3)
  /INTERCEPT = INCLUDE
  /CRITERIA = ALPHA(.05)
/outfile = design('c:\best.sav')
  /DESIGN = region gr83 gr83*region x region*x .