load lalonde_exper

res=[];

% 1. 

Sample=1;
All=1; 
M=1;             
BiasAdj=0;  
Weight=1;   
Var_calc=0;  
X=[age];             
Z=X;             
V=educ;
Y=re78/1000;
T=treat;
K=1;
Weight_matrix=eye(K);             
exact=0;
n=length(Y);
weight=ones(n,1);

w_educ=exp((educ-10.1)/2);


[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[1,est,se];
res=[res;rr];

X=[age,educ,re74,re75];             
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[2,est,se];
res=[res;rr];
rr=[3,est,se];
res=[res;rr];

M=3;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[4,est,se];
res=[res;rr];
M=1;

All=0;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[5,est,se];
res=[res;rr];
All=1;

All=2;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[6,est,se];
res=[res;rr];
All=1;

Weight=2;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[7,est,se]
res=[res;rr];
Weight=1;

Weight_matrix=eye(4);             
Weight=3;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[8,est,se];
res=[res;rr];
Weight=1;

BiasAdj=1;
Z=X;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[9,est,se];
res=[res;rr];
BiasAdj=0;

BiasAdj=1;
Z=[married age];
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[10,est,se];
res=[res;rr];
BiasAdj=0;

exact=1;
V=age;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[11,est,se];
res=[res;rr];
exact=0;

exact=1;
V=[married u74];
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[12,est,se];
res=[res;rr];
exact=0;

Sample=0;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[13,est,se];
res=[res;rr];
Sample=1;

Var_calc=3;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[14,est,se];
res=[res;rr];
Var_calc=0;

weight=w_educ;
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[15,est,se];
res=[res;rr];
weight=ones(length(Y),1);
%pause



Sample=0;
M=3;
All=0;
V=[age];
BiasAdj=1;
Z=[married age];
X=[age educ re74 re75];
Weight=3;
exact=1;
Var_calc=3;
weight=w_educ;
Weight_matrix=eye(4);             
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[75,est,se]
res=[res;rr];


Sample=1;
M=3;
All=1;
V=[married];
BiasAdj=1;
Z=[age re75];
X=[age educ re74];
Weight=2;
exact=1;
Var_calc=0;
weight=w_educ;
Weight_matrix=eye(4);             
[est,se,se_cond,em,art_data,aug_data]=match(Y,T,X,Z,V,All,M,BiasAdj,Weight,Weight_matrix,Var_calc,weight,Sample,exact);
rr=[76,est,se]
res=[res;rr];

diary lalonde_exper_04feb2.dia
res
diary off