We present a method for the analysis of radial velocity (RV) measurements, in the context of searching for planets around chromospherically active stars. We assume that the observed RV signal is induced by the reflex motion of the star as well as by distortions of spectral line profiles, measured by the Bisector Velocity Span (BVS). The RVs are fitted with a Keplerian planetary model complemented by a correction term that is linearly dependent on the BVS. The BVS represents the stellar activity contribution to the observed RV. The coefficient of this linear dependence is an additional free parameter of the model. This approach differs from methods in the literature, which make usage of the BVS measurements before or after fitting the planetary model. We test the method on a simulated data set of a 1-planet configuration. The results are compared with the outcomes of algorithms found in the literature.