data withdums;
set gandh.chapter3sim;
dum1=(FERTIL=1);
dum2=(FERTIL=2);
run;

data withparameters;
set withdums;
k3=13.5;
k1=1.2;
k2=3.4;
sigma=1.1;
run;

data;
set withparameters;
noise=normal(0);
y=k3+k1*dum1+k2*dum2+sigma*noise;
run;

proc glm;
class FERTIL;
model y=FERTIL / solution;
run;