I am trying to understand how complex numbers made their way into QM. Can we have a theory of the same physics without complex numbers? If so, is the theory using complex numbers easier?
The nature of complex numbers in QM turned up in a recent discussion, and I got called a stupid hack for questioning their relevance. Mainly for therapeutic reasons, I wrote up my take on the issue:

On the Role of Complex Numbers in Quantum Mechanics

Motivation

It has been claimed that one of the defining characteristics that separate the quantum world from the classical one is the use of complex numbers. It's dogma, and there's some truth to it, but it's not the whole story:

While complex numbers necessarily turn up as first-class citizen of the quantum world, I'll argue that our old friend the reals shouldn't be underestimated.

A bird's eye view of quantum mechanics

In the algebraic formulation, we have a set of observables of a quantum system that comes with the structure of a real vector space. The states of our system can be realized as normalized positive (thus necessarily real) linear functionals on that space.

In the wave-function formulation, the Schrödinger equation is manifestly complex and acts on complex-valued functions. However, it is written in terms of ordinary partial derivatives of real variables and separates into two coupled real equations - the continuity equation for the probability amplitude and a Hamilton-Jacobi-type equation for the phase angle.

The manifestly real model of 2-state quantum systems is well known.

Complex and Real Algebraic Formulation

Let's take a look at how we end up with complex numbers in the algebraic formulation:

We complexify the space of observables and make it into a $C^*$-algebra. We then go ahead and represent it by linear operators on a complex Hilbert space (GNS construction).

Pure states end up as complex rays, mixed ones as density operators.

However, that's not the only way to do it:

We can let the real space be real and endow it with the structure of a Lie-Jordan-Algebra. We then go ahead and represent it by linear operators on a real Hilbert space (Hilbert-Schmidt construction).

Both pure and mixed states will end up as real rays. While the pure ones are necessarily unique, the mixed ones in general are not.

The Reason for Complexity

Even in manifestly real formulations, the complex structure is still there, but in disguise:

There's a 2-out-of-3 property connecting the unitary group $U(n)$ with the orthogonal group $O(2n)$, the symplectic group $Sp(2n,\mathbb R)$ and the complex general linear group $GL(n,\mathbb C)$: If two of the last three are present and compatible, you'll get the third one for free.

An example for this is the Lie-bracket and Jordan product: Together with a compatibility condition, these are enough to reconstruct the associative product of the $C^*$-algebra.

Another instance of this is the Kähler structure of the projective complex Hilbert space taken as a real manifold, which is what you end up with when you remove the gauge freedom from your representation of pure states:

It comes with a symplectic product which specifies the dynamics via Hamiltonian vector fields, and a Riemannian metric that gives you probabilities. Make them compatible and you'll get an implicitly-defined almost-complex structure.

Quantum mechanics is unitary, with the symplectic structure being responsible for the dynamics, the orthogonal structure being responsible for probabilities and the complex structure connecting these two. It can be realized on both real and complex spaces in reasonably natural ways, but all structure is necessarily present, even if not manifestly so.

Conclusion

Is the preference for complex spaces just a historical accident? Not really. The complex formulation is a simplification as structure gets pushed down into the scalars of our theory, and there's a certain elegance to unifying two real structures into a single complex one.

On the other hand, one could argue that it doesn't make sense to mix structures responsible for distinct features of our theory (dynamics and probabilities), or that introducing un-observables to our algebra is a design smell as preferably we should only use interior operations.

While we'll probably keep doing quantum mechanics in terms of complex realizations, one should keep in mind that the theory can be made manifestly real. This fact shouldn't really surprise anyone who has taken the bird's eye view instead of just looking throught the blinders of specific formalisms.
