51Peg
After installing kima, we'll go through one of the examples to start playing with the package.
Let's find 51 Peg b!
Redoing a first¶
51 Peg b was the first exoplanet discovered around a solar-type star. The original paper, by Mayor & Queloz (1995), used data from the ELODIE spectrograph. Here we will use another dataset, obtained with the Hamilton echelle spectrograph, situated in the Lick Observatory, in California (see Butler et al. 2006).
Let's import the package and the example
Note
The name _51Peg is a bit weird simply because Python variable names cannot start with a number. 🤷
To take a quick look at the RV data, let's call the _51Peg function to just build the model, without running it:
The data consist of 256 observations over almost 6 years, from October 1995 to October 2001. Our goal is to fit a Keplerian model to these radial-velocity observations. We will assumme the number of Keplerians to be free, with a uniform prior between 0 and 1, while all other priors take default values.
Running the model¶
We could run this model using the kima.run() function, but let's use the example directly (which
runs 5000 steps by default):
# Seeding random number generators. First seed = 1777119527.
# Generating 8 particles from the prior...done.
# Sampling...
[ ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0% ] 0%1%2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%# ] 2%3%4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%# ] 4%5%6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%# ] 6%7%8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%# ] 8%9%10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%# ] 10%11%12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%# ] 12%13%14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%# ] 14%15%16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%# ] 16%17%18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%# ] 18%19%20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%# ] 20%21%22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%# ] 22%23%24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%# ] 24%25%26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%# ] 26%27%28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%# ] 28%29%30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%# ] 30%31%32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%# ] 32%33%34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%# ] 34%35%36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%# ] 36%37%38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%# ] 38%39%40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%# ] 40%41%42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%# ] 42%43%44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%# ] 44%45%46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%# ] 46%47%48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%# ] 48%49%50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%# ] 50%51%52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%# ] 52%53%54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%# ] 54%55%56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%# ] 56%57%58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%# ] 58%59%60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%# ] 60%61%62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%# ] 62%63%64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%# ] 64%65%66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%# ] 66%67%68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%# ] 68%69%70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%# ] 70%71%72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%# ] 72%73%74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%# ] 74%75%76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%# ] 76%77%78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%# ] 78%79%80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%# ] 80%81%82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%# ] 82%83%84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%# ] 84%85%86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%# ] 86%87%88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%# ] 88%89%90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%# ] 90%91%92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%# ] 92%93%94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%# ] 94%95%96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%# ] 96%97%98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%# ] 98%99%100%#] 100%
# Took 28.274 seconds
Loading files took 0.08 seconds
log(Z) = -920.71
Information = 45.26 nats
BMD = 9.43
Effective sample size = 476.5
Analysing the results¶
We also loaded the results into the res variable, which we can use to look at some posterior distributions.
For example, the posterior for the number of planets \(N_p\) is quite clear, showing the significant detection of the planet:
Np probability ratios: []
The orbital period is also very well constrained:
Note
The plot above also shows samples from the prior distribution, which is log-uniform extending from 1 day to the timespan of the data.
We can also plot the histograms of the posteriors for the systemic velocity and for the instrumental jitter
[<Axes: title={'center': 'Posterior distribution for extra white noise'}>]
and the (more interesting) phase plot using the maximum likelihood solution
Sample with the highest likelihood value (logL = -869.50)
-> might not be representative of the full posterior distribution
jitter:
[2.95397022]
number of planets: 1
orbital parameters: P K M0 e w
4.23073 55.93585 4.58921 0.01323 1.10531
vsys: -1.77699292
So, kima tells us that a planet at a period of 4.23 days and an amplitude of ~56 m/s is the best model given the RV data. If it was 1995, we'd be going about changing the history of astronomy!