R-Lum/Luminescence

New function - fit_SurfaceExposure()

Closed this issue · 1 comments

tzerk commented

As the name suggest, the new function fit_SurfaceExposure() is intended to fit eq. 1 in Sohbati et al. (2012a) to OSL surface exposure data.

Features

  • Fit 'standard' depth-intensity OSL surface exposure data using eq. 1 in Sohbati et al. (2012a)
  • (optional) Error weighted fitting (1/sigma^2)
  • Fixing main parameters (mu, sigma_phi, t(=age))
  • Global fit of multiple data sets to further constrain mu and/or sigma_phi
  • (optional) Include constant dose rate term (eq. 12 in Sohbati et al. 2012b)
  • (optional) consider variability of dose rate with sampling depth (D(x) in eq. 12) this feature may be included at a later stage

Tasks

  • Function
  • Documentation
  • Generate & include (multiple) synthetic data set(s)
  • testthat cases
  • Update NEWS

A shiny application is already available and will be included in the R package RLumShiny.

tzerk commented

A fully functional beta version is now implemented and documented. Incremental improvements are due after user feedback. So long, this issue can be closed.