Skip to contents

Bias Rose

Usage

rose_metbias(
  data,
  ws = "ws",
  wd = "wd",
  ws2 = "ws2",
  wd2 = "wd2",
  angle = 10,
  facet = NULL,
  bias_corr = TRUE,
  breaks = NA,
  width = 1,
  border_colour = NA,
  alpha = 1
)

Arguments

data

A data frame containing fields ws, wd, ws2 and wd2.

ws

Name of the column representing wind speed.

wd

Name of the column representing wind direction.

ws2

Name of the column representing wind speed (2).

wd2

Name of the column representing wind direction (2).

angle

Angle of the spokes. Ideally a number by which 360 is evenly divisible (see bias_corr).

facet

One or two faceting columns. facet determines how the data are split and then plotted. When facet is length 1 it is passed to ggplot2::facet_wrap(), and when it is length 2 it is passed to ggplot2::facet_grid() with the first element being used as columns and the second rows. Some other options (e.g., multiple pollutant columns) can limit the the number of faceting columns to 1.

bias_corr

When angle does not divide exactly into 360 a bias is introduced in the frequencies when the wind direction is already supplied rounded to the nearest 10 degrees, as is often the case. For example, if angle = 22.5, N, E, S, W will include 3 wind sectors and all other angles will be two. A bias correction can made to correct for this problem. A simple method according to Applequist (2012) is used to adjust the frequencies.

breaks

The number of break points to use when binning the data (used alongside ws_int), or a vector of specific break points. For example, breaks = 6 and ws_int = 2 generates the break points 2, 4, 6 & 8 m/s, whereas breaks = c(0, 1, 10, 100) breaks the data into four segments: <1, 1-10, 10-100, & >100.

width

Width of each bar as a fraction of its maximum width. width = 1 makes all bars meet at their edges and width = 0 makes them disappear entirely. Defaults to 0.9.

border_colour

Border colour for shaded areas. Default is no border.

alpha

The transparency of the plot. This is mainly useful to overlay it on a map.

See also

Other polar directional analysis functions: polar_annulus(), polar_cluster(), polar_diff(), polar_freq(), polar_plot(), rose_percentile(), rose_pollution(), rose_wind()