The higher order (HO) correlation beyond the coupled-cluster single double (triple) CCSD(T) level of theory, second-order spin-orbit coupling (SOC), and core-valence (CV) correlation effects on bond length, r, vibrational frequency, ω, and dissociation energy, D, are studied for a set of 17 lanthanide containing diatomics, including lanthanum, europium, ytterbium, and lutetium oxides and halides. Convergence in the magnitudes of the SOC, CV, and HO corrections with respect to basis set size is examined using a sequence of double, triple, and quadruple-ζ basis sets, with the complete basis set (CBS) limit estimates provided in most cases. The CV effects on D, r, and ω are calculated to amount up to 1.3 kcal·mol, 0.008 Å, and 5 cm, respectively. A detailed analysis of the origin of the CV effect with a particular accounting for various subvalence shells reveals that, generally, the Ln 4d correlation makes a major contribution, although in some instances the lower-lying 4sp shells contribute largely and even more substantially than 4d. The second-order SOC effect evaluated via two-component and four-component relativistic techniques proves to be non-negligible, especially for heavier species, e.g., LaI, in which it is as large as 0.8 kcal·mol in D, 0.002 Å in r, and 1.3 cm in ω. Higher order correlation effects assessed through the CCSDT(Q) level are mostly less than 0.8 kcal·mol, 0.004 Å, and 5 cm; however, in species with prominence of nondynamical correlation, e.g., EuO, YbF, and LuO, the HO correction can amount to 1.2-1.6 kcal·mol, 0.005-0.008 Å, and 8-30 cm in D, r, and ω, respectively. In general, the [CCSD(T)+CV]/CBS + SOC + HO composite results are in good agreement with the available experimental data, exhibiting a mean absolute deviation of 1.8 kcal·mol in D, 0.0023 Å in r, and 3.5 cm in ω. A significant experimental outlier, the bond length in YbI, is revealed, implying the need for re-examination of the experimental data.