One can find in the literature three main sets of estimators for the variance components in the hierarchical credibility model. This paper presents these estimators in a unified notation, studies some of their properties important for numerical evaluation and compares their relative performance by simulation. The paper also demonstrates how function cm of the R package actuar can be used to fit hierarchical models to insurance data.