Boath activitiy and background should be arranged such that the
top of the core is first.
pb210_excess(activity, background = 0)
| activity | A vector of measured lead-210 specific activities (in Bq/kg) and
associated error. These can have |
|---|---|
| background | A vector of estimated background
lead-210 specific activity (in Bq/kg) and associated error.
These can have |
A vector with errors::errors() of the excess lead-210 specific activity.
Background is determined by the first point at which activity is less than
or equal to background. All excess values are set to NA below (after)
this point.
core <- pb210_simulate_core(depth_step = rep(1, 30)) %>% pb210_simulate_counting() pb210_excess( set_errors( core$activity_estimate, core$activity_se ), background = 10 )#> Errors: 3.639490 3.077746 2.648213 2.245967 1.906458 ... #> [1] 562.222222 399.212963 292.962963 207.916667 147.013889 102.754630 #> [7] 68.425926 43.379630 28.495370 16.574074 8.842593 2.361111 #> [13] NA NA NA NA NA NA #> [19] NA NA NA NA NA NA #> [25] NA NA NA NA NA NA