Fitting loglinear Bradley-Terry models (LLBT) for paired comparisons using the R package prefmod