function [p1,p2,y,s2]=setup(file)



filename=[file '.dat'];

eval(['load ' filename]);



ab1=eval([file '(:,1)']);

h1=eval([file '(:,2)']);

ab2=eval([file '(:,3)']);

h2=eval([file '(:,4)']);



p1=h1./ab1;

p2=h2./ab2;



y=[log(p1./(1-p1)) log(p2./(1-p2))];

s2=[1./ab1./p1./(1-p1) 1./ab2./p2./(1-p2)];





