Update belkin_schofer_redo.R

plot.hazard simulates coefficient values to get uncertainty around hazard