new file.
input Program.
compute #n = 20.
loop j = 1 to 4.
loop i = 1 to #n.
compute x = j.
leave i, j.
end case.
end loop.
end loop.
end file.
end input program.
execute.




IF (j=1) x = rv.uniform(1,5) .
IF (j=2) x = rv.uniform(6,8) .
IF (j=3) x = rv.uniform(10,15) .
IF (j=4) x = rv.uniform(16,21) .

EXECUTE .



IF (j=1) y = 3*exp(-(x-1)) + rv.normal(0,.1).
IF (j=2) y = 2*exp(-1.5*(x-6)) + rv.normal(0,.1).
IF (j=3) y = 4*exp(-2.25*(x-10)) + rv.normal(0,.1).
IF (j=4) y = 6*exp(-2.75*(x-16)) + rv.normal(0,.1).
EXECUTE .
GRAPH
  /SCATTERPLOT(BIVAR)=x WITH y BY j
  /MISSING=LISTWISE .