logLik.mkinfit: Calculated the log-likelihood of a fitted mkinfit object

In mkin: Kinetic Evaluation of Chemical Degradation Data

Calculated the log-likelihood of a fitted mkinfit object

Description

This function returns the product of the likelihood densities of each observed value, as calculated as part of the fitting procedure using dnorm, i.e. assuming normal distribution, and with the means predicted by the degradation model, and the standard deviations predicted by the error model.

Usage

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

Arguments

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

Details

The total number of estimated parameters returned with the value of the likelihood is calculated as the sum of fitted degradation model parameters and the fitted error model parameters.

Value

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

Author(s)

Johannes Ranke

See Also

Compare the AIC of columns of mmkin objects using AIC.mmkin.

Examples

  ## Not run: 
  sfo_sfo <- mkinmod(
    parent = mkinsub("SFO", to = "m1"),
    m1 = mkinsub("SFO")
  )
  d_t <- subset(FOCUS_2006_D, value != 0)
  f_nw <- mkinfit(sfo_sfo, d_t, quiet = TRUE) # no weighting (weights are unity)
  f_obs <- update(f_nw, error_model = "obs")
  f_tc <- update(f_nw, error_model = "tc")
  AIC(f_nw, f_obs, f_tc)
  
## End(Not run)

mkin documentation built on Nov. 23, 2023, 3:02 p.m.