Background: In luminescence dating

dose rate is the irradiation dose a sample absorbed per unit time. The estimation of a dose rate affects the reliability of the burial age. Quadratic Propagation of Uncertainty (QPU) is routinely used to assess the standard error of a dose rate estimate. However

dose rate calculation involves lots of non-linear transformations between various parameters. This complicates the application of QPU method in dose rate error estimation. Purpose: In the current study

a detailed introduction to the calculation of the annual dose rate in case of coarse quartz sediments is given. A Monte Carlo technique (a “parametric bootstrap” method) is employed to simulate the propagation of uncertainty in dose rate calculation. Methods: An open source R program used for performing the simulation is developed. A practical application of this technique is illustrated using a measured data set. Results & Conclusion: The Monte Carlo method is more flexible and simple in comparison with the QPU approach in dose rate error assessment. The stochastic scheme described in this article can be applied to any field of the analytical sciences. Abstract