calculate_sf {CFAcoop} R Documentation

## calculate_sf

### Description

calculates the survival fraction according to the procedure presented in Brix et al. (2020), which is robust against cellular cooperation.

### Usage

```calculate_sf(par_ref, par_treat, c_range = c(5, 20, 100))
```

### Arguments

 `par_ref` `summary.lm` object or 2-column matrix for the treatment-free reference survival `par_treat` `summary.lm` object or 2-column matrix for the clonogenic survival after treatment `c_range` colony numbers for which the survival fraction is calculated (default = c(5, 20, 100))

### Value

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)

### Examples

```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)
)
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)
```

[Package CFAcoop version 0.1.0 Index]