Nous présentons une nouvelle méthode de combinaison des mesures archéomagnétiques (inclinaison, déclinaison, intensité) qui repose sur une modélisation hiérarchique bayésienne. Cette modélisation prend en compte un niveau supplémentaire d'erreurs individuelles qui viennent se combiner aux erreurs expérimentales de mesure. Cette approche s'inscrit dans un cadre statistique plus général de combinaison des mesures et fournit ainsi une nouvelle méthode de calcul d'une moyenne scalaire ou directionnelle qui trouve de nombreuses applications en archéomagnétisme et en paléomagnétisme.