BTm                     function to fit the Bradley-Terry model and
                        extensions
