Having just learned the definition of analytic functions, I found it surprising and somewhat counterintuitive that the set of functions over the complex plane whose real and imaginary parts each satisfy Laplace's Equation $\nabla^{2} = 0$ should have such incredible significance, since this is basically a set of two harmonic functions over $\mathbb{R}$. Does this extend to higher dimensions of hypercomplex numbers such as quaternions?
More formally, is there some generalization of the space of analytic functions to higher dimensions $ 2^n | n \in \mathbb{N}, n > 1$ where each component function individually satisfies Laplace's Equation? Do these functional spaces yield similar properties to that enumerated in the field of complex analysis, or are there any salient differences?