calculate_sf {CFAcoop}  R Documentation 
calculates the survival fraction according to the procedure presented in Brix et al. (2020), which is robust against cellular cooperation.
calculate_sf(par_ref, par_treat, c_range = c(5, 20, 100))
par_ref 

par_treat 

c_range 
colony numbers for which the survival fraction is calculated (default = c(5, 20, 100)) 
vector of survival fractions.
If par_ref and par_treat are summary.lm
objects,
vector is of the same length as c_range.
If par_ref and par_treat are matrices,
vector is of the same length as nrow(par_treat)
seeded < 10^(seq(1, 5, 0.5)) counted.ref < 0.4 * 10^(seq(1, 5, 0.5) + rnorm(n = 9, 0, 0.1))^1.1 counted.treat < 0.01 * 10^(seq(1, 5, 0.5) + rnorm(n = 9, 0, 0.1))^1.2 fit_ref < pwr_reg(seeded = seeded, counted = counted.ref) fit_treat < pwr_reg(seeded = seeded, counted = counted.treat) calculate_sf(par_ref = fit_ref, par_treat = fit_treat) data("CFAdata") D < subset.data.frame( x = CFAdata, subset = cell.line == levels(CFAdata$cell.line)[1] ) fit_ref < pwr_reg(seeded = D$`Cells seeded`, counted = D$`0 Gy`) fit_treat < pwr_reg(seeded = D$`Cells seeded`, counted = D$`4 Gy`) calculate_sf(par_ref = fit_ref, par_treat = fit_treat)