SLIDE 1
Use R for Climate Research
Research & Teaching Hurricane Climatology with R
James B. Elsner Thomas H. Jagger Florida State University Climatek, Inc Support: U.S. NSF, RPI
Use R for Climate Research Research & Teaching Hurricane - - PowerPoint PPT Presentation
Use R for Climate Research Research & Teaching Hurricane - - PowerPoint PPT Presentation
Use R for Climate Research Research & Teaching Hurricane Climatology with R James B. Elsner Thomas H. Jagger Florida State University Climatek, Inc Support: U.S. NSF, RPI Summary We use existing packages to help answer important
SLIDE 2
SLIDE 3
par(las=1,mgp=c(2,0.5,0),tcl=-0.3,mar=c(5,5,4,2)) x=boxplot(GTC$Wmax~as.factor(GTC$Year),plot=F) boxplot(GTC$Wmax~as.factor(GTC$Year),xaxt="n", ylim=c(10,86),xlab="Year",ylab="Maximum Wind Speed (m/s)",main="") xx=1:26 abline(lm(x$stats[5,]~xx),col="blue",lty=1,lwd=1.5) abline(lm(x$stats[4,]~xx),col="green3",lty=1,lwd=1.5) abline(lm(x$stats[3,]~xx),col="red",lwd=1.5) mtext("a",side=2,padj=-11,line=1.5,cex=1.1) axis(1,at=seq(1,26,5),labels=seq(1981,2006,5)) axis(1,at=seq(1,26,2),labels=table(GTC$Year)[seq(1,26,2)], line=-2,tick=F,cex.axis=0.7) axis(1,at=seq(2,26,2),labels=table(GTC$Year)[seq(2,26,2)], line=-1.5,tick=F,cex.axis=0.7)
SLIDE 4
model.qr=rq(GTC$Wmax~GTC$Year,tau=seq(0.1,0.9,0.1)) par(las=1,mgp=c(2,0.5,0),tcl=-0.3) plot(summary(model.qr,alpha=.05,se="iid"),xaxt="n", ylim=c(-0.1,0.4),mar=c(5,5,4,2)+0.1,parm=2,pch=16, cex=1.1,ylab="Trend (m/s per yr)",xlab="Quantile",main="") axis(3,at=seq(0.2,0.8,0.2),labels=round(quantile(GTC$Wmax,seq (0.2,0.8,0.2)))) mtext("Maximum Wind Speed (m/s)",side=3,line=2) mtext("b",side=2,padj=-11,line=1.5,cex=1.1)
SLIDE 5
Quantile regression of satellite- derived lifetime maximum wind speed
- n year
Global Global
SLIDE 6
Download the HurRisk.zip file from http://garnet.fsu.edu/~jelsner/www/forecasts.html Unzip the file to create a directory called HurRisk For all platforms: setwd("HurRisk") Make sure you have the ismev, maps and foreign packages. Copy/paste the following code require(foreign); require(ismev); require(maps) load("data/best.use.R") load("data/landGrid.R") source("funs/selectData.R") source("funs/ismevplus.R") test0=HurRisk(r=200,distunits="nm",velunits="kt",field="Wmax" ,threshold=64) test1=HurRisk(r=250,distunits="km",velunits="m/s",field="Wma xS",threshold=32)
SLIDE 7
SLIDE 8
SLIDE 9
SLIDE 10
SLIDE 11
Use R for Climate Research Tutorial Use R for Climate Research Video Tutorial HurRisk Model
SLIDE 12