We investigate how the observed large-scale surface magnetic fields of low-mass stars (̃0.1-2 M☉), reconstructed through Zeeman-Doppler imaging, vary with age t, rotation and X-ray emission. Our sample consists of 104 magnetic maps of 73 stars, from accreting pre-main sequence to main-sequence objects (1 Myr ≲ t ≲ 10 Gyr). For non-accreting dwarfs we empirically find that the unsigned average large-scale surface field is related to age as t-0.655 ± 0.045. This relation has a similar dependence to that identified by Skumanich, used as the basis for gyrochronology. Likewise, our relation could be used as an age-dating method (`magnetochronology'). The trends with rotation we find for the large-scale stellar magnetism are consistent with the t...