Create degradation functions for known analytical solutions

Source: R/create_deg_func.R

Create degradation functions for known analytical solutions

create_deg_func(spec, use_of_ff = c("min", "max"))

Arguments

spec

List of model specifications as contained in mkinmod objects
use_of_ff

Minimum or maximum use of formation fractions

Value

Degradation function to be attached to mkinmod objects

Examples


SFO_SFO <- mkinmod(
  parent = mkinsub("SFO", "m1"),
  m1 = mkinsub("SFO"))

#> Successfully compiled differential equation model from auto-generated C code.
FOCUS_D <- subset(FOCUS_2006_D, value != 0) # to avoid warnings
fit_1 <- mkinfit(SFO_SFO, FOCUS_D, solution_type = "analytical", quiet = TRUE)

#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
fit_2 <- mkinfit(SFO_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE)

#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
# \dontrun{
if (require(rbenchmark))
  benchmark(
    analytical = mkinfit(SFO_SFO, FOCUS_D, solution_type = "analytical", quiet = TRUE),
    deSolve = mkinfit(SFO_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE),
    replications = 2)

#> Loading required package: rbenchmark
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#> Warning: Shapiro-Wilk test for standardized residuals: p =  0.0165
#>         test replications elapsed relative user.self sys.self user.child
#> 1 analytical            2   0.416    1.000     0.416        0          0
#> 2    deSolve            2   0.715    1.719     0.715        0          0
#>   sys.child
#> 1         0
#> 2         0
  DFOP_SFO <- mkinmod(
    parent = mkinsub("DFOP", "m1"),
    m1 = mkinsub("SFO"))

#> Successfully compiled differential equation model from auto-generated C code.
  benchmark(
    analytical = mkinfit(DFOP_SFO, FOCUS_D, solution_type = "analytical", quiet = TRUE),
    deSolve = mkinfit(DFOP_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE),
    replications = 2)

#>         test replications elapsed relative user.self sys.self user.child
#> 1 analytical            2   0.912     1.00     0.912    0.001          0
#> 2    deSolve            2   1.724     1.89     1.723    0.000          0
#>   sys.child
#> 1         0
#> 2         0
# }