In this work, we continue large-scale ab initio computations for single ionized lanthanides. Extended atomic calculations for the set of ions from Pr II (Z = 59) to Gd II (Z = 64) have been performed in our previous work. In this study, ions from Tb II (Z = 65) to Yb II (Z = 70) are analyzed. By employing the same multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods that are implemented in the general-purpose relativistic atomic structure package GRASP2018, the energy levels and transition data of electric dipole (E1) transitions are computed. These computations are based on the strategies (with small variations) of Paper I. Accuracy of data is evaluated by comparing the computed energy levels with the data provided by the National Institute of Standards and Technology (NIST) database and with data from various methods. We obtain the average accuracy in the energy level compared with the NIST database: 6%, 5%, 4%, 5%, 3%, and 3% for Tb II, Dy II, Ho II, Er II, Tm II, and Yb II, respectively. We also provide extensive comparison of transition probabilities and wavelengths. Our results reach the average accuracy of transition wavelengths: 9%, 9%, 9%, 3%, 4%, and 11% for Tb II, Dy II, Ho II, Er II, Tm II, and Yb II, respectively.