%estimate bresnahan/reiss model for transmission shops
global pin pinplus1 trans M
load Mazzeo.csv;
data=Mazzeo;
M=size(data,1);
for k=2:5
    d(:,k-1)=(data(:,6)>=k);
    s(:,k-1)=(data(:,6)+1>=k);
end

pin=[log(data(:,2)),ones(M,1),d];
pinplus1=[log(data(:,2)),ones(M,1),s];

trans=data(:,6);
hist(trans,[0:.5:30])

clear data Mazzeo d k s;

%beta=[.1;-.1;-.1;-.1;-.1;-.1];

beta=[1;-1;-1;-1;-1;-1];

beta=fminsearch('brlikelihood',beta);

diary ps1results
disp(['log population &  ' num2str(beta(1))]);
disp(['$\alpha_1$ &  ' num2str(beta(2))]);
disp(['$\alpha_2$ &  ' num2str(beta(3))]);
disp(['$\alpha_3$ &  ' num2str(beta(4))]);
disp(['$\alpha_4$ &  ' num2str(beta(5))]);
disp(['$\alpha_5$ &  ' num2str(beta(6))]);
disp(['Log-Likelihood & ' num2str(-brlikelihood(beta))])
disp(['# Obs.  ' num2str(M)])
diary off