Calculated the log-likelihood of a fitted mkinfit object

This function simply calculates the product of the likelihood densities calc

# S3 method for mkinfit
logLik(object, ...)

Arguments

object	An object of class mkinfit.
…	For compatibility with the generic method

Value

An object of class logLik with the number of estimated parameters (degradation model parameters plus variance model parameters) as attribute.

Examples

  sfo_sfo <- mkinmod(
    parent = mkinsub("SFO", to = "m1"),
    m1 = mkinsub("SFO")
  )
#> Successfully compiled differential equation model from auto-generated C code.
  d_t <- FOCUS_2006_D
  d_t[23:24, "value"] <- c(NA, NA) # can't cope with zero values at the moment
  f_nw <- mkinfit(sfo_sfo, d_t, quiet = TRUE) # no weighting (weights are unity)
  f_obs <- mkinfit(sfo_sfo, d_t, reweight.method = "obs", quiet = TRUE)
  f_tc <- mkinfit(sfo_sfo, d_t, reweight.method = "tc", quiet = TRUE)
  d_t$err <- d_t$value # Manual weighting assuming sigma ~ y
  f_man <- mkinfit(sfo_sfo, d_t, err = "err", quiet = TRUE)
  AIC(f_nw, f_obs, f_tc, f_man)
#>       df      AIC
#> f_nw   5 204.4619
#> f_obs  6 205.8727
#> f_tc   6 143.8773
#> f_man  4 291.8000