R/plot.tpopt.R

In rodd: Optimal Discriminating Designs

plot.tpopt <- function(x, ...)
{
    res <- x
    curve(Psi(x, res$p, res$eta, res$theta.fix, res$theta.var), res$x.lb, res$x.rb, ylab = "Psi", n = 1001, type = "n")
    grid()
    curve(Psi(x, res$p, res$eta, res$theta.fix, res$theta.var), res$x.lb, res$x.rb, n = 1001, add = TRUE)
    points(res$x, Psi(res$x, res$p, res$eta, res$theta.fix, res$theta.var))
    abline(h = Psi(res$x[1], res$p, res$eta, res$theta.fix, res$theta.var), lty = 2)
}