# NOT RUN {
z <- revd(100, loc=20, scale=0.5, shape=-0.2)
fit <- fevd(z)
fit
strip( fit )
strip( fit, use.names = FALSE )
# Compare with ...
distill( fit )
distill( fit, cov = FALSE )
# }
# NOT RUN {
data( "Fort" )
fit <- fevd(Prec, Fort, threshold=0.395,
scale.fun=~sin(2 * pi * (year - 1900)/365.25) +
cos(2 * pi * (year - 1900)/365.25),
type="PP", method="Bayesian", iter=1999, use.phi=TRUE, verbose=TRUE)
fit
strip( fit )
strip( fit, burn.in = 700 )
strip( fit, FUN = "postmode" )
# }
Run the code above in your browser using DataLab