plot.profile.etasclass: plot method for profile.etasclass objects (profile likelihood of ETAS model)

Description

plot method for profile.etasclass objects (profile likelihood of ETAS model). Plots a smooth interpolation of the profile likelihood of a parameter of an ETAS model, as output from profile.etasclass.

Usage

## S3 method for class 'profile.etasclass':
plot(x,prob=c(0.9,0.95,0.99),...)

Arguments

An object of the class profile.etasclass.

prob

A vector of coverage probability for the asymptotic confidence interval computed using -2log(LR). Default value prob=c(0.9,0.95,0.99).

...

other arguments.

Value

Plots a profile likelihood (in the scale -2log(LR)), and plots horizontal lines corresponding to the percentiles of a 1df chi-square variable of levels prob; the approximate confidence intervals corresponding to the levels prob are printed. Returns a list:
spline.profileThe spline interpolation of the profile likelihood.
confThe approximate confidence intervals corresponding to the levels prob.
probThe prob values used.

Details

Plots a spline interpolation of the profile likelihood for a parameter of the ETAS model for earthquake seismicity, computed with profile.etasclass; the order of parametrs is: (mu, k0, c, p, a, gamma, d, q). A plot method is defined for profile.etasclass objects. A number of grid points nprofile of 7 (the default) usually is enough to have a good interpolation of the profile likelihood.

Examples

Run this code

## see example in profile.etasclass

Run the code above in your browser using DataLab