# 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)
# if only one of v.set and dir.set is given,
# the dataset is assigned to both
neubuerg.freq.2 <- frequency(mast=neubuerg, v.set=1)
neubuerg.freq.3 <- frequency(mast=neubuerg, dir.set=1)
# compare output
neubuerg.freq
neubuerg.freq.2
neubuerg.freq.3
# use different datasets for wind speed and direction
neubuerg.freq.4 <- frequency(neubuerg, 3) # no direction in dataset
neubuerg.freq.4 <- frequency(neubuerg, v.set=3, dir.set=2)
# compare output
neubuerg.freq
neubuerg.freq.4
# change number of direction sectors
neubuerg.freq.5 <- frequency(neubuerg, v.set=1, num.sectors=4)
neubuerg.freq.6 <- frequency(neubuerg, v.set=1, num.sectors=16)
# compare output
neubuerg.freq.5
neubuerg.freq.6
# calculate frequency for 1 m/s speed bins and without binning
neubuerg.freq.7 <- frequency(neubuerg, v.set=1, bins=1:25)
neubuerg.freq.8 <- frequency(neubuerg, v.set=1, bins=0:25) # same as above
neubuerg.freq.9 <- frequency(neubuerg, v.set=1, bins=NULL)
# compare output
# note: bins exceeding maximum wind speed of the data dropped
neubuerg.freq.7
neubuerg.freq.8
neubuerg.freq.9
Run the code above in your browser using DataLab