# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2])
set30 <- createSet(height=30, v.avg=winddata[,6])
set20 <- createSet(height=20, v.avg=winddata[,10])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set40, set30, set20)
neubuerg <- clean(mast=neubuerg)
# calculate stats
neubuerg.stats <- monthStats(mast=neubuerg, print=FALSE)
# plot all datasets
plotMonthStats(stats=neubuerg.stats)
# plot one dataset
plotMonthStats(stats=neubuerg.stats, set=1)
plotMonthStats(stats=neubuerg.stats, set="set1") # same as above
# customize plot
plotMonthStats(stats=neubuerg.stats, border="darkgray", bty="l", cex.axis=0.7,
cex.lab=0.9, col=c(gray(0.3), gray(0.7)), col.axis="darkgray",
col.box="darkgray", col.lab="darkgray", col.ticks="darkgray", las=0,
legend=FALSE, mgp=c(2, 0.7, 0), plot.names=FALSE, ylim=c(0,7),
ylab="Wind speed [m/s]")
Run the code above in your browser using DataLab