We study the Kaiser formula for biased tracers in massive neutrino cosmologies by comparing its predictions with a large set of N-body simulations. In particular, we examine the ambiguity in the definition of the peculiar velocity contribution at linear level, whether it should be expressed in terms of the total matter velocity field or only the cold-matter component, comprised of cold dark matter and baryons. We revisit and extend previous qualitative studies on the topic with larger statistics and including a full fit of halo power spectrum and halo-matter cross-power spectrum measurements. We find clear evidence that the clustering of halos in redshift space is correctly described assuming that halo velocity is tracing the velocity field of cold-matter. The opposite assumption provides a worse fit to the simulated data and can lead to a spurious running of nuisance parameters in the perturbative galaxy power spectrum model.
