pp.fitrange: Fitting the Point Process Model Over a Range of Thresholds

Description Usage Arguments Value See Also Examples


Maximum-likelihood fitting for a stationary point process model, over a range of thresholds. Graphs of parameter estimates which aid the selection of a threshold are produced.


pp.fitrange(data, umin, umax, npy = 365, nint = 10, show = FALSE, ...)



A numeric vector of data to be fitted.

umin, umax

The minimum and maximum thresholds at which the model is fitted.


The number of observations per year/block.


The number of fitted models.


Logical; if TRUE, print details of each fit.


Optional arguments to pp.fit.


Three graphs showing maximum likelihood estimates and confidence intervals of the location, scale and shape parameters over a range of thresholds are produced. A list object is returned invisibly with components: 'threshold' numeric vector of length 'nint' giving the thresholds used, 'mle' an 'nint X 3' matrix giving the maximum likelihood parameter estimates (columns are location, scale and shape respectively), 'se' an 'nint X 3' matrix giving the estimated standard errors for the parameter estimates (columns are location, scale and shape, resp.), 'ci.low', 'ci.up' 'nint X 3' matrices giving the lower and upper 95 intervals, resp. (columns same as for 'mle' and 'se').

See Also

pp.fit, mrl.plot, gpd.fit, gpd.fitrange


## Not run: data(rain)
## Not run: pp.fitrange(rain, 10, 40)

Search within the ismev package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.