A new method of analysing isothermal remanent magnetisation (IRM) acquisition curves based on cumulative log Gaussian analysis [Robertson and France, Phys. Earth Planet. Inter. 82 (1994) 223–234] is proposed. It is based on the curve fitting of the IRM acquisition curve versus the logarithm of the applied field with: (i) the acquisition curve on a linear scale, (ii) the acquisition curve expressed as a gradient, and (iii) the acquisition curve on a probability scale. Even when a sample is not saturated, its magnetic properties can be defined, although with less certainty. The number of magnetic components required for an optimal fit to a measured IRM acquisition curve is evaluated statistically. The method discriminates on the basis of different mineral coercivity. Therefore, additional rock-magnetic tests are still required to separate minerals with similar coercivities.