The quaternion normal distribution is derived and a number of characteristics are highlighted. The maximum likelihood estimation procedure in the quaternion case is examined and the conclusion is reached that the estimation procedure is simplified if the unknown parameters of the associated real probability density function are estimated. The quaternion estimator is then obtained by regarding these estimators as the components of the quaternion estimator. By means of an example attention is given to a test criterion which can be used in the quaternion model.