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(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL,
set40, set30, set20)
neubuerg <- clean(neubuerg)
# calculate stats
neubuerg.stats <- monthStats(neubuerg)
# plot all datasets
plotMonthStats(neubuerg.stats)
# plot one dataset
plotMonthStats(neubuerg.stats, set=1)
# change colour, labels and text size
plotMonthStats(neubuerg.stats, col=c(gray(0.4), gray(0.6)), cex=1.2,
ylab=list(rep("Wind speed [m/s]", 3), c("40 m", "30 m", "20 m")))
plotMonthStats(neubuerg.stats, set=1, col=c(2,3), cex=0.8,
ylab=list("C1.A1", "Velocity [m/s]"))
Run the code above in your browser using DataLab