30.09.2014 Views

Package 'extRemes' - What are R and CRAN?

Package 'extRemes' - What are R and CRAN?

Package 'extRemes' - What are R and CRAN?

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

26 distill.fevd<br />

par(mfrow=c(1,2))<br />

plot(x, devd(x, 1, 1, -0.5), type="l", col="blue", lwd=1.5,<br />

ylab="GEV df")<br />

# Note upper bound at 1 - 1/(-0.5) = 3 in above plot.<br />

lines(x, devd(x, 1, 1, 0), col="lightblue", lwd=1.5)<br />

lines(x, devd(x, 1, 1, 0.5), col="darkblue", lwd=1.5)<br />

legend("topright", legend=c("(reverse) Weibull", "Gumbel", "Frechet"),<br />

col=c("blue", "lightblue", "darkblue"), bty="n", lty=1, lwd=1.5)<br />

plot(x, devd(x, 1, 1, -0.5, 1, type="GP"), type="l", col="blue", lwd=1.5,<br />

ylab="GP df")<br />

lines(x, devd(x, 1, 1, 0, 1, type="GP"), col="lightblue", lwd=1.5)<br />

lines(x, devd(x, 1, 1, 0.5, 1, type="GP"), col="darkblue", lwd=1.5)<br />

legend("topright", legend=c("Beta", "Exponential", "P<strong>are</strong>to"),<br />

col=c("blue", "lightblue", "darkblue"), bty="n", lty=1, lwd=1.5)<br />

# Emphasize the tail differences more by using different scale parameters.<br />

par(mfrow=c(1,2))<br />

plot(x, devd(x, 1, 0.5, -0.5), type="l", col="blue", lwd=1.5,<br />

ylab="GEV df")<br />

lines(x, devd(x, 1, 1, 0), col="lightblue", lwd=1.5)<br />

lines(x, devd(x, 1, 2, 0.5), col="darkblue", lwd=1.5)<br />

legend("topright", legend=c("(reverse) Weibull", "Gumbel", "Frechet"),<br />

col=c("blue", "lightblue", "darkblue"), bty="n", lty=1, lwd=1.5)<br />

plot(x, devd(x, 1, 0.5, -0.5, 1, type="GP"), type="l", col="blue", lwd=1.5,<br />

ylab="GP df")<br />

lines(x, devd(x, 1, 1, 0, 1, type="GP"), col="lightblue", lwd=1.5)<br />

lines(x, devd(x, 1, 2, 0.5, 1, type="GP"), col="darkblue", lwd=1.5)<br />

legend("topright", legend=c("Beta", "Exponential", "P<strong>are</strong>to"),<br />

col=c("blue", "lightblue", "darkblue"), bty="n", lty=1, lwd=1.5)<br />

## End(Not run)<br />

distill.fevd<br />

Distill Parameter Information<br />

Description<br />

Usage<br />

Distill parameter information (<strong>and</strong> possibly other pertinent inforamtion) from fevd objects.<br />

## S3 method for class fevd<br />

distill(x, ...)<br />

## S3 method for class fevd.bayesian<br />

distill(x, cov = TRUE, FUN = "mean", burn.in = 499, ...)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!