Bayesian deviance information criterion (BIC) was used to estimate the goodness of fit of consecutive models [31] The BIC values for each model were compared, and the model with the lowest value was considered the better one for hierarchical models [31].