hold on
scatter(beta2015,gamma2015,Markersize,Markerstyle,'MarkerFaceColor',[211,211,211]./255,'MarkerEdgeColor',[211,211,211]./255)
plot(beta_line,gamma_line,'r-','LineWidth',2)