# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
set30 <- createSet(height=30, v.avg=winddata[,6], dir.avg=winddata[,16])
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 frequency
neubuerg.freq <- frequency(neubuerg, 1)
neubuerg.freq.2 <- frequency(neubuerg, 1, bins=NULL)
# plot frequency
plotFrequency(neubuerg.freq)
plotFrequency(neubuerg.freq.2)
# customize plot
plotFrequency(neubuerg.freq, col=rainbow(5))
plotFrequency(neubuerg.freq.2, col=gray(0.4), cex=1.5)
Run the code above in your browser using DataLab