File: vsensitivity.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 (30 lines) | stat: -rw-r--r-- 805 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
22
23
24
25
26
27
28
29
30
//#define XERR
#include "mammo.ih"

//code
// cf. Isheden, G., & Humphreys, K. (2017). Modelling breast cancer tumour
// growth for a stable disease population. Statistical Methods in Medical
// Research, 28(3), 681-702.

double Mammo::vSensitivity(size_t idx) const
{
    double diameter = d_tumor.diameter();
    double mValue = d_m[idx].value();

    double expBeta =
        exp(
            d_beta[0].value()            +
            d_beta[1].value() * diameter +
            d_beta[2].value() * mValue   +
            d_beta[3].value() * mValue / (diameter * diameter)
        );

    double ret = (1 - d_sysErr) * expBeta / (1 + expBeta);

    g_log << __FILE__ " idx = " << idx << ", diameter: " << diameter <<
            ", m: " << mValue << ", sens.: " << ret << "\n";

    return ret;
}
//=