  We report our progress in data analysis on two-point correlation functions of
the $B$ meson using sequential Bayesian method. The data set of measurement is
obtained using the Oktay-Kronfeld (OK) action for the bottom quarks (valence
quarks) and the HISQ action for the light quarks on the MILC HISQ lattices. We
find that the old initial guess for the $\chi^2$ minimizer in the fitting code
is poor enough to slow down the analysis somewhat. In order to find a better
initial guess, we adopt the Newton method. We find that the Newton method
provides a natural test to check whether the $\chi^2$ minimizer finds a local
minimum or the global minimum, and it also reduces the number of iterations
dramatically.
