B-11 NMR chemical shifts of tricoordinated BX3 and tetracoordinated BX4- compounds (X = H, CH3, F, Cl, Br, I, OH, SH, NH2, and CH=CH2) were computed, and the shielding tensors were explored not only within the non relativistic GIAO approach but also by application of both relativistic ZORA computations including spin orbit coupling as well as scalar nonrelativistic ZORA computations (BP86 level of density functional theory). The contributions of the spin orbit coupling to the overall shieldings are decisive for X = Br and I in both series. No relationship was found between the 2p orbital occupancies or 1/Delta E (difference between LUMO and suitably occupied MO that can be coupled with LUMO) with the shielding tensors (or their principal values) in the BX3 series. However, a multidimensional statistical approach known as factor analysis (frequently used in chemometrics) revealed that three factors account for 92% of the cumulative proportion of total variance. The main components of the first factor are occupancies in the 2p(x) and 2p(y)), orbitals and 1/Delta E; the second factor is mainly the occupancy in the 2p(z) orbital and the inductive substituent parameters by Taft. Finally, the third factor consists exclusively (98.4%) of the electrostatic potential (V-max), which is directly related to the so-called pi-hole magnitudes.