KatさんがMCMC法でEAのタイプの周期を求めることをすすめて下さり、この間ブログに書いたようにやってみました。 MCMC法でKaiV20の要素を改良 - muttenz's blog そのスクリプトの fph <- 0.5 - abs((((x-p[2])/p[1]) %% 1) - 0.5) yxep <- ifelse(fph < p[3]…