# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set1)
neubuerg <- clean(mast=neubuerg)
# calculate Weibull parameters
neubuerg.wb <- weibull(mast=neubuerg, v.set=1, print=FALSE)
# plot
plotWeibull(wb=neubuerg.wb)
# show parameters in legend
plotWeibull(wb=neubuerg.wb, show.ak=TRUE)
# customize plot
plotWeibull(wb=neubuerg.wb, bty="l", bty.leg="l", cex.axis=1.2,
cex.lab=1.4, cex.leg=0.9, col.axis="darkgray",
col.box="darkgray", col.lab="darkgray", col.leg="darkgray",
col.ticks="darkgray", las=0, leg.text=c("measured", "calculated"),
mar=c(3,3,0.5,0.5), mgp=c(1.8,0.5,0), pos.leg="right",
xlab="velocity [m/s]", ylab="frequency [%]", xlim=c(0,20),
ylim=c(0,15), x.intersp=1, y.intersp=1)
# customize bars
plotWeibull(wb=neubuerg.wb, border="darkgray", breaks=seq(0,21,0.5),
col="lightgray")
# customize line
plotWeibull(wb=neubuerg.wb, line="black", lty="dotdash", lwd=2)
Run the code above in your browser using DataLab