Markov Chain Monte Carlo Methods Applied to Photometric Spot Modeling

I demonstrate that Markov chain Monte Carlo methods can be used very effectively to determine best‐fit values, uncertainties, and possible correlations or degeneracies in the fitted parameters of photometric spot modeling. Details of the Markov chain Monte Carlo methods applied here are briefly described, including the Metropolis‐Hastings algorithm and the tests that are used to ensure proper convergence and mixing. This Markov chain Monte Carlo functionality is applied to recent observations of ϵ Eridani by the Microvariablity and Oscillations of Stars (MOST) satellite, and the two‐spot solution showing differential rotation, as discussed in B. Croll et al. Conclusions in the latter are largely confirmed, but a strong correlation between the inclination and other fitted parameters is noted. The Markov chain Monte Carlo functionality has been included in StarSpotz, a freely available program for photometric spot modeling.

[1]  S. Ravi.  Bayesian Logical Data Analysis for the Physical Sciences: a Comparative Approach with Mathematica® Support , 2007 .

[2]  H. Dole,et al.  The 1–1000 μm spectral energy distributions of far-infrared galaxies , 2006 .

[3]  P. Gregory.  A Bayesian Analysis of Extrasolar Planet Data for HD 73526 , 2005 .

[4]  J. Valenti,et al.  The Planet-Metallicity Correlation , 2005 .

[5]  E. I. Robson,et al.  Structure in the ϵ Eridani Debris Disk , 2005 .

[6]  E. Ford.  Quantifying the Uncertainty in the Orbits of Extrasolar Planets , 2003, astro-ph/0305441.

[7]  K. Strassmeier,et al.  Time‐series photometric spot modelling – VI. A new computer code and its application to 23 years of photometry of the active giant IM Pegasi , 2003 .

[8]  Edward J. Wollack,et al.  First-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Parameter Estimation Methodology , 2003, astro-ph/0302218.

[9]  N. Christensen,et al.  Bayesian methods for cosmological parameter estimation from cosmic microwave background measurements , 2000, astro-ph/0103134.

[10]  Mullard Space Science Laboratory,et al.  A Dust Ring around epsilon Eridani: Analog to the Young Solar System , 1998, astro-ph/9808224.

[11]  Dani Gamerman,et al.  Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference , 1997 .

[12]  A. Gelman,et al.  Weak convergence and optimal scaling of random walk Metropolis algorithms , 1997 .

[13]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[14]  William H. Press,et al.  Numerical recipes in FORTRAN (2nd ed.): the art of scientific computing , 1992 .

[15]  J. D. Dorren,et al.  A New Formulation of the Starspot Model, and the Consequences of Starspot Structure , 1987 .

[16]  E. Budding.  The interpretation of cyclical photometric variations in certain dwarf Me-type stars , 1977 .

[17]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[18]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[19]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[20]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[21]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .