% program to generate data for ch19ex1
% ch19ex2a.m
% the data is generated using the following model
% y=4.05 m +1.20d1+1.44d2+80
% x1=3.66m+1.16d1+1.27 d2+200
% x2= 5.92m+1.73d1+1.79d2+100

% first generate random numbers for m,d1 and d2
m=rand(20,1)+100;
d1=rand(20,1);
d2=rand(20,1);
yreg=4.05*(m-100)+1.20*d1+1.44*d2+90;
x1reg=3.66*(m-100)+1.16*d1+1.27*d2+200;
x2reg=5.92*(m-100)+1.73*d1+1.79*d2+100;
index=1:1:20;index=index';
zreg=[index m x1reg x2reg yreg]
% save ch19ex2_reg_data zreg -ascii
% Now create some random data for testing
m=rand(20,1)+100;
d1=rand(20,1);
d2=rand(20,1);
yreg=4.05*(m-100)+1.20*d1+1.44*d2+90;
x1reg=3.66*(m-100)+1.16*d1+1.27*d2+200;
x2reg=5.92*(m-100)+1.73*d1+1.79*d2+100;
index=1:1:20;index=index';
zreg=[index m x1reg x2reg yreg]
save ch19ex2_test_data zreg -ascii;