The family (\(\gamma_{l}\)) and tree (\(\delta_{m}\)) effects were assumed to come from normal distributions \(N(0,\tau_{f})\) and \(N(0,\tau_{t})\), with vague gamma priors \(\tau_{f}\sim\Gamma(0.01,0.01)\) and \(\tau_{t}\sim\Gamma(0.01,0.01)\) respectively. The statistical model is also presented graphically following \cite{Kruschke2014} in Figure \ref{133147}.