## S3 method for class 'evd':
profile(fitted, which = names(fitted$estimate), conf = 0.999,
mesh = fitted$std.err[which]/4, xmin = rep(-Inf, length(which)),
xmax = rep(Inf, length(which)), convergence = FALSE, method = "BFGS",
control = list(maxit = 500), ...)
"evd"
.conf
or less can be deriwhich
. The values represent the
distance between the points profiled. By default mesh
is
one quarter of the standard errors. If the fitted object does not
which
. Each value represents the theoretical
lower/upper bound of the corresponding parameter.
The arguments are needed only when a parameter has a
lower/upper bound aconvergence
.)optim
. See optim
for
details."profile.evd"
, which is a list with an
element for each parameter being profiled. The elements are
matrices. The first column contains the values of the profiled
parameter. The second column contains profile deviances. The
remaining columns contain the constrained maximum likelihood
estimates for the remaining model parameters.profile2d.evd
, plot.profile.evd
uvdata <- rgev(100, loc = 0.13, scale = 1.1, shape = 0.2)
M1 <- fgev(uvdata)
M1P <- profile(M1)
par(mfrow = c(2,2))
cint <- plot(M1P)
cint
Run the code above in your browser using DataLab