  We collected the archival data of blazar OJ~287 from heterogeneous very long
baseline interferometry (VLBI) monitoring programs at 2.3 GHz, 8.6 GHz, 15 GHz
and 43 GHz. The data reduction and observable extraction of those multi-band
multi-epoch observations are batch-processed consistently with our automated
pipeline. We present the multivariate correlation analysis on the observables
at each band. We employ the cross-correlation function to search the
correlations and the Monte Carlo (MC) technique to verify the certainty of
correlations. Several correlations are found. The foremost findings are the
correlations between the core flux density and the jet position angles on
different scales, which validated the plausible predictions of the jet with
precession characteristics. Meanwhile, there is a variation in the offset
between the core EVPA and the inner-jet position angle over time at 15~GHz and
43~GHz.
