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 questions about climate change (quantreg). We develop and maintain our own scripts from bits and pieces of other packages for our students and our clients (HurRisk). We develop teaching and research resources for our students and our clients (Use R for Climate Research).

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)

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)

Global Global Quantile regression of satellite- derived lifetime maximum wind speed on year

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)

Use R for Climate Research Tutorial Use R for Climate Research Video Tutorial HurRisk Model

More Information Google hurricane climate http://garnet.fsu.edu/~jelsner/www jelsner@fsu.edu tjagger@blarg.net