# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
set2 <- createSet(height=30, v.avg=winddata[,6], dir.avg=winddata[,16])
set3 <- createSet(height=20, v.avg=winddata[,10])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set1, set2, set3)
neubuerg <- clean(mast=neubuerg)
# calculate availability
neubuerg.avail <- availability(mast=neubuerg, print=FALSE)
# plot availability
plotAvailability(avail=neubuerg.avail)
plotAvailability(avail=neubuerg.avail, set=2) # one dataset
plotAvailability(avail=neubuerg.avail, set="set2") # same as above
plotAvailability(avail=neubuerg.avail, set=c(1,2)) # several datasets
# customize plot
plotAvailability(avail=neubuerg.avail, border="darkgray", cex.axis=0.7,
cex.lab=0.9, col=c("darkgreen", "orange", "red4"), col.axis="blue",
col.lab="blue", fill=c("lightgreen", "yellow", "red"), lwd=0.5,
mar=c(1,1,1,1), plot.names=FALSE, xlab="jour", ylab="mois")
Run the code above in your browser using DataLab