Skip to contents

Draws the decision boundary for the Rayleigh V test against a specified direction mu0. Unlike the Rayleigh test (a circle, see add_critical_r), the V test privileges one direction, so its boundary is a straight line perpendicular to mu0 at distance \(c = z_\alpha / \sqrt{2n}\) from the centre. A mean-direction arrow (add_heading_arrow) is V-significant if and only if its tip falls on the far side of this line – equivalently, if its projection onto mu0 exceeds \(c\).

Usage

add_critical_v_line(
  hd,
  mu0,
  alpha = 0.05,
  angle_col = "heading",
  facets = NULL,
  group_col = NULL,
  show_region = FALSE,
  colour = "firebrick",
  linetype = "dashed",
  linewidth = 0.6,
  region_fill = "firebrick",
  region_alpha = 0.08,
  axial = FALSE,
  n_pts = 100L
)

Arguments

hd

Data frame of headings with a heading column (radians).

mu0

Hypothesised direction in radians (unit-circle convention).

alpha

Significance level. Default 0.05.

angle_col

Heading column name. Default "heading".

facets

Character vector of column names used for faceting (panel placement). Each occupied combination gets its own boundary; the columns are attached to the returned data so facet_wrap/facet_grid routes each boundary to the correct panel. NULL (default) pools across facets.

group_col

Optional single column. Subdivides each facet cell into subgroups, one boundary per subgroup (each with its own n). All boundaries are drawn in the fixed colourgroup_col does not map colour.

show_region

Logical; shade the rejection segment. Default FALSE.

colour

Line colour. Default "firebrick".

linetype

Line type. Default "dashed".

linewidth

Line width. Default 0.6.

region_fill

Fill colour for the rejection region. Default "firebrick".

region_alpha

Fill opacity. Default 0.08.

axial

Logical. Render the V-test boundary for axial (bidirectional, mod-pi) data: the decision chord is mirrored to both poles. Default `FALSE`.

n_pts

Points approximating the rejection arc. Default 100L.

Value

A list of ggplot2 layers, or NULL if no group has a boundary inside the unit circle.

Details

The line is clipped to the unit circle (drawn as a chord). With show_region = TRUE the circular segment beyond the line – the rejection region – is shaded.

Sample size n is taken per boundary from hd: one boundary per occupied combination of facets and group_col columns, or a single pooled boundary when both are NULL.