Description Usage Arguments Details Value References Examples

Calculating the annual dose rate and burial age.

1 2 3 4 |

`dose` |
vector( |

`minGrainSize` |
numeric( |

`maxGrainSize` |
numeric( |

`Ucontent` |
vector( |

`Thcontent` |
vector( |

`Kcontent` |
vector( |

`Wct` |
vector( |

`depth` |
numeric( |

`altitude` |
numeric( |

`latitude` |
numeric( |

`longitude` |
numeric( |

`bulkDensity` |
numeric(with default): average density of bulk sample (unit, g/cm^3) |

`alphaValue` |
numeric(with default): average alpha efficiency |

`nsim` |
integer(with default): number of Monte Carlo iterations |

`rdcf` |
numeric(with default): constant relative standard error for dose-rate conversion factors (unit, 1 percent) |

`rba` |
numeric(with default): constant relative standard error for beta dose absorption fraction (unit, 1 percent) |

`plot` |
logical(with default): draw a plot or not |

Function calDA is used for calculating the annual dose rate and burial age using concentrations of radioactive nuclides (uranium, thorium, potassium) obtained from Neutron Activation Analysis (NAA), grain size, water content, average sample density, geographical parameters (depth, altitude, latitude, longitude), and an equivalent dose value. The elemental concentrations are converted into annual dose rate according to dose-rate conversion factors (Adamiec and Aitken, 1998). Beta dose absorded fractions are calculated through linear interpolation using published data (Fain et al., 1999). The cosmic ray dose rate is estimated as a function of depth, altitude and geomagnetic latitude (Prescott and Hutton, 1994).

The standard error of an annual dose rate is estimated by a "parametric bootstrap" method. To perform the simulation, dose-rate conversion factors and beta dose absorption factor are assumed to have constant relative standard errors. Arguments such as `dose`

, `Ucontent`

, `Thcontent`

, `Kcontent`

, `wct`

should be two-element vectors with the form `c(value,sd)`

, as their measurement uncertainties are major error sources. Arguments such as `depth`

, `altitude`

, `latitude`

, `longitude`

, `bulkDensity`

, `alphaValue`

can be either a scalar of a two-element vector. This means that uncertainties from these quantities can be either ignored or taken into account during the simulation.

Return a matrix that contains calculated annual dose rate and age

Adamiec G, Aitken M, 1998. Dose-rate conversion factors: update. Ancient TL, 16(2): 37-49.

Fain J, Soumana S, Montret M, Miallier D, Pilleyre T, Sanzelle S, 1999. Luminescence and ESR dating Beta-dose attenuation for various grain shapes calculated by a Monte-Carlo method. Quaternary Science Reviews, 18(2): 231-234.

Prescott JR, Hutton JT, 1994. Cosmic ray contributions to dose rates for Luminescence and Esr dating: large depths and long-term time variations. Radiation Measurements, 23(2-3): 497-500.

1 2 3 4 5 6 |

```
Pars stdPars
DoseRate 3.314279 0.2266184
Age 7.555188 0.5504003
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.