Optimize algebraic systems which describe thermodynamic binding systems

Usage

sensitivity(
  case,
  parameters,
  path,
  additionalParameters,
  percentage = NULL,
  OffsetBoundaries = NULL,
  runAsShiny = FALSE
)

Arguments

case: is a character describing which system should be investigated. Either: "dba_host_const", "dba_dye_const", "ida" or "gda".
parameters: is a numeric vector containing already optimized parameter. In case of hg the order of the parameters is: khd, I0, IHD and ID In case of ida and ga the order of the parameters is: kg, I0, IHD and ID.
path: is a filepath which contains tabular x-y data. The concentraion of dye or guest respectivly is assumed to be in the first column. Furthermore, should the corresponding signal be stored in the second column. As an alternative an already loaded data.frame can be passed to the function.
additionalParameters: are required parameters which are specific for each case. In case of hg a numeric vector of length 1 is expected which contains the concentration of the host. In case of ida a numeric vector of length 3 is expected which contains the concentration of the host, dye and the khd parameter. In case of gda a numeric vector of length 3 is expected which contains the concentration of the host, guest and the khd parameter.
percentage: is the percentage +/- from parameters in which the sensitivity should be analysed.
OffsetBoundaries: in case percentage is not suitable a numeric vector (equivalent to parameters) can be used which is added/substracted from parameters. It is only possible to set either percentage or OffsetBoundaries.
runAsShiny: is internally used when running the algorithm from shiny.

Value

a plot showing the sensitivity

Examples

path <- paste0(system.file("examples", package = "tsf"), "/IDA.txt")
res <- opti("ida", c(1, 0, 0, 0), c(10^9, 10^6, 10^6, 10^6), path, c(4.3, 6.0, 7079458))
#> Warning: coercing argument of type 'double' to logical
#> Warning: coercing argument of type 'double' to logical
#> [1] ""
#> [1] 2
#> [1] "1.619e+08" "2.982e+00" "1.042e+03" "1.477e+02"
#> [1] "9.251e+00"
#> [1] ""
#> [1] 3
#> [1] "2.793e+08" "4.672e+00" "1.626e+03" "6.817e+01"
#> [1] "3.133e+00"
#> [1] ""
#> [1] 4
#> [1] "2.793e+08" "4.672e+00" "1.626e+03" "6.817e+01"
#> [1] "3.133e+00"
#> [1] ""
#> [1] 5
#> [1] "2.793e+08" "4.672e+00" "1.626e+03" "6.817e+01"
#> [1] "3.133e+00"
#> [1] ""
#> [1] 6
#> [1] "1.847e+08" "2.529e+02" "1.259e+03" "2.208e+01"
#> [1] "1.451e+00"
#> [1] ""
#> [1] 7
#> [1] "1.847e+08" "2.529e+02" "1.259e+03" "2.208e+01"
#> [1] "1.451e+00"
#> [1] ""
#> [1] 8
#> [1] "1.847e+08" "2.529e+02" "1.259e+03" "2.208e+01"
#> [1] "1.451e+00"
#> [1] ""
#> [1] 9
#> [1] "1.847e+08" "2.529e+02" "1.259e+03" "2.208e+01"
#> [1] "1.451e+00"
#> [1] ""
#> [1] 10
#> [1] "1.847e+08" "2.529e+02" "1.259e+03" "2.208e+01"
#> [1] "1.451e+00"
#> [1] ""
#> [1] 11
#> [1] "2.210e+08" "1.000e-15" "1.382e+03" "4.666e+01"
#> [1] "1.279e+00"
#> [1] ""
#> [1] 12
#> [1] "2.210e+08" "1.000e-15" "1.382e+03" "4.666e+01"
#> [1] "1.279e+00"
#> [1] ""
#> [1] 13
#> [1] "1.683e+08" "2.384e+02" "1.344e+03" "1.000e-15"
#> [1] "7.863e-01"
#> [1] ""
#> [1] 14
#> [1] "1.683e+08" "2.384e+02" "1.344e+03" "1.000e-15"
#> [1] "7.863e-01"
#> [1] ""
#> [1] 15
#> [1] "1.090e+08" "1.514e+02" "1.206e+03" "1.000e-15"
#> [1] "6.995e-01"
#> [1] ""
#> [1] 16
#> [1] "1.090e+08" "1.514e+02" "1.206e+03" "1.000e-15"
#> [1] "6.995e-01"
#> [1] ""
#> [1] 17
#> [1] "1.090e+08" "1.514e+02" "1.206e+03" "1.000e-15"
#> [1] "6.995e-01"
#> [1] ""
#> [1] 18
#> [1] "1.090e+08" "1.514e+02" "1.206e+03" "1.000e-15"
#> [1] "6.995e-01"
#> [1] ""
#> [1] 19
#> [1] "1.090e+08" "1.514e+02" "1.206e+03" "1.000e-15"
#> [1] "6.995e-01"
#> [1] ""
#> [1] 20
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 21
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 22
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 23
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 24
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 25
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 26
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 27
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 28
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 29
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 30
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 31
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 32
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 33
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 34
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 35
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 36
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 37
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 38
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 39
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 40
#> [1] "1.258e+08" "1.000e-15" "1.286e+03" "3.439e+01"
#> [1] "6.933e-01"
#> [1] ""
#> [1] 41
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 42
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 43
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 44
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 45
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 46
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 47
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 48
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 49
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 50
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 51
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 52
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 53
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 54
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 55
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 56
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 57
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 58
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 59
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 60
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 61
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 62
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 63
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 64
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 65
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 66
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 67
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 68
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 69
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 70
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 71
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 72
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 73
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 74
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 75
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 76
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 77
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 78
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 79
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 80
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 81
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 82
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 83
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 84
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 85
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 86
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 87
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 88
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 89
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 90
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 91
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 92
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 93
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 94
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 95
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 96
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 97
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 98
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 99
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 100
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 101
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 102
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 103
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 104
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 105
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 106
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 107
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 108
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 109
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 110
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 111
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 112
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 113
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 114
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 115
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 116
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 117
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 118
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 119
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 120
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 121
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 122
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 123
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 124
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 125
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 126
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 127
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 128
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 129
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 130
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 131
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 132
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 133
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 134
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 135
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 136
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 137
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 138
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 139
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 140
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 141
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 142
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 143
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 144
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 145
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 146
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 147
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 148
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 149
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 150
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 151
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 152
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 153
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 154
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 155
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 156
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 157
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 158
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 159
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 160
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 161
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 162
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 163
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 164
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 165
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 166
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 167
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 168
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 169
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 170
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 171
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 172
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 173
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 174
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 175
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 176
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 177
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 178
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 179
#> [1] "9.950e+07" "2.330e+01" "1.359e+03" "1.000e-15"
#> [1] "6.349e-01"
#> [1] ""
#> [1] 180
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 181
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 182
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 183
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 184
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 185
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 186
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 187
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 188
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 189
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 190
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 191
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 192
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 193
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 194
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 195
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 196
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 197
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 198
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 199
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
#> [1] ""
#> [1] 200
#> [1] "1.387e+08" "1.957e+02" "1.296e+03" "7.096e-01"
#> [1] "4.673e-01"
sensitivity("ida", res[[2]], path, c(4.3, 6.0, 7079458), 20)
#> [1] 1
#> [1] 2
#> [1] 3
#> [1] 4
#> [1] 5
#> [1] 6
#> [1] 7
#> [1] 8
#> [1] 9
#> [1] 10
#> [1] 11
#> [1] 12
#> [1] 13
#> [1] 14
#> [1] 15
#> [1] 16
#> [1] 17
#> [1] 18
#> [1] 19
#> [1] 20
#> [1] 21
#> [1] 22
#> [1] 23
#> [1] 24
#> [1] 25
#> [1] 26
#> [1] 27
#> [1] 28
#> [1] 29
#> [1] 30
#> [1] 31
#> [1] 32
#> [1] 33
#> [1] 34
#> [1] 35
#> [1] 36
#> [1] 37
#> [1] 38
#> [1] 39
#> [1] 40
#> [1] 41
#> [1] 42
#> [1] 43
#> [1] 44
#> [1] 45
#> [1] 46
#> [1] 47
#> [1] 48
#> [1] 49
#> [1] 50
#> [1] 51
#> [1] 52
#> [1] 53
#> [1] 54
#> [1] 55
#> [1] 56
#> [1] 57
#> [1] 58
#> [1] 59
#> [1] 60
#> [1] 61
#> [1] 62
#> [1] 63
#> [1] 64
#> [1] 65
#> [1] 66
#> [1] 67
#> [1] 68
#> [1] 69
#> [1] 70
#> [1] 71
#> [1] 72
#> [1] 73
#> [1] 74
#> [1] 75
#> [1] 76
#> [1] 77
#> [1] 78
#> [1] 79
#> [1] 80
#> [1] 81
#> [1] 82
#> [1] 83
#> [1] 84
#> [1] 85
#> [1] 86
#> [1] 87
#> [1] 88
#> [1] 89
#> [1] 90
#> [1] 91
#> [1] 92
#> [1] 93
#> [1] 94
#> [1] 95
#> [1] 96
#> [1] 97
#> [1] 98
#> [1] 99
#> [1] 100
#> Scale for x is already present.
#> Adding another scale for x, which will replace the existing scale.