fitPS (version 1.0.1)

plot.psFit: S3 plot method for an object of class psFit

Description

S3 plot method for an object of class psFit

Usage

# S3 method for psFit
plot(
  x,
  ylim = c(0, 1),
  conf = FALSE,
  conf.level = 0.95,
  ci.type = c("wald", "prof"),
  log.scale = FALSE,
  ...
)

Value

No return value, called for side effects

Arguments

x

an object of class psFit, usually from fitDist or fitZIDist.

ylim

the limits of the y-axis.

conf

if TRUE, and the model is the the zeta model (as opposed to the zero-inflated zeta (ZIZ), then confidence intervals (based on the standard error of the shape parameter) are drawn on the plot. If the ZIZ model has been used, then this is ignored.

conf.level

the confidence level for the confidence intervals. Must be between 0.75 and 0.99.

ci.type

Specifies the type of confidence interval. If conf == TRUE, then then ci.type can be either "wald" "prof" (or an abbreviation), depending on whether the Wald interval or the profile likelihood interval should be used. Note that these are intervals on the shape parameter and not the density heights. Therefore the intervals around the probabilities should not really be thought of as confidence intervals but rather something more similar to a "sensitivity" interval.

log.scale

if TRUE the \(y\)-axis is changed to a logarithmic (base 10) axis.

...

other arguments passed to plot.

Examples

p = readData(system.file("extdata", "p.xlsx", package = "fitPS"))
fit = fitDist(p)
plot(fit)

## An example with Wald generated intervals
plot(fit, conf = TRUE)

plot(fit, conf = TRUE, ci.type = "p")