A Bayesian approach to trends in cometary magnitude data