Method based on 2001 manual by the International Atomic Energy Agency (IAEA). Cytogenetic Analysis for Radiation Dose Assessment, Technical Reports Series (2001). Retrieved from https://www.iaea.org/publications/6303/cytogenetic-analysis-for-radiation-dose-assessment.
Usage
estimate_whole_body_delta(
num_cases,
case_data,
fit_coeffs,
fit_var_cov_mat,
conf_int = 0.95,
protracted_g_value = 1,
aberr_module = c("dicentrics", "translocations", "micronuclei")
)Examples
#The fitting RDS result from the fitting module is needed. Alternatively, manual data
#frames that match the structure of the RDS can be used:
fit_coeffs <- data.frame(
estimate = c(0.001280319, 0.021038724, 0.063032534),
std.error = c(0.0004714055, 0.0051576170, 0.0040073856),
statistic = c(2.715961, 4.079156, 15.729091),
p.value = c(6.608367e-03, 4.519949e-05, 9.557291e-56),
row.names = c("coeff_C", "coeff_alpha", "coeff_beta")
)
fit_var_cov_mat <- data.frame(
coeff_C = c(2.222231e-07, -9.949044e-07, 4.379944e-07),
coeff_alpha = c(-9.949044e-07, 2.660101e-05, -1.510494e-05),
coeff_beta = c(4.379944e-07, -1.510494e-05, 1.605914e-05),
row.names = c("coeff_C", "coeff_alpha", "coeff_beta")
)
case_data <- data.frame(
ID= "example1",
N = 361,
X = 100,
C0 = 302,
C1 = 28,
C2 = 22,
C3 = 8,
C4 = 1,
C5 = 0,
y = 0.277,
y_err = 0.0368,
DI = 1.77,
u = 10.4
)
#FUNCTION ESTIMATE_WHOLE_BODY_DELTA
estimate_whole_body_delta(
num_cases = 1,
case_data = case_data,
fit_coeffs = as.matrix(fit_coeffs),
fit_var_cov_mat = as.matrix(fit_var_cov_mat),
conf_int = 0.95,
protracted_g_value = 1,
aberr_module = "dicentrics"
)
#> [[1]]
#> [[1]]$est_doses
#> lower estimate upper
#> 1 1.648182 1.931213 2.214245
#>
#> [[1]]$est_yield
#> lower estimate upper
#> 1 0.2048733 0.277 0.3491267
#>
#> [[1]]$AIC
#> [1] 7.057229
#>
#> [[1]]$conf_int
#> [1] 0.95
#>
#>