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