# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
ts <- formatTS(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL, set1)
neubuerg <- clean(neubuerg)
# calculate Weibull parameters
neubuerg.wb <- weibull(neubuerg, 1)
# plot
plotWbDir(neubuerg.wb)
# show parameters in legend
plotWbDir(neubuerg.wb, show.ak=TRUE)
# customize plot
plotWbDir(neubuerg.wb, col=c(rainbow(12), gray(0.4)),
lty=c(rep(3, 12), 1), lwd=c(rep(1, 12), 2), cex=1.2)
plotWbDir(neubuerg.wb, lty=c(2,2,3,3,3,4,4,4,5,5,5,2,1),
lwd=rep(1, 12), cex=0.9)
Run the code above in your browser using DataLab