This work describes a methodology of hypothesis testing for possible existence of a body orbiting a pulsar. The methodology based on Bayesian inference. The tested hypotheses are a single pulsar, that exhibits power-law red noise in timing residuals, and a binary pulsar with white measurement noise. The methodology was tested on a series of simulated data and was used for analysis of pulsar’s B0329+54 timing residuals.