Work with raw count data

In an ideal world, the people who count your samples also calculate the specific activity and associated error. In an imperfect world, these functions are here to calculate specific activity and error given the required information.

pb210_error_from_activity(activity, count_mass, count_time)

pb210_error_from_counts(counts, count_mass, count_time)

pb210_activity_from_counts(counts, count_mass, count_time)

pb210_counts_from_activity(activity, count_mass, count_time)

Arguments

activity	A specific activity of lead-210, in Bq / kg.
count_mass	The amount of mass each sample that was counted (kg). Can be a vector to specify different masses for each sample.
count_time	The time each sample was counted, as resolved by lubridate::as.duration(). Raw numbers are treated as a number of seconds, however we reccomend using lubridate::ddays() and/or lubridate::dhours() to make the units explicit.
counts	The number of measured decays

Value

A vector of specific activities or specific activity errors in Bq / kg.

Examples

pb210_error_from_activity(
  c(500, 100, 10, 5),
  count_mass = 0.5 / 1000,
  count_time = lubridate::ddays(1)
)
#> [1] 3.4020691 1.5214515 0.4811252 0.3402069

pb210_error_from_counts(
  c(21600, 4320, 432, 216),
  count_mass = 0.5 / 1000,
  count_time = lubridate::ddays(1)
)
#> [1] 3.4020691 1.5214515 0.4811252 0.3402069

pb210_activity_from_counts(
  c(21600, 4320, 432, 216),
  count_mass = 0.5 / 1000,
  count_time = lubridate::ddays(1)
)
#> [1] 500 100  10   5

pb210_counts_from_activity(
  c(500, 100, 10, 5),
  count_mass = 0.5 / 1000,
  count_time = lubridate::ddays(1)
)
#> [1] 21600  4320   432   216