# 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
plotWeibull(neubuerg.wb)
# show parameters in legend
plotWeibull(neubuerg.wb, show.ak=TRUE)
# customize bars
plotWeibull(neubuerg.wb, col="lightgray", border="darkgray")
plotWeibull(neubuerg.wb, col=3, border="white", cex=1.2)
# customize line
plotWeibull(neubuerg.wb, line="black", lty="dotdash", lwd=2, cex=0.7)
plotWeibull(neubuerg.wb, line="orange", lty=2, lwd=1.5)
Run the code above in your browser using DataLab