Monte Carlo algorithm. #"Add"a"planet"at"the"period"corresponding #"to"the"highest"peak kadd.planet(k,,c(period,=,p[1,,'period'])) kminimize(k) plot(k) #"Creates"a"new"model"object k,<:,knew() #"Load"new"data kadd.data(k,,"1pl.vels") #"Calculate"the"power"spectrum"of"the"data p,<:,kperiodogram(k) 0 1 2 3 4 -50 0 50 Time [JD] Radial velocity [m/s] 0 1 2 3 4 -40 -20 0 20 40 Time [JD] Residuals [m/s] 1 10 100 1000 10000 0 20 40 60 80 Period [d] Normalized power #"Run"a"Markov>Chain"Monte"Carlo"analysis #"(with"default"parameters) kmcmc(k,,chains=5)