aboutsummaryrefslogtreecommitdiff
path: root/R/lod.R
blob: 9f90d48199c88c4bc0eaa7098b4bb82d57296410 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
lod <- function(object, ..., alpha = 0.05, k = 1, n = 1, w = "auto")
{
  UseMethod("lod")
}

loq <- function(object, ..., alpha = 0.05, k = 3, n = 1, w = "auto")
{
  lod(object = object, alpha = alpha, k = k, n = n, w = w)
}

lod.default <- function(object, ..., alpha = 0.05, k = 1, n = 1, w = "auto")
{
  stop("lod is only implemented for univariate lm objects.")
}

lod.lm <- function(object, ..., alpha = 0.05, k = 1, n = 1, w = "auto")
{
  f <- function(x) {
    y <- predict(object, data.frame(x = x))
    p <- inverse.predict(object, rep(y, n), ws = w, alpha = alpha)
    (p[["Prediction"]] - k * p[["Confidence"]])^2
  }
  tmp <- optimize(f,interval=c(0,max(object$model$x)))
  return(tmp$minimum)
}

Contact - Imprint