function f=brlikelihood(beta)

global pin pinplus1 M trans

L=(trans==0).*normcdf(-pinplus1*beta)+(trans>0&trans<5).*(normcdf(-pinplus1*beta)-normcdf(-pin*beta))+(trans>=5).*(ones(M,1)-normcdf(-pin*beta));

f=-sum(log(L));