File: beir7err1.cc

package info (click to toggle)
simrisc 16.06.00-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 4,568 kB
  • sloc: cpp: 6,889; fortran: 669; makefile: 112; ansic: 112; sh: 107
file content (21 lines) | stat: -rw-r--r-- 696 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
//#define XERR
#include "screening.ih"

// called from (deprecated) errRisk1404.cc
// index is the biRad column index to use with BC
// and 0 when used with LC, as LC has no dose variation.

// static
double Screening::beir7Err(Round const &round, uint16_t index, double beta,
                           double eta, Modalities const &modalities)
{
    double risk = 1;
    double factor = beta * pow(round.age() / 60, eta) / 1000;

                            // for modalities not defining 'dose' (MRI)
                            // 0 is returned
    for (ModBase const *modBase: modalities.use(round.modalityIndices()))
        risk *= (1 + factor * modBase->dose(index));

    return risk;
}