Analytical four-loop result for the β-function in QED in MS and MOM schemes