In this paper, we estimate the Stark and Zeeman shifts in the transition frequencies of the 16 O2 + molecular ion, as a step for the search for the variation in the proton-to-electron mass ratio µ. The X 2 Π v = 21 − a 4 Π v = 0 or the X 2 Πv = 21 − a 4 Π v = 1 transition frequencies (THz region) of the 16 O2 + molecular ion have particularly high sensitivity to the variation in µ. Note also that the Stark shift in the 16 O2 + transition frequencies is expected to be much smaller than that for heteronuclear diatomic molecules. However, the actual systematic uncertainties for the 16 O2 + transition frequencies have never been estimated. We estimated the Stark and Zeeman shifts in the different 16 O2 + transition frequencies. When the molecular ions in a string crystal formed in a linear trap (trap electric field < 0.1 V/cm, and Stark shift < 10 -20 ) are used, the X 2 Π 1/2 (v, J) = (0, 1/2) − (v ′ , 1/2)(v ′ ≥ 1) transition frequencies are most advantageous for the search for the variation in µ ∆µ/µ < 10 −17 because the Zeeman shift is easily suppressed to lower than 10 -18 and the electric quadrupole shift is zero. On the other hand, the X 2 Π 1/2 (v, J) = (21, 1/2) − a 4 Π 1/2 (v, J) = (0, 1/2) transition frequency has another merit that the positive Stark shift induced by the trap electric field can be canceled by the quadratic Doppler shift. Therefore, the measurement using molecular ions in a Coulomb crystal broadened in the radial direction is also possible, when the Zeeman shift is effectively eliminated. PACS numbers: 32.60.+i, 06.20.Dk, 33.20Ea