uvdata <- rgev(100, loc = 0.13, scale = 1.1, shape = 0.2)
M1 <- fgev(uvdata)
if (FALSE) M1P <- profile(M1)
if (FALSE) par(mfrow = c(2,2))
if (FALSE) cint <- plot(M1P, ci = c(0.95, 0.99))
if (FALSE) cint
Run the code above in your browser using DataLab