Let me give an answer i find very good and based on a small textbook on relativistic quantum mechanics (Relativistic Quantum Mechanics, S Trahanas, in greek)
First one can start with the known data of particles with integral spin (refered as bosons, obeying Bose-Einstein statistics) and particles with half-integral spin (refered as fermions obeying Fermi-Dirac statistics, or Pauli's exclusion principle)
One then can try to find a relativistic quantum equation (effectively a Schrödinger equation compatible with special relativity). One can note that the ordinary Schrödinger (free) equation
$$\hat{H}\psi=ih\frac{\partial \psi}{\partial t}$$
with a hamiltonian of the form $H=\hat{p}^2/2m$ reproduces the classical energy-momentum relation $E=p^2/2m$ which is not relativistic.
One then can use the same association of canonical variables into operators but in the relativistic energy-momentum relation $E^2=m^2 + p^2$ (1) (taking $c=1$), in order to derive a relativistic version of Schrödinger equation.
If one does that in relation (1), one gets the Klein-Gordon equation. This equation in good for bosons and the energy is bounded from below when bosons are described by this equation.
However if one tries to apply the KG equation to fermions, problems arise. For starters, fermions which obey Pauli's principle (thus have to be quantised with anti-commutators), make the KG equation have unbounded energy from below (which is something equivalent to a perpetuum mobile).
The problem is that the relation (1) is second order in the Energy (thus in equation form, is second order in time evolution $\partial/\partial t$).
Dirac tried to solve this by factoring the relation (1) into this form (2):
$$E=\mathbf{\alpha} m + \mathbf{\beta} p$$
but in order for the relation (2) to be equal to relation (1) when squared, $\mathbf{\alpha}$ and $\mathbf{\beta}$ are not ordinary numbers but some kind of matrices (specificaly pauli matrices). Thus we enter what is refered as formalism of spinors in the Dirac equation which is a quantum equation derived from relation (2).
Indeed the Dirac equation can describe fermions (aka Pauli exclusion principle, anti-commutators) and the energy is in fact bounded from below.
One artifact is that the dirac equation describes 2 particles and not one, since matrices/spinors are involved. Eventualy this led to the discovery of the positron (and anti-particles/anti-matter).
Finaly, the synthesis of quantum mechanics and relativity (along with anti-particles) led to the Quantum Field Theory formalism and the Spin-Statistics theorem, which relates theoreticaly the spin of a particle with the type of statistics it follows (integral spin->Bose-Einstein statistics, half-integral spin->Fermi-Dirac statistics)
References:
- PAM Dirac, The Quantum Theory of the Electron
- PAM Dirac, The Quantum Theory of the Electron. Part II
- PAM Dirac, A Positive-Energy Relativistic Wave Equation
- PAM Dirac, A Positive-Energy Relativistic Wave Equation. II