The fully marginalized likelihood, or Bayesian evidence, is of great importance in probabilistic data analysis, because it is involved in calculating the posterior probability of a model or re-weighting a mixture of models conditioned on data. It is, however, extremely challenging to compute. This paper presents a geometric-path Monte Carlo method, inspired by multi-canonical Monte Carlo to evaluate the fully marginalized likelihood. We show that the algorithm is very fast and easy to implement and produces a justified uncertainty estimate on the fully marginalized likelihood. The algorithm performs efficiently on a trial problem and multi-companion model fitting for radial velocity data. For the trial problem, the algorithm returns the correct fully marginalized likelihood, and the estimated uncertainty is also consistent with the standard deviation of results from multiple runs. We apply the algorithm to the problem of fitting radial velocity data from HIP 88048 ($\nu$ Oph) and Gliese 581. We evaluate the fully marginalized likelihood of 1, 2, 3, and 4-companion models given data from HIP 88048 and various choices of prior distributions. We consider prior distributions with three different minimum radial velocity amplitude $K_{\mathrm{min}}$. Under all three priors, the 2-companion model has the largest marginalized likelihood, but the detailed values depend strongly on $K_{\mathrm{min}}$. We also evaluate the fully marginalized likelihood of 3, 4, 5, and 6-planet model given data from Gliese 581 and find that the fully marginalized likelihood of the 5-planet model is too close to that of the 6-planet model for us to confidently decide between them.

Markov Chain Monte Carlo (MCMC) proves to be powerful for Bayesian inference and in particular for exoplanet radial velocity fitting because MCMC provides more statistical information and makes better use of data than common approaches like chi-square fitting. However, the non-linear density functions encountered in these problems can make MCMC time-consuming. In this paper, we apply an ensemble sampler respecting affine invariance to orbital parameter extraction from radial velocity data. This new sampler has only one free parameter, and it does not require much tuning for good performance, which is important for automatization. The autocorrelation time of this sampler is approximately the same for all parameters and far smaller than Metropolis-Hastings, which means it requires many fewer function calls to produce the same number of independent samples. The affine-invariant sampler speeds up MCMC by hundreds of times compared with Metropolis-Hastings in the same computing situation. This novel sampler would be ideal for projects involving large datasets such as statistical investigations of planet distribution. The biggest obstacle to ensemble samplers is the existence of multiple local optima; we present a clustering technique to deal with local optima by clustering based on the likelihood of the walkers in the ensemble. We demonstrate the effectiveness of the sampler on real radial velocity data.

We report the detection of three new exoplanets from Keck Observatory. HD 163607 is a metal-rich G5IV star with two planets. The inner planet has an observed orbital period of 75.29 $\pm$ 0.02 days, a semi-amplitude of 51.1 $\pm$ 1.4 \ms, an eccentricity of 0.73 $\pm$ 0.02 and a derived minimum mass of \msini = 0.77 $\pm$ 0.02 \mjup. This is the largest eccentricity of any known planet in a multi-planet system. The argument of periastron passage is 78.7 $\pm$ 2.0$^{\circ}$; consequently, the planet's closest approach to its parent star is very near the line of sight, leading to a relatively high transit probability of 8%. The outer planet has an orbital period of 3.60 $\pm$ 0.02 years, an orbital eccentricity of 0.12 $\pm$ 0.06 and a semi-amplitude of 40.4 $\pm$ 1.3 \ms. The minimum mass is \msini = 2.29 $\pm$ 0.16 \mjup. HD 164509 is a metal-rich G5V star with a planet in an orbital period of 282.4 $\pm$ 3.8 days and an eccentricity of 0.26 $\pm$ 0.14. The semi-amplitude of 14.2 $\pm$ 2.7 \ms\ implies a minimum mass of 0.48 $\pm$ 0.09 \mjup. The radial velocities of HD 164509 also exhibit a residual linear trend of -5.1 $\pm$ 0.7 \ms\ per year, indicating the presence of an additional longer period companion in the system. Photometric observations demonstrate that HD 163607 and HD 164509 are constant in brightness to sub-millimag levels on their radial velocity periods. This provides strong support for planetary reflex motion as the cause of the radial velocity variations.