P. C. Gregory 2005 ApJ 631 1198 doi:10.1086/432594
P. C. Gregory1
Show affiliationsA new Bayesian tool for nonlinear model fitting has been developed that employs a parallel tempering Markov chain Monte Carlo algorithm with a novel statistical control system. The algorithm has been used to reanalyze precision radial velocity data for HD 73526. For a single-planet model, three possible orbits were found with periods of 127.88
, 190.4
, and 376.2
days. The 128 day orbit, with an eccentricity of 0.56 ± 0.03, has a maximum value of prior × likelihood that is 16 times larger than for the next highest solution at 376 days. However, the 376 day orbit, with an eccentricity of 0.10
, is formally more probable because for this sparse data set there is a much larger volume of parameter space with a significant probability density in the vicinity of the 376 day peak. The previously reported orbit (Tinney et al. 2003) of 190.5 ± 3.0 days corresponds to our least probable orbit. The analysis highlights the need for measurements around phase 0.5 for the 376 day period.
methods: data analysis; planetary systems; stars: individual (HD 73526)
Issue 2 (2005 October 1)
Received 2004 December 22, accepted for publication 2005 June 7
P. C. Gregory 2005 ApJ 631 1198
A. Finoguenov and C. Jones 2001 ApJ 547 L107
Anatoly A Gitelson et al 2009 Environ. Res. Lett. 4 045003
Christian Reisswig et al 2007 Class. Quantum Grav. 24 S327
O. Godet et al 2009 ApJ 705 L109
T. Bandyopadhyay 2008 EPL 81 16003
Erwin T. Lau et al. 2009 ApJ 705 1129
Shilong Piao et al 2009 Environ. Res. Lett. 4 044012
Shi-Jie Xiong and Ye Xiong 2007 EPL 80 60008
Joseph J Fins, MD, FACP 2009 J. Neural Eng. 6 050201