HD 189733 RM fit. Top: best fit model (solid red) of the combined observations of HD 189733b with HARPS. The gray area corresponds to a random sample of the posterior about the highest likelihood. The data points present two error bars, the green corresponding to the original RVs’ uncertainties computation from the CCF method and the black after quadratically adding the white-noise level fitted to the data. Bottom: residuals after subtraction of the data from the best-fit. The lighter yellow represents the ingress and egress and the darker the region where the planet is fully in front of the stellar disk. At the top of the residual’s plot, σ_res represents the average value of the RVs’ residual dispersion, σ_W is the estimate of the jitter velocity, and σ is the average uncertainty in the residual RVs.