https://github.com/cran/qpcR
Raw File
Tip revision: 50b4350d7c8be8e5917ae3d0580ab19842dacf4e authored by Andrej-Nikolai Spiess on 08 August 1977, 00:00:00 UTC
version 1.1-1
Tip revision: 50b4350
evidence.Rd
\name{evidence}
\alias{evidence}

\title{Evidence ratio for model comparison with AIC or AICc}

\description{
 The evidence ratio \deqn{\frac{1}{exp(-0.5 * (AIC2 - AIC1))}} is calculated for fitted models either from two 'drc' objects or two numerical values.
 Models can be compared that are not nested and where the f-test on residual-sum-of-squares is not applicable.
}

\usage{
  evidence(x, y, type = c("AICc", "AIC"))
}

\arguments{
  \item{x}{a 'drc' object or numerical value.}
  \item{y}{a 'drc' object or numerical value.}
  \item{type}{Bias-corrected (\code{\link{AICc}}, default) or original (\code{\link{AIC}}) versions of the Akaike Information Criterion.} 
}

\details{
Small differences in AIC values can mean substantial more 'likelihood' of one model over the other.
For example, a model with AIC = -130 is nearly 150 times more likely than a model with AIC = -120.
Unfortunately, not yet a generic function, so works only with 'drc' objects.
}

\value{
 A value of the first model \code{x} being more likely than the second model \code{y}.
 If large, first model is better. If small, second model is better.
}

\author{
  Andrej-Nikolai Spiess
}

\examples{
### compare two models, four-parameter and five-parameter
m1 <- multdrc(F1.1 ~ Cycles, data = reps, fct = l4())
m2 <- multdrc(F1.1 ~ Cycles, data = reps, fct = l5())
evidence(m2, m1)
### ratio of two AIC's
evidence(-120, -123)
}

\keyword{models}
\keyword{nonlinear}
back to top