R/plotKinSpecEst.R

"plotKinSpecEst" <-
  function(t, plotoptions, multimodel) {
    oldpar <- par(no.readonly = TRUE)
    on.exit(par(oldpar))
    if (dev.cur() != 1) {
      dev.new()
    }
    resultlist <- multimodel@fit@resultlist
    par(mfrow = c(1, 2), mar = c(5, 2, 2, 2))
    plotClp(multimodel, t, plotoptions,
      newplot = FALSE,
      kinspecerr = plotoptions@kinspecerr
    )
    ## plotEst(multimodel, plotoptions)
    if (length(plotoptions@title) != 0) {
      mtext(plotoptions@title,
        side = 3, outer = TRUE,
        line = 1
      )
      par(las = 2)
    }
    if (dev.interactive() && length(plotoptions@makeps) != 0) {
      if (plotoptions@output == "pdf") {
        pdev <- pdf
      } else {
        pdev <- postscript
      }

      dev.print(device = pdev, file = paste(plotoptions@makeps,
        "_kinspecrates.", plotoptions@output,
        sep = ""
      ))
    }
  }

Try the TIMP package in your browser

Any scripts or data that you put into this service are public.

TIMP documentation built on Dec. 28, 2022, 3:06 a.m.