High accuracy timing over packet networks is achieved by generating correction factors from multiple separation intervals and timing information contained in packets in both directions between a master and a slave. The methods are based on evaluating the weighted average of short-term, medium-term, and long-term measurements of local clock offset. Weighted averages are used to develop robust correction terms that are modified with an α-shaping factor to provide additional immunity to packet network instabilities.