data = xlsread('wage2.csv');

brthord = data(:,3);

data = data(brthord>0,:);

lwage = data(:,1);
brthord = data(:,3);
edu = data(:,2);
sibs = data(:,4);

n = length(lwage);

Y = lwage;
X = [ones(n,1) edu];

disp('OLS Estimation')
rr1 = ols(Y,X,[],1);

W = [ones(n,1) sibs];
disp('IV Estimation')
rr2 = iv(Y,X,W,[],1);

W = [ones(n,1) brthord sibs];
disp('TSLS Estimation')
rr3 = tsls(Y,X,W,[],1);
